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Abstract 

This  thesis  addresses  the  application  of  optimal,  multiobjective  control  theory  control 
theory  to  flight  control  design  for  the  approach  and  landing  phase  of  flight.  Five  flight  control 
systems  were  designed  using  classical,  H2,  Hx,  and  Mixed  H2/H^  methods.  The  MATLAB™ 
MUTOOLS™  and  AFIT  MXTOOLS  toolboxes  were  used  to  produce  the  optimal,  multiobjective 
designs.  These  designs  were  implemented  for  flight  test  on  the  Calspan  VSS  I  Learjet, 
simulating  the  unstable  longitudinal  dynamics  of  an  F-16  type  aircraft.  A  limited  handling 
qualities  investigation  was  performed.  Model  following  was  used  in  the  design  phase  to  meet 
handling  qualities  specifications.  The  designs  were  successfully  implemented  and  verified  on 
the  Calspan  Learjet  prior  to  flight  test.  An  unmodeled  aircraft  mode  was  discovered  just  prior  to 
flight  test  that  made  three  of  the  designs  slightly  unstable.  However,  all  of  the  designs  achieved 
Level  II  or  better  Cooper-Harper  handling  qualities  ratings  for  the  landing  tasks  performed 
illustrating  that  the  optimal  multiobjective  methods  used  can  give  acceptable  or  better  handling 
qualities. 
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Flight  Test  and  Handling  Qualities  Analysis  of  a  Longitudinal  Flight 
Control  System  Using  Multiobjective  Control  Design  Techniques 


1.  Introduction 


1. 1  Motivation/Background 

All  aircraft  control  system  designs  are  based  on  a  mathematical  model  of  the  aircraft  at  a 
specific  flight  condition  (altitude,  airspeed,  weight,  load  factor,  and  gear/flap  configuration). 
Usually,  this  model  is  not  perfect.  There  will,  in  general,  be  many  underlying  assumptions  made 
to  develop  the  aircraft  model.  These  assumptions  include  such  things  as:  approximation  of 
stability  derivatives  not  easily  measured,  linearization  of  the  aircraft  equations  of  motion,  and 
various  other  simplifications  that  make  the  control  design  problem  tractable.  The  end  result  is  a 
model  which  may  be  far  from  reality  in  its  representation  of  the  aircraft. 

Classical  control  design  techniques,  such  as  root  locus,  are  limited  in  the  sense  that  they 
offer  no  way  to  incorporate  allowances  for  the  assumptions  made  about  the  aircraft  model.  If  the 
actual  aircraft  model  is  significantly  different  from  the  design  model,  the  designer  is  forced  to  go 
back  and  adjust  control  gains  or  possibly  redesign  the  entire  control  law  if  the  response  is  not  as 
desired.  This  tuning  process  becomes  cumbersome  and  time  consuming  for  any  sort  of  complex, 
multiple  loop  control  scheme.  Even  at  a  specific  flight  condition  where  a  good  aircraft  model 
exists,  there  are  in-flight  factors  such  as  wind  changes,  pressure  changes,  temperature  changes, 
and  weight  change  due  to  fuel  consumption  which  will  alter  the  aircraft  model.  Finally,  the 
controller  will  receive  measurements  that  are  corrupted  by  different  types  and  levels  of  noise. 

All  of  these  things  can  be  looked  at  as  uncertainties  that  affect  the  aircraft  model.  If  there  was  a 
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way  of  incorporating  these  aircraft  model,  environmental,  and  noise  uncertainties  into  the 
problem  setup,  it  should  be  possible  to  design  a  more  capable  controller. 

Multiobjective  control  design  techniques  offer  the  promise  of  addressing  these 
uncertainties  that  root  locus  techniques  cannot.  These  techniques  allow  the  designer  to 
incorporate  different  types  of  weightings  to  account  for  things  like  wind  disturbances,  sensor 
noise  distortion,  and  unmodeled  or  changing  aircraft  dynamics  into  the  design.  Therefore 
multiobjective  controllers  should  produce  better  performance  for  a  larger  range  of  flight 
conditions  (i.e.  they  are  more  robust)  than  the  classical  controller  . 

Two  classes  of  multiobjective  control  techniques  that  allow  the  designer  to  incorporate 
the  types  of  weightings  previously  discussed  are  //j  optimal  control,  and  optimal  control. 
H2  optimal  control  was  derived  to  design  for  noise  rejection.  H2  optimization  minimizes  the 
energy  of  a  system's  response  with  respect  to  exogenous  noise  disturbances  that  are  modeled  as 
zero-mean,  unit  intensity  white  Gaussian  noises,  by  minimizing  the  specified  transfer  function 
two-norm.  For  example,  an  H2  control  design  in  an  aircraft  might  minimize  the  noise  on  a  pitch 
rate  (q)  signal  due  to  turbulence  and  pitch  rate  sensor  noise  inputs.  optimization  minimizes 
the  energy  output  of  a  system  due  to  an  unknown,  but  bounded  energy,  input  by  minimizing  the 
specified  transfer  function  infinity-norm.  An  example  here  would  be  to  include  a  weight  in  the 
design  to  account  for  unmodeled  aircraft  dynamics  but  still  require  tight  tracking  of  the  control 
command.  There  is  a  third  class  of  optimal  control  called  optimal  control  that  was  considered 
in  this  research.  It  is  closely  related  to  optimal  control.  The  difference  between and 
optimization  is  that  fj.  optimal  control  handles  multiple  uncertainty  constraints  better  than 
methods.  Designs  that  included  fj.  optimal  control  techniques  did  not  make  it  to  Phase  II  of 
this  research  for  reasons  that  will  be  apparent  in  Chapter  IV. 


1-2 


In  a  realistic  aircraft  control  design  problem,  there  are  both  stochastic  noise  sources  and 
uncertainties  related  to  unmodeled  aircraft  dynamics.  Therefore,  there  might  be  an  advantage  to 
designing  a  controller  that  leverages  the  benefits  of  both  H2  and  optimization. 

7.2  The  Mixed-Norm  Problem 

One  of  the  earliest  formulations  of  a  mixed  H2  /  problem  was  derived  by  Bernstein 
and  Haddad  [1]  in  1989.  Their  approach  minimized  an  overbound  to  the  i/j 
transfer  function  and  also  satisfied  an  constraint  for  a  related  transfer  function.  The  inputs 
to  each  transfer  function  were  required  to  be  the  same  while  the  outputs  could  be  different.  In 
other  words,  there  existed  a  controller  that  provided  a  tradeoff  between  robustness  and  noise 
rejection.  Ridgely,  et  al,  [2,  3]  showed  in  1992  that  the  general  mixed  H2  /  problem,  where 
no  overbound  to  the  H2  norm  is  used  and  the  subject  transfer  function  need  not  be  related,  could 
not  be  solved  analytically.  They  also  proposed  a  numerical  solution  to  the  problem  for  a  fixed 
order  controller.  This  solution  required  unique  solutions  to  Lyapunov  equations  and  stabilizing 
solutions  to  Riccati  equations.  Although  Ridgely,  et  al,  developed  a  numerical  solution  to  the 
problem,  the  implementation  still  did  not  lend  itself  to  efficient  controller  design  because  of  the 
large  number  of  iterations  required  to  solve  the  Lyapunov  and  Riccati  equations.  Walker  [4] 
improved  on  Ridgely's  method  in  1994  .  He  re-cast  Ridgely's  original  problem  formulation  as  a 
convex  optimization  problem.  However,  this  formulation  did  not  address  the  problem  of  actually 
finding  the  optimal  solution.  Further  work  by  Megretski  [5]  in  this  area  indicates  that  the 
optimal  solution  is  not  rational.  It  can  only  be  expressed  as  an  infinite  dimensional  compensator. 
Thus,  using  the  actual  optimal  compensator  would  not  be  practical  from  an  applications 
standpoint. 
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For  this  reason.  Walker  reformulated  the  problem  as  a  Lagrange  multiplier  problem  so 
that  it  could  be  solved  numerically.  The  ultimate  result  of  the  work  begun  by  Walker  was  the 
creation  of  MXTOOLS  for  MATLAB™  by  Jacques,  Canfield,  and  Ridge ly  [6],  This  is  a  public 
domain  toolbox  that  efficiently  solves  the  mixed  norm  control  problem  for  a  fixed  order 
controller.  The  toolbox  is  geared  toward  quick  design  of  a  mixed  norm  controller,  where  an  H2 
objective  function  is  constrained  by  an  ,  // ,  or  constraint  (  £,  constraints  will  not  be 
considered  in  this  thesis).  The  program  also  allows  for  multiple  constraints.  The  user  is  only 
required  to  provide  the  plant  model  and  the  constraint  model(s)  to  be  included  in  the  problem 
formulation.  The  toolbox  then  forms  the  problem,  checks  whether  or  not  any  constraints  are 
active,  and  iterates  until  a  solution  is  found.  The  program  uses  a  Sequential  Quadratic 
Programming  (SQP)  algorithm  to  find  the  optimal  solution.  An  advantage  of  this  type  of 
algorithm  is  that  it  is  efficient  at  solving  constrained  optimization  problems. 

Edwards  [7]  was  the  first  to  use  MXTOOLS  to  create  a  mixed  Hj  /  //„  control  design 
that  was  flight  tested.  He  also  flight  tested  three  other  designs:  an  optimal  control  design,  an 
optimal  control  design,  and  a  classic  control  law  designed  using  root-locus  methods.  All  of 
the  designs  used  a  4th  order,  unstable  F-16  short  period  and  phugoid  aircraft  model  provided  by 
Calspan  SRL  Corporation.  The  mixed  Hj  /  control  was  a  combination  (via  MXTOOLS)  of 
the  flight  test  //j  and  optimal  control  designs.  The  optimal  control  design  included 
weights  for  turbulence  rejection  and  command  tracking.  The  optimal  control  design 
included  weights  to  account  for  model  uncertainty  through  maximizing  stability  margins  and 
tracking  a  command  input.  Because  of  the  setup  of  the  H2  and  designs,  the  mixed  H2  / 
design  included  good  command  tracking,  stability  margins,  and  turbulence  rejection.  The  classic 
control  law  was  designed  using  angle  of  attack  (a)  and  pitch  rate  (q)  feedback  gains. 
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Handling  qualities  predictions  for  all  of  these  designs  were  checked  using  Hoh’s 
Bandwidth  Criteria  for  the  landing  phase  of  flight.  This  criteria  is  defined  in  MIL-STD-1 797A, 
Flying  Qualities  of  Piloted  Aircraft  [8].  All  control  designs  are  desired  to  have  Level  1  handling 
qualities  predictions  based  on  the  Cooper-Harper  aircraft  handling  qualities  rating  scale.  This 
rating  scale  is  presented  in  Appendix  A.  All  of  the  designs  used  by  Edwards  in  flight  test  had 
Level  1  handling  qualities  predictions. 

Edwards’  flight  test  program  was  named  HAVE  INFINITY  by  the  USAF  Test  Pilot 
School.  Each  HAVE  INFINITY  flight  control  system  was  evaluated  by  four  pilots  during  the 
approach  and  landing  phase  of  flight  in  the  Calspan  Variable  Stability  Learjet  24.  Several 
straight-in  and  horizontal  offset  landings  were  made  by  each  pilot  to  evaluate  all  four  flight 
control  systems.  The  horizontal  offset  landings  were  performed  primarily  to  raise  the  pilot’s 
gains,  which  would  expose  any  potential  handling  qualities  “cliffs”  in  the  flight  control  systems. 
HAVE  INFINITY  used  all  of  the  handling  qualities  ratings  for  each  type  of  landing  to  come  up 
with  an  overall  level  rating  for  each  flight  control  design.  The  evaluation  pilots  also  gave  a 
Pilot-In-the-Loop-Oscillation  (PIO)  rating  using  the  scale  defined  in  [8].  The  PIO  rating  gives  a 
measure  of  the  flight  control  system’s  susceptibility  to  PIO  for  a  given  task.  This  rating  scale 
along  with  an  explanatory  decision  tree  are  also  presented  in  Appendix  A.  All  PIO  ratings  were 
between  2  and  5.  The  classical  controller  was  rated  Level  II  on  the  Cooper-Harper  (C-H)  Rating 
Scale,  and  the  H2,  Hco,  and  mixed  H2/H^  designs  all  received  Level  III  C-H  ratings.  The  Level 
III  C-H  ratings  were  primarily  attributed  to  an  uncommanded  pitch  up  in  the  landing  flare.  The 
HAVE  INFINITY  test  team  concluded  that  the  pitch  up  problem  may  have  been  related  to 
problems  in  the  implementation  of  the  HAVE  INFINITY  optimal  control  laws.  It  turned  out  that 
the  H2,  Hoo,  and  mixed  H2/Hao  designs  contained  high  frequency  and  unstable  modes  internal  to 
the  control  laws.  This  made  the  designs  very  difficult  to  implement  on  the  Calspan  Learjet. 
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Subsequent  model  verification  and  validation  testing  revealed  that  predicted  time  and  frequency 
responses  did  not  match  the  actual  time  and  frequency  responses.  The  HAVE  INFINITY  test 
team  felt  that  the  dubious  nature  of  the  implementation  and  the  potential  negative  impact  on 
handling  qualities  rendered  the  results  inconclusive  as  to  the  potential  benefits  of  the 
multiobjective  design  techniques  used. 

Due  to  the  HAVE  fNEINITY  test  results,  the  decision  was  made  to  continue  the 
multiobjective  flight  control  research  in  hopes  of  producing  a  definitive  answer  about  their  utility 
with  regard  to  handling  qualities.  This  decision  was  in  lieu  of  testing  some  of  the  potential 
advantages  of  the  methods  such  as  turbulence  rejection  and  stability  robustness,  as  had  been 
originally  planned.  It  was  felt  that  getting  good  handling  qualities  using  the  optimal  methods  had 
to  be  accomplished  before  more  advanced  testing  could  take  place. 

1.3  Research  Objectives 

The  overall  objective  of  this  thesis  was  to  perform  an  inflight  handling  qualities 
evaluation  of  several  multiobjective,  optimal  flight  control  designs.  Specifically,  this  research 
focused  on  the  development,  simulation,  and  flight  test  of  several  //j  mixed-norm  control 
designs  during  the  approach  and  landing  phase  of  flight.  The  designs  were  flight  tested  on  the 
Calspan  Variable  Stability  Learjet  simulating  the  F-16  handling  characteristics.  All  designs  were 
evaluated  using  Hoh's  Bandwidth  Criteria  for  predicting  pilot  opinion  of  aircraft  handling 
qualities.  This  research  was  conducted  in  two  phases.  Phase  I  was  conducted  at  AFIT.  Phase  I 
assumed  the  F-16  longitudinal  dynamics  could  be  accurately  approximated  on  the  Calspan 
Learjet  by  a  second  order  short  period  transfer  function.  Phase  I  specific  objectives  were  to 
formulate,  synthesize,  and  simulate: 
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1 .  The  model-following  subproblem  to  determine  the  optimal  H2  controller  that 
provided  the  best  possible  noise  rejection  while  following  the  ideal  model. 

2.  The  model-following  subproblem  to  determine  the  optimal  controller  that 
provided  good  command  tracking  in  addition  to  stability  robustness. 

3.  The  ju  subproblem  to  determine  the  optimal  controller  to  handle  multiple 
uncertainties. 

4.  The  mixed  H2  with  multiple  constraints  problem(s)  to  tradeoff  noise 
rejection  with  separate  tracking,  input  model  uncertainty,  and  output  stability 

margin  constraints. 

5.  The  mixed  H2  !  fJ^  problem  to  tradeoff  noise  rejection  and  robust 
performance  for  multiple  uncertainties. 

6.  Determine  the  best  designs  to  use  in  phase  II  based  on  analytical  analysis  and 
simulation. 

Phase  II  specific  objectives  of  this  research  were; 

1 .  Perform  flight  test  evaluation  of  selected  control  designs  from  Phase  I  that 
satisfies  Objective  6. 

2.  Obtain  and  evaluate  qualitative  and  quantitative  pilot  opinion  and  Cooper- 
Harper  Pilot  Ratings  of  the  flight  control  designs  selected  for  flight  test  evaluation. 

3.  Compare  flight  test  handling  qualities  ratings  to  those  predicted  in  Phase  I. 

4.  Draw  conclusions  with  regards  to  apparent  trends  in  Objective  3,  Phase  II. 

5.  Collect  and  archive  data  for  future  use  by  AFIT  and  USAF  TPS. 

1.4  Thesis  Organization 

Chapters  Il-IV  cover  the  analytical  phase  (Phase  I)  of  this  research  and  chapters  V  and 
VI  cover  the  flight  test  phase  (Phase  II).  Chapter  VII  is  a  summary  chapter.  Specifically, 
Chapter  II  contains  the  theory  for and  mixed  H2I  control  design.  This  chapter  II 
also  describes  the  complex  structured  singular  value,  fj. ,  in  terms  of  analysis,  synthesis,  and 
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mixed  !  \3i  ■  Chapter  III  contains  the  setup  and  synthesis  for  the  various  mixed  norm  control 
designs.  The  design  examples  use  a  short  period,  longitudinal  state  space  model  for  the  Calspan 
Variable  Stability  Learjet  simulating  an  F-16  in  the  landing  phase.  The  results  of  the  design 
synthesis  and  simulation  accomplished  during  Phase  I  are  then  presented  in  Chapter  IV.  Flight 
control  design  implementation,  flight  test  techniques  used,  flight  test  structure,  and  model 
verification/validation  are  described  in  Chapter  V.  Handling  qualities  results  for  each  flight 
control  design  flight  tested  are  presented  in  Chapter  VI.  Finally,  Chapter  VII  brings  the 
analytical  and  flight  test  phases  together  and  makes  conclusions  and  recommendations  for 
further  research. 
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II.  Theoretical  Preliminaries 


This  chapter  first  focuses  on  the  single  norm  optimization  methods  that  will  be  used: 
H2,H„,  cind  fj. .  The  last  portion  of  the  chapter  provides  a  brief  overview  of  fixed-order  mixed 
H2  /  control  optimization  and  fixed-order  mixed  H2  !  It  control  optimization. 


2. 1  H2  Optimal  Control 

2.1.1  The  H2  space. 

is  defined  as  the  space  of  all  transfer  function  matrices  which  are  stable  (all 
eigenvalues  in  the  open-left  half  complex  plane)  and  have  a  bounded  two-norm.  The  two- 
norm,  a ,  is  defined  as: 


a  = 


00  ^  K 

\tr[T^^{j(o)T^„{jco)]d(o 

-1  J 


(2.1) 


The  subspace  iR/Zj  defined  as  the  space  of  rational  functions  with  real  coefficients  in  H2 . 
Equation  2.1  is  not  easy  to  compute.  Fortunately,  there  is  an  easier  method  for  computing  the 
two-norm  of  a  transfer  function  (see  [9]  for  a  more  detailed  explanation).  Consider  the  transfer 
function 


G{s)  = 


'A 

B' 

_C 

0_ 

(2.2) 


where  the  above  notation  represents  G{s)  =  C(sl  -A)  ^  B  +  0.  Assuming  the  real  part  of  all 
the  eigenvalues  of  A  are  less  than  zero,  we  can  find  the  controllability  and  observability 
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grammians,  and  ,  which  are  the  positive  semidefmite  solutions  to  the  Lyapunov  equations 

AL^  +  L^A^  +  BB^  =0  (2.3) 

L„A  +  A^L„+C^C  =  0  (2.4) 

From  here  we  can  easily  compute  the  two-norm,  which  is  given  by 

||G(s)||“  =  ft'CCi.C’')  =  (2.5) 

2.J.2  Optimization. 

The  H2  optimal  control  problem  is  described  by  the  block  diagram  in  Figure  2.1 .  This 
notation  is  known  as  a  lower  Linear  Fractional  Transformation  (LFT)  of  P  and  K,  denoted 
Ff{P,K).  The  exogenous  input  to  the  system,  w,  is  zero-mean,  unit  intensity,  white  Gaussian 
noise.  Typically,  these  noise  inputs  would  model  wind  gust  disturbances  and  sensor  noise  in  an 
aircraft.  The  output,  z,  is  whatever  the  control  designer  wishes  to  have  minimized  with  respect  to 
w.  The  plant  model  is  designated  P  and  can  include  weights  on  parameters  that  the  designer 
wishes  to  emphasize  or  de-emphasize  in  certain  frequency  ranges.  An  example  would  be  a 
weight  on  control  use  to  account  for  actual  aircraft  limits  of  control  deflection.  The  goal  then  is 
to  determine  the  feedback  controller,  K,  using  the  method  described  in  this  section.  The 
measurements,  y,  are  the  outputs  of  P  that  are  fed  back  to  K.  The  control  law,  u,  is  the  output  of 
K  that  is  fed  into  P. 
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Figure  2.1  Block  Diagram  for  H2  Problem 

The  objective  of  //j  control  optimization  is  to  find  a  stabilizing  controller  that 
minimizes  the  rms  energy,  or  two-norm,  of  z,  given  the  noise  inputs,  w,  described  previously. 

H2  optimization  is  performed  by  solving  two  Algebraic  Riccati  Equations  (ARE's)  and  is  a 
generalization  of  the  standard  Linear  Quadratic  Gaussian  (LQG)  problem.  This  generalization 
allows  noise  inputs  and  desired  outputs  to  be  placed  anywhere  in  the  problem  formulation. 
Dynamic  weightings  of  outputs  and  inputs  are  also  allowed.  See  Section  B.l,  Appendix  B  for 
further  mathematical  development. 

2.2  Optimal  Control 

2.2.1  The  //„  Space. 

H„  is  defined  as  the  space  of  all  transfer  function  matrices  which  are  stable  (all 
eigenvalues  in  the  open-left  half  complex  plane)  and  have  a  bounded  infinity-norm.  The  infinity- 
norm  can  be  looked  at  as  the  maximum  possible  energy  gain  of  a  system.  Thus,  to  minimize  the 
energy  to  a  deterministic,  but  unknown,  bounded  energy  input  (e.g.  a  pulse),  we  must  minimize 
the  infinity-norm  of  the  associated  transfer  function.  Analogous  to  the  Hj  space,  9^/4,  is  the 
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subspace  of  all  real,  rational  functions.  The  infinity-norm  is  an  induced  norm  and  is  defined 
as 


_ 8 

II 

(2.6) 

IIML 

=  sup"„  , 

(2.7) 

=  supllelL 

(2.8) 

=  sup(T[7;^C7'ty)] 

(2.9) 

(0 


where  2^^  is  a  closed  loop  transfer  function  whose  meaning  will  be  clarified  later. 

An  important  characteristic  of  the  infinity-norm  is  that  since  it  is  an  induced  operator  norm,  it 
has  a  sub-multiplicative  property  [10];  given  F,G  g  then 

MLsWLIIgil  (2.10) 

The  easiest  way  to  calculate  the  infinity-norm  of  a  transfer  function  is  to  calculate  its  maximum 
singular  values  over  a  sufficiently  large  range  of  frequencies  and  select  the  maximum  value. 
However,  this  approach  may  not  be  numerically  practical  since  the  frequency  at  which  the 
maximum  singular  value  will  be  attained  is  not  known  a  priori.  A  more  precise  approach  is 
based  on  the  eigenvalues  of  a  Hamiltonian  matrix  associated  with  a  state  space  realization  of  a 
proper  stable  transfer  function  [10]. 

Consider  the  transfer  function 
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(2.11) 


G(s)  = 


~A 

B' 

D 

The  associated  Hamiltonian  is 


A  +  BR-'D^C 

-C''\l  +  DD^y^C  -(A  +  BR-'D'^Cf 


where  R:-  y^I  —  D^D .  The  infinity-norm  of  the  transfer  function  is  the  smallest  value  of  y 
such  that  JT has  no  eigenvalues  on  the  imaginary  axis. 

2.2.2  Optimization. 

Consider  the  block  diagram  in  Figure  2.2.  This  is  the  same  setup  as  Figure  2.2  except 
that  the  inputs  and  outputs  have  been  defined  differently.  The  input,  d,  is  an  unknown  but 
deterministic  bounded  energy  signal.  The  task  is  then  to  find  the  controller,  K,  that  minimizes 
the  output,  e,  with  respect  to  d. 


Figure  2.2  Block  Diagram  for  the  Ha,  Problem 
The  plant  P  can  be  partitioned  as 


(2.13) 
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such  that 


e  =■  P  jd  +  P  u 

^  ^  ed^  eu 


(2.14) 


(2.15) 

This  is  equivalent  to  minimizing  the  infinity-norm  of  the  closed  loop  transfer  function,  Tg(ji. 
Further  mathematical  development  on  control  optimization  can  be  found  in  Section  B.2, 
Appendix  B. 


2.3  The  Complex  Structured  Singular  Value 

The  next  optimal  control  tool  of  interest  is  the  structured  singular  value,  /a  .  This 
section  will  introduce  /a  and  how  to  find  an  upper  bound  on  it  through  optimization 
techniques.  The  application  of  /i  -synthesis  to  the  F-16  landing  problem  as  a  means  of  ensuring 
robust  stability  and  robust  performance  is  presented  in  Chapter  III.  The  reader  is  invited  to  see 
[11]  for  a  tutorial  on  the  complex  structured  singular  value.  Consider  the  problem  represented 
in  Figure  2.3. 


i  A 


Figure  2.3  System  with  Multiple  Uncertainties 
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In  this  figure,  dr  is  the  reference  input,  Cp  is  the  controlled  output,  Act  is  the  actuator,  G  is  the 
core  plant,  and  K  is  the  controller.  The  system  has  an  uncertainty  at  the  output  of  the  plant. 
These  uncertainties  include  frequency  weights  for  such  things  as  neglected  high  frequency 
dynamics  or  uncertain  parameters  in  the  system.  The  structure  in  the  uncertainty  is  apparent 
from  the  figure.  There  is  an  uncertainty  for  the  actuator  and  another  for  the  plant  output  and 
these  uncertainties  occur  at  different  places  in  the  system,  ju  allows  us  to  use  this  structure  to 
combine  the  two  uncertainties  into  a  single  structured  uncertainty,  and  come  up  with  a  robust 
design  that  is  less  conservative  than  the  standard  design  would  be  in  this  case.  The 
mathematical  development  for  the  structured  singular  value  is  found  in  Section  B.3,  Appendix  B. 

2.4  Mixed  H2  /  Control 

2.4.1  Background. 

Recall  that  both  H2  and  optimization  minimize  the  energy  of  the  chosen  output(s). 
However,  the  input(s)  for  H2  are  characterized  by  white  Gaussian  noise  (wind  gusts  and  sensor 
noise,  for  example),  and  the  inputs  for  are  unknown,  but  bounded  energy  deterministic  inputs 
(unmodelled  dynamics,  etc.).  Since  each  type  of  control  optimization  is  designed  for  a  different 
type  of  input,  there  may  be  some  advantage  to  mixing  the  two  in  a  problem  where  both  types  of 
inputs  are  present.  However,  we  must  also  realize  that  mixing  these  two  different  types  of 
control  designs  will  often  produce  competing  objectives.  If  our  H2  design  is  optimized  for  noise 
rejection  and  our  design  is  optimized  to  provide  good  tracking  by  minimizing  output 
sensitivity,  we  can  expect  to  tradeoff  noise  rejection  performance  for  tracking  performance. 
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Ridgely  [2,3]  and  Walker  [4]  have  both  developed  techniques  to  handle  mixed  norm 
control  optimization  using  the  Hi  subproblem  as  the  objective  function  to  be  minimized,  and 
append  one  or  more  subproblem(s)  as  constraints.  Note  that  we  are  not  doing  control 
optimization  here  but  rather  constraining  the  infinity-norm  to  some  level  specified  by  the 
designer.  This  process  can  be  formally  stated  as  follows 


H^zwIL’  subject  to  \\Tjl<r 

-^adm 


(2.16) 


The  details  of  the  mathematical  development  of  Walker’s  mixed  control  method 

is  found  in  Section  B.4,  Appendix  B.  The  results  of  Walker’s  work  provided  the  basis  for  a 
reliable  numerical  method  for  mixed  control  optimization. 

2.4.2  Numerical  Approach. 

The  method  developed  by  Jacques,  Canfield,  Ridgely,  and  Spillman  [6]  uses  a  gradient 
based  Sequential  Quadratic  Programming  (SQP)  numerical  optimization  routine  to  solve  the 
mixed  I  problem.  This  approach  uses  inequality  rather  than  equality  constraints.  The 
main  reason  for  this  is  that  inequality  constraint  approaches  overcome  the  numerical  drawbacks 
associated  with  using  equality  constraints.  The  approach  can  be  expressed  mathematically  as 


min 

Kstabilizing' 


2 

2 


subject  to 


(2.17) 


Now  define  the  parameter  vector  X  as 


where  a^,  ,b^, ,  and  c^,  are  the  columns  of  the  closed  loop  state  space  matrices  Aq,  Bq,  and  Q  , 

respectively.  We  can  now  use  the  SQP  algorithm  to  minimize  the  performance  index  of  this 
program  defined  above.  A  feature  of  this  algorithm  is  that  it  is  able  to  search  over  both  feasible 
and  infeasible  solutions.  Since  the  only  solutions  of  interest  are  stable  ones,  a  stability  constraint 
and  penalty  function  was  added  to  keep  the  algorithm  from  getting  lost  in  an  unstable  region. 

The  stability  constraint  is  stated  as 


g,,(^)  =  lmax 


<0 


(2.19) 


where  denotes  the  zth  eigenvalue  of  the  closed-loop  system.  The  penalty  function  added 

to  the  objective  function  is  simply  the  square  of  the  stability  constraint.  This  assures  continuous 
derivatives  at  the  stability  boundary.  The  analytical  gradients  of  the  objective  as  well  as  the 
constraints  are  derived  in  Jacques,  et  al  [6].  The  normal  approach  is  to  compute  the  H2  optimal 
controller  which  gives  us  a  starting  point  on  the  Pareto  curve  (see  App.  B,  page  B-26).  We  then 
step  to  the  left  on  the  curve  by  decrementing  y  until  we  approach  y  .  At  each  y  step,  SQP  will 

determine  the  controller  which  minimizes  a  subject  to  the  constraint  y  .  When  a  is  minimized 

as  much  as  possible  the  controller  is  saved  and  ,  ct  )  is  plotted  on  the  Pareto  curve.  It 
should  be  noted  here  that  this  method  is  not  guaranteed  to  converge  to  a  global  minimum  at  each 
y  step.  This  approach  is  implemented  in  MATLAB^*'^  through  the  MXTOOLS™  Toolbox  and 
is  available  via  anonymous  ftp  at  the  Air  Force  Institute  of  Technology. 
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2.5  Mixed  H2 1  \x  Control  Optimization. 

This  is  just  an  extension  of  the  mixed  H2 1  theory  already  discussed.  It  follows  the 
theoretical  framework  for  H2 1  //„  control  optimization  that  was  developed  in  section  B.4, 
Appendix  B.  H2  optimization  is  done  as  before,  but  now  the  part  of  the  problem  is  modified 
by  the  D  matrices  discussed  in  Section  B.3,  Appendix  B.  First,  p,  synthesis  is  done  on  the 
portion  of  the  problem  dealing  with  the  uncertainties.  Instead  of  being  interested  in  the  resulting 
controller,  we  are  now  interested  in  the  D  scales.  They  will  be  absorbed  into  the  open  loop  P  of 
the  //„  problem,  and  that  will  become  our  new  problem.  Our  new  P  is  calculated  using  the 
formula 

P  =  £)?„£>■’  (2.20) 

Now  [|7^^||^  in  Figure  B6  becomes  .  Other  than  this  modification,  there  is  no 

difference  from  the  H2 1  optimization  techniques  previously  discussed  in  Appendix  B. 
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III.  Problem  Formulation 


Now  that  the  mathematical  theory  has  been  established,  the  design  application  that  is 
used  in  the  rest  of  this  thesis  can  be  outlined.  However,  before  the  details  of  the  design  process 
are  discussed,  it  is  instructive  to  review  the  longitudinal  parameters  of  interest  during  the  landing 
phase  of  flight.  Section  3.1  contains  a  brief  overview  of  how  a  pilot  normally  accomplishes  a 
landing.  The  rest  of  the  chapter  details  all  of  the  problem  formulations  that  were  used  to  produce 
the  design  results  discussed  in  the  next  chapter.  The  initial  system  definitions  are  those  common 
to  all  of  the  problem  formulations.  Then  the  single  norm  subproblem  setups  will  be  discussed 
(/fj ,  .  and  /j)  followed  by  the  H.^  /  and  H^l  ^  problems.  The  last  part  of  the  chapter  will 

summarize  significant  findings  relating  to  the  setup  phase  as  a  prelude  to  the  design  results. 

3. 1  Background  on  Landing  Phase 

MIL-STD-1797A  [8]  defines  the  landing  phase  as  those  maneuvers  which  require 
precise  flight  path  control  using  gradual  maneuvers  during  the  terminal  phases  of  flight.  Precise 
tracking  tasks  generally  require  high  open  loop  (here,  open  loop  refers  to  the  aircraft  and  flight 

control  system  without  the  pilot)  system  stability  and  high  short  period  damping,  .  This 
enables  the  pilot  to  track  high  frequency  inputs  and  reject  disturbances  without  unacceptable 
oscillations  due  to  low  ^ in  the  system  [15].  In  the  initial  stage  of  the  approach,  the  pilot 

adjusts  pitch  angle,  9 ,  to  control  the  flight  path  angle,  y  .  Throttle  position,  is  used  to  control 
airspeed.  As  the  pilot  nears  the  touchdown  point,  the  throttle  is  moved  to  idle.  The  control 
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inputs  now  affect  the  flight  path  angle,  y  ■  The  pilot's  goal  is  to  smoothly  transition  y  to  zero  at 
touchdown[16].  Figure  1  show  the  definitions  of  pitch  attitude  angle  and  flight  path  angle. 


3.2  System  Definitions 

The  plant  model  used  throughout  this  thesis  is  a  short  period  model  of  the  Calspan 
Variable  Stability  Learjet  simulating  an  F-16.  The  state  space  model  is  defined  as 
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where  the  pitch  rate,  q,  is  in  deg/sec,  and  angle  of  attack,  a ,  elevator  input,  ,  are  in  degrees. 
This  plant  system  can  also  be  represented  in  terms  of  the  following  transfer  function 


GM  - 


g 

IA1 

- 1 

(3.2) 
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A  short  period  only  model  was  used  because  Calspan’s  experience  indicated  that  the  phugoid 
mode  did  not  have  a  significant  effect  on  the  handling  qualities.  The  second  order  plant  model 
also  offered  the  advantage  of  keeping  the  resulting  control  laws  lower  order.  The  second  order 
model  for  the  actuator  dynamics  was  also  supplied  by  Calspan  and  the  state  space  model  is  given 
by 


■-17.915 

-49.154' 

■-4.9575' 

49.154 

-66.085 

4.9575 

=[-4.9575  -4.9575]x„„ +[0]^„ 


(3.3) 


The  actuator  system  is  also  referred  to  as  G^ct  is  defined  as 


'A. 

J^act 

- 1 

(3.4) 


The  next  system  definition  common  to  all  of  the  design  problems  is  the  weight  on  the 
reference  input,  .  It  was  not  desired  to  shape  the  frequency  content  of  the  reference  input  in 
the  design  process  so  was  set  equal  to  1 . 

The  final  element  that  is  common  to  all  designs  is  the  "ideal"  closed-loop  model  that  is 
used  in  the  model  following  portion  of  the  problem.  The  reference  input  to  the  controller  for  this 
problem  is  a  pitch  rate  command.  So  we  will  need  to  choose  a  desired  closed  loop  model  of  the 
pitch  rate  to  commanded-pitch-rate  transfer  function  to  follow.  In  this  case,  good  handling 
qualities  predictions  and  good  tracking  of  a  pitch  rate  input  were  the  primary  figures  of  merit  in 
the  selection  of  an  ideal  model.  It  was  also  desired  to  keep  the  order  of  the  model  as  low  as 
possible  so  the  order  of  the  resulting  controller  would  not  get  too  large  for  practical 
implementation  in  the  flight  test  phase.  A  first  order  model  was  found  that  satisfied  all  of  these 
requirements  and  is  defined  as 
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(3.5) 
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3.3  i/j  Problem 

The  problem  is  an  LQG  design  that  minimizes  the  output  states  (q  and  a ),  elevator 
control  rate  usage,  elevator  usage,  and  pitch  rate  tracking  error  to  wind  disturbances  and  sensor 
noise.  The  problem  setup  is  depicted  in  Figure  3.2. 


Figure  3.2.  H2  Problem  Setup 

It  is  important  to  now  spend  some  time  discussing  the  various  weight  selections  that 
were  used  in  this  H2  problem.  These  weights  determine  what  performance  measures  are 
emphasized  in  the  resulting  controller  (tracking  vs.  noise  rejection,  for  instance).  It  turns  out  that 
the  choices  of  the  control  rate  weight,  Wcr,  and  the  tracking  weight,  Wp,  have  a  significant 
impact  on  handling  qualities  predictions.  The  control  rate  weight  selection  directly  impacted  the 
bandwidth  of  the  resulting  closed-loop  system.  In  general,  larger  control  rate  weights  produced 
closed-loop 


3-4 


systems  with  smaller  bandwidth  that  rolled  off  rapidly  in  phase  at  frequencies  above  the 
bandwidth  frequency.  This  means  that  the  pilot  has  less  ability  to  increase  the  bandwidth  of  the 
system  through  increasing  his  gain  (i.e.  he  can’t  aggressively  track  without  causing  instability  in 
the  system).  Continuing  to  increase  the  control  weight  rate  will  eventually  decrease  the 
bandwidth  and  increase  the  phase  rolloff  to  the  point  where  level  1  handling  qualities  predictions 
are  no  longer  possible.  This  will  be  illustrated  further  in  Chapter  4  when  looking  at  the  Hi 
problem  results.  A  potential  solution  to  this  problem  would  be  to  select  a  very  small  weight  on 
control  rate  use.  There  are  two  problems  with  this  approach.  First,  a  very  small  weight  on 
control  use  will  usually  cause  rate  limiting  of  control  actuators,  which  can  lead  to  instability. 
Second,  widening  the  bandwidth  of  the  system  will  increase  its  susceptibility  to  sensor  noise 
corruption. 

The  tracking  weight  had  a  more  indirect  effect  on  the  bandwidth  of  the  closed-loop 
system.  The  higher  the  DC  gain  of  this  weight,  the  better  the  tracking  performance  of  the  closed- 
loop  system.  However,  the  price  to  be  paid  for  this  good  tracking  was  rate  limit  saturation  of  the 
actuators.  In  order  to  rid  the  system  of  the  rate  limit  saturation,  the  control  rate  use  weight  had  to 
be  increased.  This  narrowed  the  bandwidth  of  the  system,  and  increased  phase  rolloff,  which 
produced  a  degradation  in  handling  qualities  predictions.  The  bottom  line  is  that  there  ends  up 
being  a  trade-off  between  tracking  performance  and  handling  qualities  in  the  selection  of  these 
weights.  With  this  in  mind,  the  tracking  weight  is  defined  as 


"  s  +  4 


(3.6) 


The  smallest  corresponding  rate  limit  weight  that  just  prevents  the  system  from  rate  limiting  for 
a  5  deg,  4  sec  pulse  input  is 
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0,000667  0 

^ in 


(3.7) 


The  control  usage  weight  had  very  little  effect  on  the  output  of  the  overall  system  because 
the  rate  limit  constraint  was  more  restrictive  for  this  flight  condition  and  control  input. 

However,  this  weight  needs  to  be  included  because  the  rate  limit  may  not  always  be  the 
dominant  constraint  depending  on  flight  condition  and  control  input.  The  saturation  limit  of  the 
elevator  on  an  F-16  is  +/-  25  degrees.  Therefore  the  control  use  weight  was  set  at  the  inverse  of 
this  and  is  given  by 

pr„=l/25  (3.8) 

The  next  weight  that  was  chosen  was  the  state  weighting  on  the  outputs  q  and  a .  Since 
one  of  the  objectives  here  is  to  minimize  the  states  given  wind  and  sensor  noise  inputs,  we  want 
to  just  minimize  the  outputs  of  the  states  themselves.  For  this  reason  Wgi  was  chosen  as  an 
identity  matrix. 

ri  01 

The  next  set  of  weights  to  be  discussed  are  those  associated  with  the  wind  turbulence  and 
sensor  noise.  Wind  turbulence  manifests  itself  as  an  angle  of  attack  disturbance  in  an  aircraft. 

For  this  reason  it  makes  sense  to  weight  the  angle  of  attack  states  in  the  plant  A-matrix.  It  was 
found  through  iteration  (while  holding  other  weights  fixed)  that  a  representative  turbulence 
noise  is  produced  by  introducing  the  following  weight  as  an  angle  of  attack  disturbance: 

r  =0,0316x^^(:,2)  (3.9) 
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The  sensor  noise  weightings  corrupted  the  state  feedback  measurements,  q  and  a ,  and 
effectively  simulates  the  noise  that  is  present  on  those  measurements.  A  static  weight  was  used 
here  and  was  again  found  by  iteration  (again,  while  holding  other  weights  fixed).  The  sensor 
noise  weight  is  given  by 

.1  01 

^■-[o  lJU 

This  says  that  the  noise  affects  the  pitch  rate  measurement,  q,  by  10%  over  all  frequencies  and 
the  angle-of-attack  measurement,  a ,  by  100%  over  all  frequencies.  It  is  worth  spending  a 
moment  discussing  these  extremely  conservative  weights.  Although  sensor  noise  is  mainly  a 
high  frequency  disturbance,  static  weights  were  found  to  work  well  and  had  the  added  benefit  of 
keeping  the  order  of  the  H2  problem  down  (This  will  become  important  later  in  the  mixed 

problem).  The  magnitude  of  the  sensor  noise  weights  were  given  the  values  above  to 
optimize  the  noise  rejection  performance  of  the  H2  design.  One  way  to  get  better  noise  rejection 
performance  from  an  LQG  design  is  to  input  a  higher  level  (magnitude  and  frequency)  of  noise 
than  the  system  would  actually  expect  to  experience.  However,  the  improvement  in  noise 
rejection  usually  comes  at  the  expense  of  other  performance  specifications  like  command 
tracking.  For  this  design  the  values  of  provided  excellent  noise  rejection  with  no  appreciable 
degradation  in  tracking. 

3.4  Problem 

The  problem  formulation  included  three  different  objectives;  command  tracking, 
maximizing  output  stability  margins,  and  making  the  system  robust  to  uncertainty  at  the  input  to 
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the  plant  by  accounting  for  neglected  high  frequency  plant  dynamics  and  other  modeling  errors 
that  may  be  present  in  the  plant.  A  block  diagram  of  this  setup  is  shown  in  Figure  3.3 


Figure  3.3.  Problem  Setup 

The  tracking  objective  is  controlled  by  minimizing  the  transfer  function  of  the  weighted 
output  63  and  input  ds.  The  output  stability  margins  are  addressed  by  the  minimizing  the 
weighted  complimentary  sensitivity  transfer  function,  e2d2.  Finally,  the  uncertainty  robustness 
objective  is  contained  in  the  eidi  transfer  function.  The  weights  associated  with  each  of  these 
outputs  will  now  be  discussed. 

The  tracking  weight  for  the  problem  was  much  more  stringent  than  in  the 
problem.  As  mentioned  previously,  the  tracking  weight  had  an  effect  on  the  handling  qualities 
predictions  in  the  problem  by  affecting  the  bandwidth  of  the  closed- loop  system.  This  is  not 

the  case  in  this  work  because  the  problem  is  used  only  as  a  constraint  in  the  mixed 
problem.  It  was  found  in  this  work  that  the  handling  qualities  predictions  were  almost 
completely  influenced  by  the  problem.  Therefore,  no  tradeoff  between  handling  qualities 
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and  the  various  weight  selections  had  to  be  made  in  the  problem.  Thus  the  tracking  weight 


is  defined  as 


5+.004 


(3.11) 


The  output  vector  stability  margins  are  maximized  by  minimizing  the  weighted 
complementary  sensitivity  of  the  measurement  outputs  of  the  plant.  The  reader  is  invited  to  see 
Franklin,  et  al.  [17]  for  a  detailed  discussion  of  vector  stability  margins.  There  was  no  need  to 
frequency  weight  the  output  measurements  in  order  to  minimize  complementary  sensitivity,  so 
this  weight  was  set  to  identity  and  is  expressed  as 


1  0 
0  1 


(3.12) 


The  input  uncertainty  weight  was  chosen  to  account  for  neglected  high  frequency 
dynamics  and  modeling  errors  in  the  actuator  and  plant.  By  minimizing  this  objective,  the 
designer  should  be  able  to  improve  the  controller  performance  for  an  “off-nominal”  plant.  The 
uncertainty  weight  is  expressed  as 


_  3;y  +  30 
”s  +  250 


(3.13) 


The  final  item  that  needs  to  be  addressed  with  regards  to  this  problem  is  the  absence 
of  control  rate  and  control  use  weights.  Recall  from  Chapter  2/Appendix  B,  that  the  constraints 
in  the  mixed  problem  come  only  from  the  portion  of  the  problem.  Therefore  the 

rate  limit  and  control  use  constraints  do  not  need  to  be  included  in  the  setup  because  they  are 
taken  eare  of  by  the  problem. 
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3.5  n  Problem 


The  /u  problem  is  setup  exactly  the  same  as  the  problem  discussed  in  section  3.3 
other  than  the  “D”  scales  associated  with  the  uncertainty  structure,  A  (see  Figure  3.2).  Recall 
from  Chapter  2,  that  ju  has  the  advantage  of  handling  multiple  uncertainties  less  conservatively 
than  the  problem,  ju  also  has  the  advantage  of  guaranteeing  robust  performance  for  values  of 
/i  <  1,  see  [1].  This  means  that  performance  specifications  can  be  met  while  remaining  robustly 
stable  for  all  perturbations  contained  in  the  structure  of  A.  For  this  problem  setup,  the  A 
structure  contains  the  tracking,  output  stability  margins,  and  input  uncertainty  objectives. 

3.6  Mixed  H,  /  Problem 

Recall  the  mixed  problem  can  be  expressed  as 

i"!  11^- 11^  ’  P'«'  IL  -  y 

where  ^adm  is  the  set  of  stabilizing  controllers  of  some  fixed  order.  This  problem  can  have  a 
single  constraint  transfer  function  or  multiple  constraint  transfer  functions.  Allowing  the 
problem  to  be  set  up  with  a  single  or  multiple  constraints  gives  the  designer  great  flexibility  to 
tailor  the  problem  to  whatever  objective(s)  are  important  for  a  particular  application.  For 
example,  a  control  engineer  has  an  LQG  tracking  design  that  produces  good  command  following 
and  noise  rejection,  but  insufficient  output  stability  margins.  An  constraint  minimizing 
output  sensitivity  or  output  complementary  sensitivity  could  be  employed  in  a  mixed 
problem  formulation  to  correct  the  stability  margin  deficiency  of  the  problem.  The  reader  is 

reminded  that  these  constraints  are  allowed  to  be  singular  as  all  of  the  regularity  constraints  (i.e. 
certain  matrices  must  have  a  nonzero  determinant)  are  handled  by  the  portion  of  the  problem. 
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The  setup  of  the  mixed  problem  is  provided  in  Figure  3.4.  The  inputs  to  the  system  are 

w  (unit  intensity  white  Gaussian  noise),  and  d,  an  unknown  but  bounded  energy  input.  The 
outputs  z  and  e  are  the  signals  whose  two-norms  are  either  minimized  or  constrained.  All  of  the 
weights  used  in  the  mixed  problem  are  carried  over  from  the  single  norm  problem 

formulations. 


3. 7  Mixed  Problem 

Analogous  to  the  single  norm  case,  the  mixed  !  ju  is  setup  the  same  way  as  the  mixed 
problem.  However  //-synthesis  must  now  be  performed  to  determine  the  optimal  “D” 
scales  to  include  in  the  //  constraint.  We  are  not  interested  in  the  controller  generated  by  //- 
synthesis,  only  the  “D”  scales.  In  essence  we  still  have  a  mixed  problem  that  now 

includes  scalings  associated  with  uncertainty  structure  of  the  constraint.  The  problem  setup, 
minus  the  “D”  scales  is  shown  in  Figure  3.3. 

3.8  Summary 

This  chapter  detailed  the  single  norm  subproblem  setups  and  ju)  followed  by 

the  mixed  and  H^l  ^  problems.  The  transfer  function  matrices  common  to  all  problems 
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and  the  weights  on  various  objectives  of  each  problem  formulation  were  also  presented.  These 
will  be  used  to  build  the  system  matrices  that  are  utilized  in  the  design  and  simulation  results 
discussed  in  Chapter  4. 

The  most  significant  finding  from  this  chapter  is  that  the  control  rate  usage  and  tracking 
performance  weights  appears  to  influence  the  handling  qualities  predictions  of  a  control  design 
by  controlling  the  bandwidth  of  the  closed  loop  system. 

Also,  of  note,  is  that  the  LQG  problem  can  be  “tricked”  into  providing  better  noise 
rejection  by  increasing  the  weight  on  the  noise  signals  being  input  to  the  system.  However,  this 
comes  at  the  price  of  reduced  performance  from  other  objectives  and  a  tradeoff  must  be  made  by 
the  designer. 

In  order  to  allow  the  reader  to  easily  track  the  evolution  of  the  control  designs  (and 
methods)  through  the  flight  test,  design  summary  tables  will  be  provided  at  the  end  of  each 
chapter  from  this  point  forward. 


TABLE  3.1 


SUMMARY  OF  DESIGN  METHODS  AND  NUMBERS  TO  BE  PRODUCED 
FOR  FLIGHT  TEST  CONSIDERATION 


Control  Design  Method 

Number  of  Designs 

1 

1 

JU 

1 

H,  ///. 

1 

HJm 

1 
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IV.  Design  and  Simulation  Results 


This  chapter  presents  the  phase  I  control  designs  that  are  candidates  for  flight  testing  in 
phase  II.  There  are  five  pitch  rate  designs,  and  two  angle  of  attack  control  designs.  There  was 
also  a  modification  made  to  the  simulation  that  converted  the  five  pitch  rate  designs  to  angle  of 
attack  designs.  This  will  be  discussed  further  when  the  angle  of  attack  designs  are  presented. 
Some  of  the  twenty  other  designs  that  were  not  good  enough  to  be  flight  test  candidates  will  also 
be  mentioned  where  appropriate  to  help  illustrate  and  support  the  analytical  findings  contained 
herein.  It  should  also  be  noted  here  that  time  and  financial  constraints  may  only  allow  flight 
testing  of  3-4  of  the  flight  test  candidate  designs.  However,  all  seven  are  included  for  the  sake  of 
analytical  completeness.  The  angle  of  attack  designs  were  added  very  late  in  Phase  I. 

Discussions  with  Calspan  engineers  indicated  that  the  pitch  rate  command  following  designs 
may  exhibit  deficient  handling  qualities  when  the  pilot  is  required  to  be  tightly  in  the  control 
loop  (e.g.,  turbulent  air),  and  that  an  angle  of  attack  system  would  be  preferable.  However,  a  test 
program  conducted  by  the  F-16  Combined  Test  Force  (CTF)  indicates  that  pilots  preferred  a 
pitch  rate  system  for  landing  [18],  Because  rejection  of  wind  turbulence  is  included  directly  in 
the  pitch  rate  designs,  it  is  hypothesized  that  the  presence  of  turbulence  will  not  greatly  affect 
these  designs.  For  these  reasons,  both  the  pitch  rate  and  angle  of  attack  designs  will  be 
represented  in  Phase  II.  It  is  hoped  that  there  will  be  some  discrimination  in  phase  II  as  to  which 
type  of  design  is  preferable.  The  pitch  rate  designs,  simulations,  and  handling  qualities 
predictions  will  be  discussed  first.  Next,  the  angle  of  attack  designs,  simulations,  and  handling 
qualities  will  be  discussed.  Finally,  a  summaiy  of  Phase  I  findings  along  with  a  summary  of  the 
designs  picked  for  flight  test  will  be  given  as  a  prelude  to  Phase  II  and  Chapter  V. 
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4. 1  Pitch  Rate  Design  Results 


The  five  pitch  rate  control  designs  that  will  be  discussed  consist  of:  a  single  norm 

design,  three  mixed  designs,  and  one  mixed  H^l  jJ.!  design.  The 

following  table  defines  the  constraints  and  objectives  used  in  the  mixed-norm  designs 


TABLE  4. 1 

SUMMARY  OF  CONSTRAINTS  AND  OBJECTIVES  FOR  THE  MIXED 

NORM  PROBLEMS 


Constraint 

Objective 

Command  tracking  performance 

Plant  Input  Uncertainty  Robustness 

T 

Output  Stability  Margins 

Combination  of  and 

The  objectives  in  Table  4.1  were  handled  separately  in  the  mixed  I design 
because  the  resulting  design  was  better  than  trying  to  include  all  three  objectives  in  one 
constraint.  There  are  two  of  reasons  for  this.  First,  the  reader  will  recall  that  problem 
formulations  are  guaranteed  to  only  handle  one  uncertainty  in  a  nonconservative  manner.  The 
objectives  in  Table  4.1  can  be  looked  at  as  uncertainties.  So  by  virtue  of  including  all  three 
objectives  in  one  constraint,  the  resulting  design  may  be  too  restrictive.  Second,  the 
objectives  may  have  the  same  frequency  ranges  in  which  they  compete.  By  breaking  our  multiple 
objective  constraint  into  three  single  objectives,  we  can  see,  and  control,  where  the 
conservativeness  is  introduced  by  tightening  or  relaxing  the  individual  constraints.  In  the  case 
where  all  the  objectives  are  together,  the  most  difficult  one  to  achieve  will  limit  the  performance 
of  the  others.  A  good  example  of  this  will  be  illustrated  in  the  mixed  //j  /  //  /  design. 
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Tables  4.2  and  4.3  summarize  the  results  of  the  five  pitch  rate  designs.  Table  4.2 
compares  the  2-norms  and  oo  -norms  of  the  various  designs  while  Table  4.3  compares  input  and 
output  vector  gain  and  phase  stability  margins  [18].  Following  the  tables,  there  are  discussions 
of  each  design. 


TABLE  4.2 


PITCH  RATE  DESIGN  DATA-lr^pf^l^  vs.  TRACKING, 

INPUT  ROBUSTNESS,  AND  OUTPUT  STABILITY  CONSTRAINTS 


Case 

Design  Type 

Il^.»ll2 

IL 

|5,',L 

1-1 

2.518 

528.970 

13.280 

2.950 

1-2 

5.022 

1.400 

4.000 

1.650 

1-3 

^2  !  ^|</i  ^  ^ 

5.810 

1.400 

2.800 

1.650 

1-4 

H2lT,^dJT,^dJT,^d, 

3.940 

1.400 

9.925 

1.650 

1-5 

8.410 

34.220 

.990 

1.613 

TABLE  4.3 

PITCH  RATE  DESIGN  DATA  -  STABILITY  MARGIN  COMPARISON 


Case 

Complement. 
Sensitivity 
Vector  Gain 
Margins(dB) 
Input  of  Plant 

Sensitivity 
Vector  Gain 
Margins  (dB) 
Input  of  Plant 

Phase 
Margins 
(deg) 
Input  of 
Plant 

Complement. 
Sensitivity 
Vector  Gain 
Margins  (dB) 
Output  of 
Plant 

Sensitivity 
Vector  Gain 
Margins  (dB) 
Output  of  Plant 

Phase 
Margins 
(deg) 
Output  of 
Plant 

1-1 

[-6.53, 3.68] 

[-5.88,  29.72] 

57.85 

[-3.62,2.55] 

[-2.84, 4.24] 

1-2 

[-20.11,5.58] 

[-5.13, 14.20] 

53.57 

[-14.59,5.17] 

[-5.11,  14.02] 

1-3 

[-21.62,5.65] 

[-4.79, 11.61] 

54.59 

[-14.94,  5.21] 

[-4.78,  11.52] 

48.47 

1-4 

[-16.25,  5.32] 

[-5.46, 18.06] 

51.89 

[-14.18,5.13] 

[-4.99,  13.07] 

47.44 

1-5 

[-10.56,4.63] 

IBKsMililigfl 

45.42 

[-8.35,4.18] 

[-4.78,  11.52] 

4. 1. 1  The  H2  Design  (Case  1-1) 

This  design  is  included  as  a  flight  test  candidate  because  it  is  the  basis  for  all  of  the 
mixed  and  !  /a  designs  that  follow.  The  large  value  of  |fe^3|L  Table  4.2 
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indicates  that  this  design  does  not  track  pitch  rate  command  inputs  well  at  all.  Additionally,  the 
value  of  L  indicates  that  this  control  design  is  not  robustly  stable  to  input  perturbations  in 
the  class  included  in  the  problem  formulation  because  the  infinity-norm  is  not  less  than  1 . 
We  would  not  expect  this  as  we  did  not  include  this  type  of  weighting  in  the  problem 
formulation.  The  final  deficiency  of  this  design  can  be  seen  in  Table  4.2.  The  output  stability 
margins  are  unacceptably  low.  Good  stability  margins  for  aircraft  applications  are  [-6.0,12.0]dB 
of  Gain  Margin  and  +/-  30  deg  of  Phase  Margin  [19].  The  tracking  and  output  stability  margin 
problems  must  be  addressed  in  order  to  have  a  good  design  while  the  input  uncertainty  weighting 
is  an  artificial  constraint  that  falls  into  the  "nice  to  have"  category  for  this  problem.  The  first 
question  to  ask  is  how  to  improve  tracking  and  output  stability  margin  performance?  Tracking 
performance  might  be  improved  by  simply  tightening  the  tracking  constraint  in  the  original 
problem  formulation.  Likewise,  the  sensor  noise  weights  could  be  adjusted  to  improve  our 
output  stability  margins.  Both  of  these  approaches  were  evaluated  and  resulted  in  a  significant 
degradation  in  noise  rejection  performance  which  led  to  unsatisfactory  results.  Another 
approach  would  be  to  include  these  performance  objectives  in  constraints  that  are  part  of  a 
multiobjective  tradeoff  such  as  the  problem  discussed  next. 


4. 1.2  The  /T;^^  Designs  (Cases  1-2  through  1-4). 

These  designs  represent  varying  tradeoffs  between  noise  rejection  performance 
quantified  by  ^nd  input  uncertainty  robustness  quantified  by  ||  .  The  reader  will 

note  that  the  norms  for  the  tracking  and  output  stability  margin  constraints  are  all  the  same  for 
each  case.  It  was  found  that  these  norm  values  produced  good  tracking  and  output  stability 
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margins.  These  constraints  were  then  "frozen"  while  trying  to  minimize  the  input  uncertainty 
constraint.  As  in  the  H2  design,  robust  input  stability  is  not  guaranteed  for  any  of  these  designs. 
The  singular  value  plots  of  each  closed  loop  constraint  (where  the  //j  optimal  controller  is  used 
to  close  the  loop)  in  Figure  4.1  show  why  there  are  problems  optimizing  the  robust  stability 


Figure  4.1.  Singular  Value  Plots  of  Constraints 

constraint.  The  tracking  and  output  stability  margin  constraints  (  and  )  are  composed  of 

low  frequency  dynamics  while  the  input  robustness  constraint  (  )  is  composed  of  high 

frequency  dynamics.  This  competition  between  objectives  indicates  that  we  cannot  get  robust 
input  stability,  at  90  rad/sec  and  above,  without  sacrificing  tracking  and/or  output  stability 
margin  performance.  A  technique  to  potentially  improve  these  designs  is  discussed  next. 
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4.1.3  The  Mixed  H2  /  /  fj.  Design  (Case  1-5). 

This  design  was  developed  to  try  to  improve  the  results  obtained  in  cases  1-2  through  1- 
4.  The  approach  was  to  combine  the  two  low  frequency  constraints  in  a  //  -  problem  formulation 
while  leaving  the  robust  input  stability  constraint  separate.  Then  a  mixed  design  was 

developed  that  guaranteed  robust  input  stability.  The  value  of  ||  =  0.2  resulted  from  the 
fi  design.  This  is  well  below  the  0.99  requirement  to  guarantee  robust  input  stability.  Now  this 
//-  constraint  was -appended  to  the  H2  /  problem  to  account  for  the  command  tracking  and 
output  stability  margin  objectives  (The  reader  is  reminded  that  this  is  not  an  explicit  // 
constraint,  but  rather  an  constraint  that  is  optimally  scaled  through  //-synthesis).  The 
approach  was  then  to  relax  the  robust  input  stability  constraint  co-norm  to  .99  while  minimizing 
the  //-constraint  00-norm  to  produce  a  design  that  meets  the  tracking,  output  stability  margin,  and 
robust  input  stability  objectives.  The  potential  advantage  of  this  formulation  vs.  the  three 
individual  constraints  was  that  the  //  constraint  would  trade  off  the  included  performance 
objectives  in  an  optimal  manner,  rather  than  requiring  the  designer  to  tradeoff  the  objectives 
individually  as  had  to  be  done  on  the  previous  cases.  The  design  results  did  support  the 
hypothesis  of  this  problem  formulation.  Table  4.2  shows  that  while  we  have  robust  input 
stability,  our  command  input  tracking  is  not  as  good  as  what  we  were  able  to  get  in  the  previous 
designs.  This  goes  back  to  the  fact  that  no  matter  how  you  set  up  this  problem,  the  robust 
stability  objective  competes  with  tracking  and  output  stability  margin  objectives.  The  bottom 
line  is  that  tracking  and/or  output  stability  margin  performance  is  going  to  have  to  be  sacrificed 
to  guarantee  robust  input  stability.  However,  the  simulations  show  that  this  design  does  indeed 
produce  good  command  tracking.  Unfortunately,  this  controller  had  one  pole  in  the  Right-Half 
Complex  Plane  (RHP)  and  was  disqualified  as  a  flight  test  design  based  on  the  problems  HAVE 
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INFINITY  had  with  implementing  control  laws  with  poles  in  the  RHP.  An  option  here  would  be 
to  pick  a  controller  that  is  "lower"  on  the  Edgeworth-Pareto  curve  (see  Appendix  B)  in  Figure  4.2 
that  has  all  Left-Half  Complex  Plane  (LHP)  poles. 


Hinf  norm  of  mu  constraint 


Figure  4.2.  Edgeworth-Pareto  Curve  for  mixed  H2I  /d  problem 
The  point  that  corresponds  to  this  type  of  controller  is  the  marked  point  on  the  left  side 
of  the  Pareto  curve  in  Figure  4.2.  Note  that  the  curve  shown  in  Figure  4.2  is  not  monotonically 
increasing  in  as  the  theory  dictates.  Since  the  stabilizing  control  law  is  found  iteratively,  the 
SQP  algorithm  can  get  ‘lost’  from  one  solution  point  to  the  next  and  not  necessarily  produce  the 
best  solution  from  a  theoretical  point  of  view.  Multiple  runs  of  the  algorithm  would  be  required 
to  eliminate  this  problem.  Since  the  control  laws  of  interest  in  this  case  were  not  contained  in 
the  non-monotonic  region,  the  optimization  was  not  re-accomplished. 
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The  norms  and  vector  margins  of  this  controller  are 


r„||,  =5.4891,  I?;,,,,  I,  =48.99.  =  .99,  ai.d|7;^JL  =  '•™0 


Case 

Complement. 
Sensitivity 
Vector  Gain 
Margins  (dB) 
Input  of  Plant 

Sensitivity 
Vector  Gain 
Margins  (dB) 
Input  of  Plant 

Phase 

Margins  (deg) 
Input  of  Plant 

Complement, 
Sensitivity 
Vector  Gain 
Margins  (dB) 
Output  of 
Plant 

Sensitivity 
Vector  Gain 
Margins  (dB) 
Output  of 
Plant 

Phase 

Margins  (deg) 
Output  of 
Plant 

l-5a 

[-22.61, 5.69] 

[-4.89,  12.24] 

55.16 

[-11.55, 4.79] 

[-4.88, 12.16] 

44.27 

The  stability  margins  of  this  controller  are  excellent,  but  there  is  degradation  in  tracking 
performance  as  indicated  by  the  increase  in  The  next  step  would  be  to  evaluate  this 

design's  tracking  performance  to  see  if  it  is  an  acceptable  design.  This  is  accomplished  in  the 
following  section. 

4.1.4  Simulations. 

Several  simulations  were  carried  out  to  validate  the  short  period  longitudinal  controller 
designs.  Unfortunately,  the  nonlinear  aircraft/atmospheric  model  of  the  Calspan  LearF-16  was 
not  available  to  perform  dynamic  simulations  and  validation  of  the  controller  designs. 

Therefore,  Matlab's™  simulation  toolbox  SIMULINK™  [20]  was  used.  The  simulations  that 
were  performed  were  "static"  in  that  the  controllers  were  tested  only  at  the  design  points.  For 
small  changes  in  pitch  rate  (which  we  would  expect  in  the  approach  and  landing  phase,  and  is  an 
assumption  of  the  linearized  equations  of  motion),  these  simulations  should  provide  a  realistic 
transient  response  as  long  as  the  state  and  input  variables  remain  relatively  small.  Figure  4.3 
illustrates  the  SIMULINK'^'*'  setup  used  to  test  the  control  designs  developed  in  cases  1-1 
through  l-5a.  As  illustrated,  the  simulation  included  wind  and  measurement  noises,  as  well  as 
rate  and  saturation  limiters  on  elevator  deflection.  A  Dryden  wind  model  [21]  was  fed  by  white 
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noise  to  produce  moderate  turbulence  conditions  in  approach  and  landing  as  specified  in  [8]. 

The  measurement  noise  inputs  consisted  of  both  filtered  additive  and  multiplicative  components 
that  were  fed  by  white  noise.  The  breakdown  was  98%  multiplicative  and  2%  additive.  The 
magnitude  of  multiplicative  noise  on  the  measurement  increased  as  the  magnitude  of  the  output 
signal  increased.  However,  when  there  was  no  output  signal,  there  was  no  multiplicative 
component  of  measurement  noise.  There  was  always  additive  noise  present  in  the  system 
representing  the  steady  state  inaccuracies  of  measurement  devices.  The  high-pass  Butterworth 


0- 

Clock 


Tbiiti  I 
Time 


Figure  4.3.  SIMULINK™  Setup  for  Pitch  Rate  Control  Design  Evaluation 

filters  were  included  because  measurement  noise  is  substantially  a  high  frequency  phenomenon. 
The  "gain"  blocks  immediately  following  the  Butterworth  Filters  contain  the  various  gains  used 
to  produce  the  correct  amount  of  sensor  noise  on  each  measurement  and  a) .  Very  little  data 
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is  available  in  academic  texts  on  typical  sensor  measurement  accuracies.  The  data  used  in  these 
simulations  was  obtained  from  McClean  [22]  and  USAF  TPS  flight  test  data  for  their  aircraft. 
The  final  elements  included  in  the  simulation  to  add  realism  are  a  pure  transport  delay  of  0.016 
seconds  and  a  stick  dynamics  model.  The  delay  accounts  for  overall  system  delays  in  signal 
processing,  etc.,  while  the  stick  dynamics  model  accounts  for  pilot  interface  delays  with  the 
flight  control  system.  The  delay  was  chosen  based  on  an  sensor  sample  rate  of  63  Hz  for  an 
operational  Block  30+  F-16.  Finally,  the  simulation  model  includes  only  short  period  dynamics. 
This  was  done  on  the  advice  of  Calspan  engineers,  because  a  good  phugoid  model  for  the 
LearF16  was  not  available  and  flight  test  experience  with  this  model  has  shown  that  the  phugoid 
mode  is  an  insignificant  factor  in  influencing  handling  qualities.  The  handling  qualities 
prediction  for  each  design  is  plotted  in  Figure  4.16  following  the  simulation  plots.  Comments 
regarding  handling  qualities  will  be  made  as  the  simulations  are  discussed  for  each  design. 

The  time  responses  in  Figures  4.4  and  4.5  show  that  the  controller  does  not  track  the 
command  exactly,  as  predicted,  however  the  error  is  only  on  the  order  of  7%  and  the  noise 
rejection  performance  is  excellent.  This  design  predicts  level  1  handling  qualities  even  with  its 
tracking  deficiencies  because  it  has  good  phase  rolloff  and  bandwidth  characteristics.  Figures 
4.6  and  4.7  represent  case  1-2  and  show  that  this  design  provides  excellent  tracking  with  only  a 
small  sacrifice  in  noise  rejection  performance.  However  the  handling  qualities  are  predicted 
level  2  for  this  design.  In  order  to  get  this  excellent  tracking,  the  bandwidth  has  decreased  and 
hence  induced  excessive  phase  rolloff.  Figures  4.8  and  4.9  are  the  simulations  for  case  1-3.  A 
small  amount  of  tracking  performance  has  been  sacrificed  here  to  move  the  handling  qualities 
back  to  the  level  1/2  border.  Figures  4.10  and  4.1 1  show  the  simulations  for  case  1-4.  Tracking 
performance  has  been  further  sacrificed.  But  this  small  sacrifice  has  resulted  in  solid  level  1 
handling  qualities  predictions.  Finally,  Figures  4.12-4.15  cover  cases  1-5  and  l-5a.  Case  l-5a  is 
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not  shown  on  the  handling  qualities  plot  because  it  results  in  the  same  prediction  as  case  1-5. 
However  these  figures  illustrate  a  fairly  noticeable  loss  in  tracking  performance  if  we  are 
constrained  to  a  controller  that  has  all  of  its  poles  in  the  LHP.  Overall,  the  results  for  the  pitch 
rate  command  designs  were  very  good.  However,  discussions  with  Calspan  test  pilots  and 
engineers  close  to  the  end  of  phase  I  revealed  that  these  types  of  designs  may  exhibit  some 
undesirable  handling  qualities  in  the  approach  and  landing  phase  based  on  their  flight  test 
experience  [23].  This  will  be  addressed  next  in  the  angle  of  attack  design  results. 

4.2  Angle  of  Attack  Design  Results 

The  angle  of  attack  designs  are  included  to  address  potential  handling  qualities 
deficiencies  in  the  pitch  rate  designs  [23].  The  main  problem  with  the  pitch  rate  designs  for  the 
approach  and  landing  phase  can  be  seen  by  looking  at  the  pitch  angle  (0) ,  and  flight  path  angle 
(y)  plots  in  Figures  4.4  -  4.15.  Note  that  there  is  a  steady  state  value  for  both  0  and  y  after  the 
control  is  removed.  This  may  give  pilots  the  sensation  they  are  "floating",  and  will  require  them 
to  push  forward  on  the  stick  with  a  pulse  type  input  in  the  flare  to  land  the  airplane.  This 
characteristic  of  these  designs  is  caused  by  an  integral  closed  loop  pole  at  the  origin.  This 
suggests  that  putting  a  closed  loop  zero  at  the  origin  might  be  a  way  to  fix  this  problem.  Two 
approaches  to  doing  this  are  discussed.  The  first  involves  adding  a  washout  filter  and  a  lead 
filter  to  the  pitch  rate  simulations.  Note  that  the  control  designs  are  the  same  ones  from  the 
previous  section.  The  command  signal  is  being  filtered  to  produce  a  specific  type  of  response. 
The  specifics  of  the  filter  structure  and  utility  will  be  discussed  in  the  next  section.  This 
approach  was  implemented  mainly  because  of  severe  time  constraints  remaining  in  Phase  I  of 
this  research.  Figure  4.17  shows  the  simulation  setup  for  this  approach.  The  second  approach 
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involves  actually  designing  a  new  controller  that  eliminates  the  steady  state  characteristics  of  the 
pitch  and  flight  path  responses.  This  will  produce  an  angle  of  attack  following  system  which 
should  be  more  desirable  for  the  landing  task.  Again,  because  of  the  severe  time  limitations, 
only  two  designs  were  completed.  However  these  designs  are  promising  and  will  be  discussed 
along  with  the  necessary  changes  in  the  problem  setup  to  produce  them. 
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Figure  4.6.  Short  Period  Simulations  for  Case  1-2  -  Noises  off 


Figure  4.8.  Short  Period  Simulations  for  Case  1-3  -  Noises  off 
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Figure  4.13.  Short  Period  Simulations  for  case  1-5  -  noise  on 


Figure  4.16.  Handling  Qualities  Prediction  for  Cases  1-1  through  1-5 
4.2. 1  Simulation  Changes. 

As  was  mentioned  above,  the  only  change  required  to  turn  the  pitch  rate  designs  into 
angle  of  attack  following  designs  was  the  addition  of  a  washout  filter  and  lead  filter  to  the 
command  path  of  the  simulation  (see  Fig.  4.17).  Since  these  are  placed  outside  the  feedback 
loop,  there  are  no  concerns  with  regards  to  changing  the  overall  stability  of  the  system  because 
it  is  handled  exclusively  in  the  feedback  path.  These  filters  will  reshape  the  frequency  (and  time) 
responses  of  the  simulations.  Therefore  the  handling  qualities  predictions  may  not  be  the  same 
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as  they  are  without  the  filters. 
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Figure  4. 17.  SIMULINK  Structure  for  Angle  of  Attack  Following  Designs 

Figure  4.30  shows  the  Hoh's  Bandwidth  handling  qualities  predictions  for  the  angle  of 
attack  following  designs.  These  will  be  discussed  as  comments  are  made  on  each  design. 
However,  before  the  simulation  results  are  presented,  and  each  design  is  discussed,  it  is 
necessary  to  discuss  the  structure  of  the  filters. 

The  washout  filter  removes  the  integral  action  on  pitch  rate  and  pitch  angle  by  placing  a 
zero  at  the  origin.  It  also  gives  the  pilot  the  ability  to  use  monotonic  stick  forces  by  placing  a 
real  pole  at  typical  phugoid  magnitudes.  The  transfer  function  used  the  washout  filter  in  Figure 
4.17  is  given  by 
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(4.1) 


wqfilt  — 


20s 
205  +  1 


The  lead  filter  restores  pitch  rate  overshoot  which  results  in  a  faster  angle  of  attack 
response.  This  is  desirable  since  this  an  angle  of  attack  following  system.  The  transfer  function 
of  the  lead  filter  in  Figure  4.17  is 


lefilt  = 


1.1925+1 
.6595  +  1 


(4.2) 


The  values  for  the  zero  and  pole  in  equation  4.2  move  a  pole  in  the  q/qc  transfer  function  from 
-0.839  to  -1.517  which  produces  the  characteristics  mentioned  above. 
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Figure  4.18.  Angle  of  Attack  Simulations  for  Case  1-1  -  noises  off 


Figure  4.21.  Angle  of  Attack  Simulations  for  Case  1-2  -  Noises  on 


Figure  4.22.  Angle  of  Attack  Simulations  for  Case  1-3  -  noise  off 
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Figure  4.23.  Angle  of  Attack  Simulations  for  Case  1-3  -  noise  on 


10  15  20  25  30 

time  (seconds) 


0  5  10  15 

time  (seconds) 


20 


Figure  4.26.  Angle  of  Attack  Simulations  for  Case  1-5  -  noise  off 


Figure  4.  27.  Angle  of  Attack  Simulations  for  Case  1-5  -  noise  on 
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Figure  4.28.  Angle  of  Attack  Simulations  for  Case  l-5a  -  noise  off 


Figure  4.30.  Angle  of  Attack  Handling  Qualities  Predictions  for  Cases  1-1  through  1-5 

The  simulation  histories  are  not  as  straight  forward  to  evaluate  as  in  the  pitch  rate 
designs.  This  is  because  we  have  removed  the  integral  action  on  the  command,  and  therefore, 
the  controller  is  not  tracking  the  command  as  well.  The  important  factor  in  these  designs  is  that 
they  restore  conventional  response  characteristics  to  the  system  which  should  result  in  good 
handling  qualities  ratings.  However,  Figure  4.30  does  not  support  this  hypothesis  very  well  as 
compared  to  the  results  for  the  pitch  rate  designs  in  Figure  4.16.  In  general,  the  washout  and 
lead  filters  have  increased  the  bandwidth  of  all  of  the  design  cases.  But,  it  is  important  to 
remember  that  there  is  evidence  that  the  Hoh’s  Bandwidth  handling  qualities  prediction  criteria 
do  not  correlate  very  well  with  flight  test  results,  especially  in  the  case  of  the  angle  of  attack 
designs.  However,  a  criteria  has  been  proposed  by  Berry  [24]  that  shows  excellent  results.  It  is 
based  on  the  percentage  the  flight  path  overshoots  when  the  command  is  terminated.  Flight  test 
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data  indicates  that  a  flight  path  overshoot  of  less  than  20%  would  be  a  good  metric  to  predict 
level  1  pilot  ratings.  All  of  the  cases  in  Figures  4.18-4.29  have  between  a  14-20%  flight  path 
overshoot.  Case  1-1  is  the  worst  while  cases  1-3, 4,  and  5  are  the  best.  It  is  important  to  note 
that  the  flight  path  overshoot  parameter  can  easily  be  adjusted  by  changing  the  numbers  in  the 
washout  and  lead  filters.  However,  that  will  be  left  to  Phase  II. 

4.2.2  Angle  of  Attack  Designs. 

An  angle  of  attack  design  (case  1-6)  and  a  mixed  angle  of  attack 

design  (case  1-7)  were  developed.  The  only  two  changes  required  to  convert  from  a  pitch  rate 
problem  formulation  to  an  angle  of  attack  problem  formulation  were  to  change  the  command  and 
the  ideal  model  we  are  attempting  to  follow.  Changing  the  command  does  not  involve  changing 
the  problem  formulation  as  the  stick  dynamics  are  the  same  for  both  pitch  rate  and  angle  of 
attack  commands.  The  ideal  model  for  these  designs  is  given  by 


Wm  = 


6.25 

+  "iSs  +  6.25 


(4.3) 


This  model  is  second  order  to  include  some  pitch  rate  overshoot  and  therefore  eliminate  the  need 
for  the  lead  and  washout  filter  in  the  simulation.  The  only  other  change  made  to  the  original 
pitch  rate  problem  formulations  was  to  tighten  the  control  rate  usage.  This  is  given  by 


Wcr  = 


(4.4) 


This  was  required  to  prevent  destabilizing  rate  limiting  of  the  system.  The  stability  margin 
design  results  are  listed  in  Table  4.4.  Norm  data  is  not  included  because  the  primary 
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performance  metric  other  than  stability  margins  will  be  the  flight  path  angle  overshoot  from  the 
simulations. 


TABLE  4.4 


ANGLE  OF  ATTACK  DESIGN  DATA  -  STABILITY  MARGIN  COMPARISON 


Case 

Complement. 
Sensitivity 
Vector  Gain 
Margins  (dB) 
Input  of  Plant 

Sensitivity 
Vector  Gain 
Margins  (dB) 
Input  of  Plant 

Phase 

Margins  (deg) 
Input  of  Plant 

Complement. 
Sensitivity 
Vector  Gain 
Margins  (dB) 
Output  of 
Plant 

Sensitivity 
Vector  Gain 
Margins  (dB) 
Output  of 
Plant 

Phase 

Margins  (deg) 
Output  of 
Plant 

1-6 

[-2.81,2.12] 

[-2.28,  3.10] 

17.28 

[-2.76,2.09] 

[-2.26, 3.07] 

1-7 

[-15.61,5.27] 

[-5.57,  19.98] 

53.47 

[-11.74, 4.82] 

[-4.95,  12.70] 

45.19  1 

The  stability  margins  for  the  problem  are  unacceptable  at  both  the  input  and  output 
of  the  plant.  The  controller  for  this  case  also  has  a  RHP  pole.  However,  by  mixing  this  design 
with  a  tracking  constraint  and  an  input  uncertainty  constraint  )  >  we  improve  the 

stability  margins  tremendously  and  get  a  controller  with  all  LHP  poles.  The  simulations  for 
these  two  designs  will  be  presented  next  followed  by  a  brief  discussion  of  each  case's  handling 
qualities. 

4.2.3  Angle  of  Attack  Design  Simulations. 

The  simulations  are  conducted  as  before  with  the  only  difference  being  that  we  have  no 
lead  filter  as  shown  in  Figure  4.30  because  these  characteristics  are  accounted  for  in  the  ideal 
model.  The  washout  filter  is  still  used  and  is  the  same  as  equation  4.1 . 
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Figure  4.3 1 .  Angle  of  Attack  Simulations  for  Case  1-6  -  noise  off 
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Figure  4.34.  Angle  of  Attack  Simulations  for  Case  1-7  -  noise  on 
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The  biggest  difference  between  the  simulations  in  Figures  4.31  through  4.34  is  that  there 
is  much  less  overshoot  of  the  flight  path  angle  after  the  command  is  terminated  than  in  the 
previous  simulations  (Figures  4.18-4.29).  Both  cases  here  exhibit  about  7%  overshoot  while  the 
minimum  before  was  12%.  This  should  result  in  better  handling  qualities  ratings  if  the  flight 
path  overshoot  hypothesis  is  correct.  We  have  gained  a  fairly  significant  improvement  here  in 
our  main  handling  qualities  metric  by  designing  an  angle  of  attack  system  rather  than  converting 
a  pitch  rate  design. .  The  Hoh’s  Bandwidth  analysis  of  Case  1-7  predicted  Level  1  handling 
qualities  (bandwidth=2.9,  tau  p=.05  sec.,  where  tau  p  is  a  measure  of  delay  in  the  closed  loop 
response.  See  [15]  for  details). 


4.3  Summary 

Because  of  the  many  different  flight  control  designs  covered  in  this  chapter,  a  summary 
of  those  that  were  finally  selected  for  phase  II  is  provided  in  Table  4.5. 


TABLE  4.5 


SUMMARY  OF  FLIGHT  CONTROL  DESIGNS  SELECTED  FOR  FLIGHT  TEST 


Case 

Method  of 
Design 

Type  of  Design 

Hoh’s  Bandwidth 
HQ  Prediction 

Berry  Flight  Path 
Overshoot  Criteria  HQ 
Prediction 

1-1 

H2 

Pitch  Rate 
Cmd 

Level  1 

Level  1 

1- 

i(f)i 

Hi 

AoA  Cmd 

Level  1 

Level  1 

1-4 

H2/H^ 

Pitch  Rate 
Cmd 

Level  1 

Level  1 

1- 

4(f)l 

H2/H„ 

AoA  Cmd 

Level  1 

Level  1 

1-7 

H2/H„ 

AoA  Cmd 

Level  1 

Level  1 

(f)  -  denotes  the  pitch  rate  designs  that  were  changed  to  AoA  by  command  path  filters 
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All  of  designs  selected  for  flight  test  had  the  best  characteristics  in  terms  of  having  all 
stable  poles,  handling  qualities  predictions,  and  command  tracking.  The  H2  design  was  included 
along  with  the  mixed  norm  designs  to  test  the  following  hypothesis;  Given  an  accurate  aircraft 
model  (which  was  the  assumption  made  with  the  LearF16),  a  simple  LQG  design  with  good 
turbulence/noise  rejection  and  command  tracking  characteristics  would  be  sufficient  (eliminating 
the  need  to  spend  the  extra  time  and  work  to  develop  a  mixed  norm  design).  The  accuracy  of  the 
findings  and  hypotheses  stated  in  these  first  four  chapters  will  be  evident  in  the  flight  test  phase 
which  will  begin  with  the  flight  test  setup  in  the  next  chapter. 
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V.  Flight  Test  Preliminaries  : 

This  chapter  begins  the  discussion  of  the  flight  test  phase  of  this  research.  Areas 
covered  will  include;  a  brief  discussion  of  the  logistics  in  setting  up  the  flight  test,  methodology, 
procedures  of  the  flight  test,  implementation  and  verification,  and  finally,  ground  and  airborne 
validation  of  the  control  laws.  Three  new  control  laws  designs  that  were  generated  based  on 
implementation  and  verification  results  are  discussed  in  the  last  two  sections. 

5.1  Flight  Test  Logistics 

All  flight  test  programs  are  required  to  have  a  program  name  that  can  be  referred  to  in 
the  documentation  required  to  set  up,  execute,  and  report  on  the  program.  This  research  was 
assigned  the  name  HAVE  INFINITY  II  by  the  U.S.  Air  Force  Test  Pilot  School.  This  point  is 
made  to  source  the  references  to  HAVE  INFINITY  II  in  this  chapter. 

The  initial  work  in  setting  up  this  test  program  dealt  with  securing  the  resources, 
facilities,  and  aircraft  required  to  perform  the  flight  test  program.  The  HAVE  INFINITY  II  test 
team  (3  pilots,  1  EWO,  and  2  Engineers)  was  responsible  (through  the  USAF  TPS  Commandant) 
for  all  phases  of  test  program  management  and  execution. 

5. 1. 1  Flight  Test  Resources. 

The  HAVE  INFINITY  II  program  was  funded  by  USAF  TPS.  The  total  program  budget 
was  $109,000.  The  major  items  funded  in  the  program  budget  included:  a  two  day  trip  to 
Calspan  facilities  in  Buffalo,  NY;  nine  1 .0  hour  support  sorties  (3  C-23  and  6  T-38)  for  the  flight 
test  program;  and  twelve  1.2  hour  sorties  in  the  Calspan  Variable  Stability  (VSS)  Learjet  24  for 
the  flight  test  program. 
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5. 1.2  Flight  Test  Facilities  and  Procedures 


Model  verification  and  validation  testing  was  performed  at  Calspan  facilities  in  Buffalo, 
NY  and  at  Edwards  AFB,  CA.  All  test  program  landings  were  accomplished  at  Palmdale,  CA 
Regional  Airport  under  a  Memorandum  of  Understanding  between  USAF  TPS  and  AF  Plant  42. 
Landings  were  performed  at  Palmdale  to  avoid  the  congestion  of  the  Edwards  AFB  traffic 
pattern.  The  logistics  of  getting  the  runway  marked  for  scoring  test  program  landings  were  also 
much  easier  at  Palmdale. 

5.1.3  Support  and  Flight  Test  Aircraft. 

The  six  T-38  support  sorties  were  used  to  practice  and  standardize  the  landing  tasks  that 
would  later  be  performed  in  the  Calspan  Learjet.  The  C-23  was  used  as  a  target  aircraft  for  the 
Calspan  Learjet  to  evaluate  any  PIO  tendencies  of  the  HAVE  INFINITY  II  flight  control 
systems  during  various  tracking  tasks.  The  Calspan  VSS  Learjet  served  as  the  test  aircraft.  The 
VSS  Learjet  is  a  highly  modified  Learjet  Model  24  that  fiinctions  as  a  three  axis  in-flight 
simulator.  The  cockpit  has  two  sets  of  side  by  side  controls.  The  controls  at  the  left  seat,  for  the 
safety  pilot,  maintain  the  Learjet’s  conventional  flight  control  system.  The  evaluation  pilot 
occupies  the  right  seat.  Those  controls  use  a  fly-by-wire,  response  feedback,  variable  stability 
and  variable  control  system,  consisting  of:  variable  feel  system,  aircraft  motion  sensors  and 
associated  signal  conditioning,  control  system  simulation  computer,  control  surface  servos, 
digital  configuration  control  system,  engage/disengage  and  safety  monitor  logic,  and 
recording/playback  capability. 
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5.2  Flight  Test  Methodology 


Specific  methods  and  procedures  had  to  be  established  in  the  planning  phase  to 
effectively  execute  the  HAVE  INFINITY  II  flight  test  program.  This  section  starts  by  detailing 
the  methodology  that  was  used  in  the  program.  The  second  section  covers  specific  flight  test 
procedures  that  support  the  methodology. 

5.2. 1  Flight  Test  Methodology. 

The  flight  iest  methodology  can  be  thought  of  as  a  road  map  that  specifies  what  the  test 
team  is  going  to  do  to  successfully  complete  the  flight  test.  It  starts  with  an  overall  purpose  and 
is  followed  by:  specific  test  objectives,  Measures  of  Performance  (MOP’s)  that  outline  what  is 
required  to  meet  an  objective,  success  criteria,  evaluation  criteria,  and  test  methodology. 

Purpose.  The  purpose  of  the  HAVE  INFINITY  II  limited  flight  test  was  to  evaluate  the 
handling  qualities  of  H2,  mixed  H2/H^,  and  classical  longitudinal  flight  control  system  designs 
during  the  approach  and  landing  phase  of  flight. 

The  classical  flight  control  system  will  be  discussed  in  detail  in  section  5.4.  The  test  was 
conducted  at  Edwards  AFB,  California  from  1  October  1997  to  9  October  1997.  The  overall  test 
point  matrix  is  shown  in  Appendix  C.  The  evaluation  order  of  the  various  flight  control 
configuration  was  randomized  by  the  TC  and  referred  by  the  code  letter  in  Table  Cl  to  preserve 
the  integrity  of  evaluation  amongst  the  test  team  pilots. 

The  following  paragraphs  outline  the  HAVE  INFINITY  II  specific  test  objectives, 

MOPs,  success  criteria,  evaluation  criteria,  and  test  methodology  that  was  used  to  successfully 
complete  this  test  program. 
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Test  Objective  1  -  Flight  Control  Configuration  Evaluatidns.  Evaluate  the  longitudinal 
handling  qualities  of  H2,  mixed  H2/H„  and  classical  longitudinal  flight  control  configurations. 
Measures  of  Performance  (MOP)  were: 

1 .  Pilot-in-the-loop  Oscillation  (PIO)  ratings  for  Handling  Qualities  During  Tracking 

(HQDT). 

2.  Cooper-Harper  (C-H)  ratings,  PIO  ratings  and  pilot  comments  for  approach  to  main 

wheel  touchdown  (MWTD). 

MOP  1  -  PIO  Ratings  for  Handling  Qualities  During  Tracking  (HQDT).  Susceptibility 
to  PIO  of  the  HAVE  INFINITY  II  flight  control  configurations  was  evaluated  by  performing 
high  bandwidth  tracking  on  an  airborne  target.  A  build-up  approach  was  accomplished  by 
performing  low  bandwidth  tracking  prior  to  high  bandwidth  tracking. 

Success  Criteria.  Collect  PIO  ratings  from  3  pilots  for  each  flight  control  configuration. 

Evaluation  Criteria.  Each  flight  control  configuration  was  evaluated  against  the  PIO 
rating  scale  (Appendix  A).  Ratings  of  1, 2,  and  3  were  considered  satisfactory.  A  rating  of  4 
was  considered  unsatisfactory,  but  tolerable.  A  rating  of  5  or  6  was  considered  unacceptable. 

Test  Methodology.  This  test  was  flown  with  a  C-23  airborne  target.  The  Learjet  was  in 
the  power  approach  configuration.  Each  flight  control  configuration  was  evaluated  using  a 
buildup  approach  starting  with  low  bandwidth  tracking  from  a  position  1 000  feet  in  trail  of  the 
target  and  progressing  to  high  bandwidth  HQDT.  Low  bandwidth  tracking  is  non-aggressive, 
gentle  maneuvering  and  open-loop  pulses  and  steps.  Handling  qualities  during  tracking  (HQDT) 
is  aggressive  and  assiduous  tracking  of  a  precision  aimpoint  to  zero  error.  During  both  low 
bandwidth  and  high  bandwidth  tracking,  the  evaluation  pilot  attempted  to  change  the  desired  aim 
point  by  10  mils.  If  the  flight  control  configuration  did  not  exhibit  divergent  oscillations  during 
normal  control  inputs  (PIO  6)  or  divergent  oscillations  during  tight  control  (PIO  5),  close 
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formation  was  performed.  In  close  formation,  the  evaluation  pilot  started  with  low  bandwidth 
tracking  and  proceeded  to  HQDT  in  a  buildup  fashion.  During  close  formation  HQDT,  the 
evaluation  pilot  attempted  to  correct  the  aircraft  to  a  desired  position  from  10  feet  below  the 

desired  position.  Close  formation  was  accomplished  on  a  30°  line  (line  up  the  main  wheel  of  the 
C-23  with  the  front  antenna)  with  10  ft  wingtip  clearance  and  nose-tail  separation  with  the  C-23. 
A  separate  PIO  rating  was  given  for  the  trail  position  and  for  the  close  formation  position,  if 
accomplished.  Landing  tasks  were  not  attempted  with  any  flight  control  configuration  assigned 
a  PIO  rating  of  5  of  6  during  HQDT.  The  Gibson  and  Ralph  Smith  PIO  criterion  did  not  predict 
PIO  for  any  of  the  HAVE  INFINITY  II  flight  test  candidates  designed  during  phase  I. 

MOP  2  -  Cooper-Harper  (C-H)  ratings,  PIO  ratings  and  pilot  comments  for  approach 
to  main  wheel  touchdown  (MWTD).  The  handling  qualities  of  the  HAVE  INFINITY  II  flight 
control  configurations  were  evaluated  by  performing  a  series  of  landings  with  each 
configuration.  Each  configuration  started  with  a  straight-in  approach  to  a  spot  landing  and 
proceeded  to  offset  landing  tasks  with  spot  landings  if  controllability  was  not  in  question. 

Success  Criteria.  Collect  C-H  ratings,  PIO  ratings  and  pilot  comments  from  three  pilots 
for  each  flight  control  configuration  that  did  not  receive  a  PIO  rating  of  5  or  6  during  MOP  I . 
The  number  of  ratings  given  by  a  pilot  for  each  configuration  was  determined  by  the  Landing 
Task  Decision  Tree,  Figure  5.1. 
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Figure  5.1 .  Landing  Task  Decision  Tree 

Evaluation  Criteria.  The  handling  qualities  of  each  flight  control  configuration  were 
categorized  into  one  of  three  levels  dependent  on  landing  performance  and  pilot  workload. 
Handling  qualities  were  rated  using  the  Cooper-Harper  rating  scale  contained  in  Appendix  A. 
Level  I  was  defined  as  satisfactory,  no  improvement  is  necessary.  Cooper-Harper  ratings  of  1,  2 
and  3  are  Level  I.  Level  II  was  unsatisfactory,  but  tolerable,  deficiencies  warrant  improvement. 
Cooper-Harper  ratings  of  4,  5  and  6  are  Level  II.  Level  III  was  unacceptable,  deficiencies 
require  improvement.  Cooper-Harper  ratings  of  7,  8  and  9  are  Level  III.  Histograms  of  C-H 
ratings  were  analyzed  along  with  PIO  ratings  and  pilot  comments.  An  appropriate  handling 
qualities  level  was  assigned  if  the  C-H  ratings  are  grouped  within  one  level.  If  a  large  dispersion 
in  ratings  was  encountered  or  the  ratings  were  split  between  levels,  it  was  not  possible  to  assign  a 
handling  qualities  level.  In  these  cases  pilot  comments  and  engineering  judgment  were  used  to 
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determine  an  appropriate  handling  qualities  level  if  possible.  PIO  ratings  were  evaluated  as 
described  in  MOP  1  Test  Methodology. 

Test  Methodology.  The  goal  of  this  test  was  to  perform  one  straight-in  landing  and  a 
horizontal  offset  landing  for  each  flight  control  configuration.  If  controllability  was  in  question 
(C-H  rating  of  8,  9,  or  1 0)  on  any  landing,  testing  of  that  configuration  by  that  pilot  was 
terminated.  If  lateral-directional  dynamics  effect  performance  during  the  horizontal  offset 
landing,  the  pilot  performed  a  vertical  offset  landing. 

This  test  started  by  flying  a  straight-in  approach  to  a  spot  landing  using  the  3°  Instrument 
Landing  System  (ILS)  glideslope  as  described  in  the  next  section.  The  evaluation  pilot  assigned 

C-H  and  PIO  ratings  and  made  comments  addressing  longitudinal  handling  qualities.  If  the 
flight  control  configuration  was  not  assigned  a  C-H  rating  of  8, 9,  or  10,  a  horizontal  offset 
landing  was  performed.  This  test  technique  is  described  in  the  next  section.  If  the  pilot  felt  the 
lateral-directional  handling  qualities  effected  the  C-H  rating  for  the  horizontal  offset  landing,  a 
vertical  offset  landing  was  performed.  This  test  technique  is  also  described  in  the  next  section. 
The  vertical  offset  was  only  performed  if  the  pilot  did  not  assign  a  C-H  rating  of  8,  9,  or  1 0  for 
the  horizontal  offset,  and  the  lateral-directional  handling  qualities  effect  the  C-H  rating.  This 
process  was  completed  for  each  flight  control  configuration. 

Test  Objective  2  -  Pilot  Comments  After  Main  Wheel  Touchdown.  Collect  pilot 
comments  on  the  longitudinal  handling  qualities  of  the  flight  control  configuration  after  main 
wheel  touchdown. 

MOP  1.  Collect  pilot  comments  on  the  longitudinal  handling  qualities  of  the  flight 
control  configuration  after  main  wheel  touchdown. 


5-7 


Success  Criteria.  Collect  pilot  comments  from  all  landings. 

Evaluation  Criteria.  None. 

Test  Methodology.  After  MWTD,  the  aircraft  stability  derivatives  are  drastically 
different  than  the  design  point  of  125  KJAS  and  1,000  feet  PA.  Therefore,  this  portion  of  the 
landing  task  was  evaluated  separately.  The  evaluation  pilot  lowered  the  nose  wheel  to  the 
runway  after  main  wheel  touch  down.  The  evaluation  continued  until  the  pilot  was  assured  a  full 
stop  landing  could  be  accomplished.  The  pilot  made  qualitative  comments  on  the  longitudinal 
handling  qualities  of  the  flight  control  configuration  during  this  phase. 

This  objective  was  added  due  to  HAVE  INFINITY’S  experience  with  an  uncommanded 
pitch  up  in  the  landing  flare  for  the  optimal  flight  control  laws. 

5.2.2  Flight  Test  Procedures 

General.  Ground  testing  was  accomplished  on  every  flight  control  configuration  prior  to 
in-flight  evaluation.  Ground  testing  provided  a  means  of  software  implementation  verification, 
of  the  designed  configuration,  on  the  Learjet  VSS. 

Model  validation  data  of  the  aircraft  design  model  and  each  flight  control  configuration 
was  collected  at  140  KIAS  with  the  aircraft  configured  as  in  the  landing  task  (landing  gear  down 
and  flaps  20%).  Using  Programmed  Test  Inputs,  at  least  three  pitch  doublets  and  two  frequency 
sweeps  were  performed  for  each  configuration. 

Each  evaluation  pilot  accomplished  high  bandwidth  HQDT  on  every  flight  control 
configuration,  in  both  1 000  ft  trail  and  close  formation,  prior  to  advancing  to  the  landing  task. 
HQDT  was  accomplished  at  140  KIAS  and  at  a  minimum  altitude  of  5,000  ft  AGL.  The  airspeed 
parameter  was  established  to  ensure  quality  of  data  between  the  HQDT  phase  accomplished  up- 
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and-away  and  the  actual  landing  phase.  Each  flight  control  configuration  was  evaluated  using  a 
buildup  approach  starting  with  low  bandwidth  tracking,  from  a  position  1000  feet  in  trail  of  the 
target,  and  progressing  to  high  bandwidth  tracking  (HQDT)  as  described  in  the  previous  section. 

The  evaluation  pilot  performed  a  minimum  of  one  approach  and  landing  with  the  VSS 
configured  with  the  baseline  Learjet  flight  control  system.  This  was  accomplished  during  each 
flight  prior  to  any  actual  testing  with  the  HAVE  D'lFINITY  II  flight  control  configurations 
installed. 

The  test  conductor  configured  the  Learjet  flight  control  system  with  the  required 
configuration  parameters,  and  the  safety  pilot  engaged  the  Variable  Stability  System.  The 
evaluation  pilot  took  control  of  the  aircraft  for  the  pattern  and  landing.  The  safety  pilot  took 
control  of  the  aircraft  upon  the  evaluation  pilot’s  determination  that  a  full  stop  landing  could 
have  been  accomplished  from  the  current  landing.  The  safety  pilot  then  took  control  for  the 
remainder  of  the  ground  roll,  take  off,  and  the  climb  to  pattern  altitude.  The  evaluation  pilot 
provided  comments  on  atmospheric  conditions  which  affected  the  approach. 

Straight-in  Landing.  Evaluation  pilots  flew  a  3°  glide  path,  using  the  ILS  glideslope  to 
assist  with  glide  path  determination,  until  the  decision  height  of  200  ft  AGL,  then  transitioned 
visually  to  consistently  flare  and  touchdown  in  the  desired  zone.  The  approach  airspeed  was 
between  125  and  135  KIAS  depending  on  the  aircraft  weight.  The  touchdown  point  should  have 
been  in  the  desired  zone  described  in  Figure  5.2.  The  target  airspeed  at  touchdown  was  10  kts 
less  than  approach  speed  (+10/-5  kts).  For  quality  of  data,  only  landings  within  this  touchdown 
airspeed  window  were  evaluated.  The  ground  test  team  made  a  radio  call  to  notify  the  aircrew  of 
the  actual  touchdown  point,  and  if  desired  or  adequate  performance  was  attained.  The  evaluation 
pilot  provided  qualitative  comments,  Cooper-Harper  and  PIO  ratings  for  the  approach  to  main 
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wheel  touchdown  in  accordance  with  Appendix  A.  After  main  wheel  touchdown,  the  evaluation 
pilot  provided  qualitative  comments  through  nosewheel  touchdown  and  to  a  point,  determined  by 
the  evaluation  pilot,  that  a  full  stop  landing  could  be  completed.  In  accordance  with  Figure  5.1, 
if  the  straight-in  landing  received  a  Cooper-Harper  rating  of  other  than  8,  9,  or  10,  this  flight 
control  configuration  would  be  evaluated  with  a  Horizontal  Offset  Landing. 

Horizontal  Offset  Landing.  This  task  forced  the  pilot  to  raise  his  gains  by  having  to 
concentrate  on  both  longitudinal  and  lateral  control  inputs.  The  Horizontal  Offset  Landing  task 
was  accomplished  by  flying  a  visual  pattern  with  a  lateral  offset  of  300  ft  from  runway 
centerline.  At  200  ft  AGL,  the  evaluation  pilot  aggressively  corrected  to  the  centerline.  The 
correction  used  an  initial  bank  angle  between  30  and  45  degrees  within  three  seconds,  and  the 
initial  aggressive  lateral  corrections  were  completed  by  100  ft  AGL.  A  simultaneous  correction 
was  made  to  intercept  a  visual  glide  path  to  touchdown.  The  desired  touchdown  point  was  at  the 
center  of  the  desired  zone.  Again,  the  target  airspeed  at  touchdown  was  10  kts  less  than 
approach  speed  (+10/-5  kts),  and  only  landings  within  this  touchdown  airspeed  window  were 
evaluated.  Pilot  comments  and  ratings  were  identical  to  the  straight-in  landing.  In  accordance 
with  Figure  5.1,  if  the  Horizontal  Offset  Landing  received  a  Cooper-Harper  rating  of  other  than 
8,  9,  or  1 0  and  the  lateral-directional  handling  qualities  affected  the  Cooper  Harper  rating,  this 
flight  control  configuration  could  be  evaluated  with  a  Vertical  Offset  Landing. 

Vertical  Offset  Landing.  This  task  was  designed  to  minimize  any  lateral-directional 
handling  qualities  effects  and  focus  attention  on  the  longitudinal  axis.  The  Vertical  Offset 
Landing  was  accomplished  by  flying  straight  and  level  at  the  published  Minimum  Descent 
Altitude  (MDA  -  397’  AGL),  as  if  flying  a  localizer  approach.  At  glideslope  intercept  (Visual 
Descent  Point),  the  pilot  would  aggressively  correct  to  a  3°  glide  path  by  100  ft  AGL. 
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Comments  and  ratings  were  identical  to  those  discussed  in  the  straight-in  landing.  Only  landings 
within  the  touchdown  airspeed  window  were  evaluated. 

A  landing  could  be  repeated,  at  the  discretion  of  the  TC,  based  on  improper  setup, 
extenuating  atmospheric  conditions  (windshear,  etc.),  or  any  factor  where  a  biased  rating  may 
have  occurred.  Test  conditions  such  as  minor  turbulence  or  wind  variations  were  documented, 
but  not  repeated. 

Excessive  Tailwind  Procedures  or  Runway  Non-Availability.  If  the  tailwind  component 
for  runway  25  at  Palmdale  Airport  was  greater  than  10  kts,  runway  07  was  used  for  landing.  The 
precision  instrument  runway  markings  depicted  in  Figure  5.2  were  used  to  help  the  pilot  define 
the  desired  and  adequate  landing  zones.  Ground  spotters  were  positioned  1200  feet  and  1600 
feet  from  the  runway  threshold  to  verify  landings  in  the  desired  zone.  If  the  Palmdale  runway 
was  closed  or  unusable,  any  other  runway  in  the  local  flying  area  could  have  been  used  with  a 
published  instrument  precision  approach  using  the  same  procedures  described  for  runway  07. 


Figure  5.2.  Touchdown  Zone 


5.3  Control  Law  Implementation  and  Verification 

The  HAVE  INFINITY  II  control  laws  were  implemented  by  converting  the  state  space 
representation  to  a  transfer  function  format  for  each  channel  of  the  control  laws  (which  is  the 
standard  implementation  for  Calspan).  This  implementation  differed  from  the  HAVE  INFINITY 
implementation.  The  HAVE  INFINITY  program  was  forced  to  use  a  state  space  implementation 
due  to  instabilities  in  their  control  laws.  The  Calspan  VSS  computer  used  MATLAB™  and 
SIMULINK™  to  implement  and  simulate  the  HAVE  INFINITY  II  control  configurations.  All  of 
the  control  laws  had  three  input  channels  (command,  along  with  q  and  a  feedback)  and  one 
output  channel  (elevator  command).  This  is  shown  in  the  SIMULINK™  block  diagram  in 
Figure  5.3.  The  command  channel  implementation  model  follows  in  Figure  5.4. 
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Figure  5.3.  Three  Channel  Control  Law  Implementation  Model 


Figure  5.4.  Command  TF  Channel  Implementation  Model 
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The  q  and  a  transfer  fiinction  implementation  models  are  very  similar.  The  structure  of 
the  implementation  models  did  not  change  between  HAVE  INFINITY  II  flight  control 
configurations.  Differences  in  controller  order  and  structure  of  the  dynamics  of  the  pole/zeros 
between  control  configurations  were  accounted  for  by  activating/deactivating  the  necessary 
transfer  function  loops  within  each  channel.  This  was  done  with  the  software  switches  that  are 
apparent  in  Figure  5.4.  Initial  HAVE  INFINITY  II  control  law  verification  and  simulation  was 
accomplished  when  integration  and  debugging  of  the  control  implementation  architecture,  with 
the  rest  of  the  Learjet  VSS,  was  completed.  The  first  verification  work  dealt  with  making  sure 
the  aircraft  model  was  correctly  implemented.  Verification  of  the  aircraft  poles  and  zeros  is  in 
Table  Dl,  Appendix  D.  The  implementation  transfer  functions  for  each  of  the  control  laws  from 
Table  4.5  are  shown  in  Tables  5. 1-5.3.  The  format  of  the  transfer  functions  is  described  in  the 
first  paragraph  of  Appendix  D. 


TABLE  5.1 

COMMAND  CHANNEL  IMPLEMENTATION  TRANSFER  FUNCTION 


Case 

Command  Channel 

1-1 

-  229 1.2359[0.7000,60.000](2 1.4843)(1.5 1 7 1)(0.2990) 

[0.7 1 89,277.1 040]  [0.9999,3 1.6083]  [6.0000)(0.8405) 

1-4 

- 1 693.5630[0.836 1,73.6640](1 8.3649)[0.7060,1.0943] 
[0.7250,276.6304](35.2238)(32.6059)[0.3659,1.4619] 

1-5 

-  36.8823(53.0559)[0.7064,24.8286](1.3 192)(0.8195)(0.4175) 

(85.4428)[0.6456,25.3192][0.9375,5.1344][0.9454,0.3811] 
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TABLE  5.2  : 

PITCH  RATE  CHANNEL  IMPLEMENTATION  TRANSFER  FUNCTION 


Case 

q  Feedback  Channel 

1-1 

1252.4538[0.7000,60.0000](6.0000)(4.0809)(1.4402) 

[0.7 1 89,277.1 040]  [0.9999,3 1.6083](6.0000)(0.8405) 

1-4 

1223.7806[0.9171,89.9112](16.5553)[0.7593, 0.9092] 
[0.7250, 276.6304](35.2238X32.6059)[03659,1.4619] 

1-5 

43.444 1[0.7898,34.0126][0.8267,5.3258](0.8936)(0.4391) 
(85.4428)[0.6456,25.3 1 92]  [0.93753.1344][0.9454,0.38 1 1] 

TABLE  5.3 

AOA  CHANNEL  IMPLEMENTATION  TRANSFER  FUNCTION 


Case 

a  Feedback  Channel 

1-1 

10.6641[0.7000,60.000](6.0000)(4.0809)(1.4402) 

[0.7 1 89,277.1 040]  [0.9999,3 1.6083](6.0000)(0.8405) 

1-4 

147.8755(168.6614)(37.9513)(22.9369)[0.7248, 1.9879] 

[0.7250,276.6304](35.2238X32.6059)[03659,1.4619] 

1-5 

33357(603737)[0.6775,24i839](4.9107)[0.9705,0.4375] 

(85.4428)[0.6456,25.3192][0.9375,5.1344][0.9454,0.3811] 

Recall  that  case  1-1  and  1-4  also  had  AOA  versions.  This  was  done  by  placing  the  filters 
described  in  chapter  IV  in  the  command  path  prior  to  the  controller.  A  software  switch  was  also 
placed  in  this  path  so  these  filters  could  be  bypassed  when  the  pitch  rate  versions  and  case  1-5 
were  being  tested.  Initial  verification  and  simulation  took  place  at  Calspan  facilities  in  Buffalo, 
NY  from  3  to  5  August  1997.  The  HAVE  INFINITY  II  test  team  felt  it  was  important  to  be  able 
to  address  any  implementation  problems  early  in  the  program.  This  would  allow  time  to  make 
any  necessary  design  adjustments  prior  to  flight  test.  It  turns  out  that  both  case  1-1  and  case  1-4 
had  implementation  problems. 

The  problems  in  both  cases  dealt  with  high  frequency  modes.  Calspan  had 
recommended  using  Euler  Integration  with  a  0.01  sec  step  time  (100  Hz  sample  rate)  for  both 
desktop  and  airborne  SIMULINK™  simulations.  This  integration  routine  was  used  by  Calspan 
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because  of  its  speed  and  simplicity.  While  the  Euler  method  was  the  simplest  and  fastest 
integration  routine  offered  in  SIMULINK™,  it  was  also  the  most  inaccurate.  Cases  1-1  and  1-4 
(both  AoA  and  pitch  rate  versions)  were  numerically  unstable  when  simulated  with  Euler 
Integration.  However,  they  were  numerically  stable  when  simulated  using  the  ‘linsim’  and 
‘rk45’  integration  routines.  Upon  consultation  with  Calspan  engineers  during  phase  I,  they 
stated  that  the  ‘rk45’  integration  routine  had  been  used  successfully  in  the  Learjet.  It  was  agreed 
that  ‘rk45’  could  also  be  used  on  this  flight  test  program.  However,  when  the  control  laws  were 
implemented  on  the  Learjet,  it  was  discovered  that  ‘rk45’  could  not  be  used  as  an  integration 
routine.  It  introduced  too  much  computational  delay  into  the  system  and  caused  the  VSS  to 
continually  trip  off  when  any  control  inputs  were  made.  Consequently  the  decision  was  made  to 
use  Euler  Integration  and  try  to  reconfigure  cases  1-1  and  1-4  to  work  with  this  integration 
method. 

Case  1-1.  Reconfiguring  this  case  to  work  with  Euler  Integration  was  successful. 
Analysis  of  the  transfer  functions  for  each  channel  revealed  that  there  were  zeros  in  this  control 
law  that  were  canceling  the  actuator  poles.  There  was  also  a  set  of  very  high  frequency  poles  at 

277  radians/sec.  It  appeared  that  the  control  law  had  canceled  the  actuator  dynamics  and  put  in 
its  own  to  meet  the  tight  tracking  constraint  on  the  first  order  6/(s+6)  model.  This  hypothesis 
was  confirmed  by  pole-zero  canceling  these  high  frequency  dynamics  from  each  channel  and 
making  the  corresponding  gain  adjustment.  The  Euler  SIMULINK™  time  responses  were  now 
stable,  and  still  matched  very  closely  to  those  in  Figures  4.4/4. 5.  The  only  difference  noted  was 
that  the  rise  time  on  the  pitch  rate  response  increased  slightly.  This  was  considered  an 
improvement  as  the  Calspan  test  pilots  felt  the  original  response  would  have  been  rated  as  too 
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abrupt.  The  stability  margins,  and  handling  qualities  predictions  were  unchanged.  The  2-norm 
increased  from  2.518  to  5.1,  reflecting  the  “slower”  response. 

It  is  appropriate  to  note  here  that  there  was  of  some  confusion  amongst  HAVE 
INFINITY  II  team  members  and  Calspan  engineers  regarding  the  naming  conventions  used  for 
the  control  laws.  Cases  1-1/1 -1(f),  and  1-4/1 -4(f)  were  easily  confused  when  collecting 
verification  and  initial  simulation  data.  Hence,  the  names  were  changed  from  “cases”  to  the 
script  file  name  that  was  used  to  create  each  design.  The  new  names  are  shown  in  Table  5.4  and 
will  be  used  from  this  point  forward. 


TABLE  5.4 

NEW  HAVE  INFINITY  II  CONTROL  LAW  NAMING  CONVENTIONS 


Old  Name 

New  Name 

Case  1-1 

H2INI 

Case  l-l(f) 

H2AIN 

Case  1-4 

MXINI 

Case  l-4(f) 

MXAIN 

Case  1-5 

MXAOA 

Verification  of  the  implementation  of  the  “new”  poles  and  zeros  of  the  H2INI/H2AIN 
control  law  is  in  Table  D2,  Appendix  D.  Verification  of  the  implementation  of  the  MXAOA 
control  law  poles  and  zeros  is  in  Table  D3,  Appendix  D. 

MXINI/MXAIN.  Reconfiguring  this  transfer  function  was  not  successful.  The  high 
frequency  dynamics  that  caused  the  numerical  instability  with  Euler  Integration  could  not  be 
pole-zero  canceled  in  the  AOA  channel  without  making  this  transfer  function  improper  (see 
Table  5.3).  This  rendered  the  MXINI/MXAIN  control  laws  unfit  for  implementation  and 
disqualified  it  from  further  consideration. 
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At  this  point,  the  HAVE  INFINTY II  test  team  felt  it  had  time  to  replace  the  two  control 
configurations  lost  with  the  disqualification  of  MXINI/MXAIN,  in  addition  to  designing  a 
classical  feedback  gain  configuration.  The  classical  configuration  was  suggested  by  Calspan  as  a 
simple  design  that  would  have  known  level  I  handling  qualities.  This  design  could  then  be  used 
as  a  baseline  to  compare  the  handling  qualities  of  the  multiobjective  designs.  The  test  team  also 
decided  to  produce  an  H2  AOA  command  design  and  a  mixed  H2/H„  pitch  rate  command  design. 
The  H2  AOA  design  would  be  fast  to  produce  and  add  a  single  norm  AOA  command  design  to 
the  filtered  pitch  rate  configuration  (H2AIN).  The  mixed  H2/H„  pitch  rate  design  replaced 
MXINI.  MXAIN  was  not  replaced  because  MXAOA  already  represented  a  mixed  H2/H„  AOA 
command  design. 

Classic  Feedback  Gain  Design.  This  design  was  produced  using  only  simple  q  and  a 
feedback  gains.  The  aircraft  dynamics  The  pitch  rate  feedback  gains,  Kq,  =  0.441,  and  the  AOA 
feedback  gain,  K<x,  =  1.60,  were  provided  by  Calspan.  These  gains  were  known  to  produce  Level 
I  handling  qualities  for  the  Learjet  (simulating  an  F-16)  based  on  previous  flight  test  experience. 

This  control  design  was  given  the  name  CLASSIC.  H2  and  norms  were  not 
computed  for  CLASSIC.  Vector  margins  are  shown  in  Table  5.5. 


TABLE  5.5 

CLASSIC  STABILITY  MARGINS 


Name 

Complement. 
Sensitivity 
Vector  Gain 
Margins  (dB) 
Input  of  Plant 

Sensitivity 
Vector  Gain 
Margins  (dB) 
Input  of  Plant 

Phase 
Margins 
(deg) 
Input  of 
Plant 

Complement. 
Sensitivity 
Vector  Gain 
Margins  (dB) 
Output  of 
Plant 

Sensitivity 
Vector  Gain 
Margins  (dB) 
Output  of 
Plant 

Phase 

Margins  (deg) 
Output  of 
Plant 

CLASSIC 

[-00, 6.02] 

[-6.02,  C»] 

60.00 

[-5.02,  3.16] 

[-3.02, 4.66] 

25.37 

While  the  stability  margins  at  the  output  of  the  plant  are  lower  than  desired,  they  are  better  than 
the  H2INI/H2AIN  configurations.  This  design  had  also  been  proven  in  previous  flight  test  by 
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Calspan,  so  the  lower  stability  margins  were  considered  acceptable.  Time  histories  for  this 
design  are  shown  in  Figure  5.5  and  5.6. 
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Input  Command  Pulse 


Pitch  Rate,  5  deg,  4  sec  pulse 


Handling  qualities  predictions  for  this  design  were  calculated  using  the  Hoh’s  Bandwidth  and 
Berry  Flight  Path  Overshoot  Criteria.  The  Hoh’s  bandwidth  criteria  predicted  level  I  with  a 
delay  parameter  (tau  p)  of  0.02  and  a  bandwidth  of  2.9  rad/sec.  The  Berry  Criteria  also  predicted 
level  I  with  a  19%  overshoot. 


Implementing  the  CLASSIC  control  law  was  very  easy  since  no  dynamics  were  involved 
in  the  control  law.  Simple  gain  blocks  were  placed  at  the  appropriate  place  in  the  control  law 
path  of  the  VSS  software  architecture.  The  gains  were  then  simply  defined  as  the  values 
selected  when  needed,  and  zero  when  evaluating  a  different  control  law.  Verification  of  the 
implementation  of  the  CLASSIC  control  gains  is  in  Table  D4,  Appendix  D. 

H2  Angle  of  Attack  Command  Design.  This  design  was  produced  so  a  single  norm  AOA 
command  system  would  be  tested  in  addition  to  the  existing  pitch  rate  design  (H2rNI).  It  was 
hoped  that  some  handling  qualities  preferences  might  be  distinguished  between  the  AOA  and 
pitch  rate  designs.  The  problem  setup  for  this  AOA  design  was  the  same  as  described  in  section 
3.2,  except  that  the  output  variable  being  differenced  with  the  ideal  model  is  a  rather  than  q. 
Some  of  the  constraint  and  performance  weightings  have  also  been  changed  from  previous 
designs. 

The  ideal  model  used  in  this  design  was  chosen  to  have  the  same  dynamics  as  the 
CLASSIC  design  since  they  had  known  level  I  handling  qualities.  The  dynamics  of  this  ideal 
model  worked  out  to  be 


Wm  = 


5.29 

5^  +2.43J  +  5.29 


(5.1) 


Other  constraint  and  performance  weight  changes  include  the  tracking  performance 
weight,  Wp,  the  control  use  constraint,  Wcu,  and  the  control  rate  usage  constraint,  Wcr.  A  static 
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weight  of  20  was  chosen  for  Wp.  Recall  that  this  weighting  will  actually  be  the  inverse  of  the 
value  chosen.  Therefore  it  will  allow  up  to  a  5%  tracking  error  over  all  frequencies.  The  control 
rate  and  control  use  weights  were  set  equal  to  1 .  This  removed  any  control  use/rate  restrictions 
from  the  problem.  Since  the  ideal  model  is  a  second  order  system  with  relatively  slow  rise  time, 
rate  limiting  and  running  out  of  control  authority  was  not  a  problem.. 

The  plant,  actuator,  turbulence,  state  weighting,  and  sensor  noise  models  all  remained 
unchanged  from  previous  designs. 

This  control  design  was  given  the  name  H2AOA.  H2  and  H„  norms  are  shown  in  Table 
5.6  and  vector  stability  margins  are  shown  in  Table  5.7. 


TABLE  5.6 


NORM  DAI 

rA  FOR  H2AOA  DESIGN 

Name 

mm 

Kl 

K4. 

H2AOA 

5.90 

2490 

0.31 

2.52 

TABLE  5.7 

H2AOA  STABILITY  MARGINS 


Name 

Complement. 
Sensitivity 
Vector  Gain 
Margins  (dB) 
Input  of  Plant 

Sensitivity 
Vector  Gain 
Margins  (dB) 
Input  of  Plant 

Phase 
Margins 
(deg) 
Input  of 
Plant 

Complement. 
Sensitivity 
Vector  Gain 
Margins  (dB) 
Output  of 
Plant 

Sensitivity 
Vector  Gain 
Margins  (dB) 
Output  of 
Plant 

Phase 

Margins  (deg) 
Output  of 
Plant 

H2AOA 

[-8.75,4.27] 

[-5.27,  15.57] 

49.36 

[-4.57, 2.98] 

[-3.28,  5.34] 

26.54 

The  high  norm  for  the  tight  tracking  constraint  (  )  indicates  this  design  does  not 


tightly  track  command  inputs.  This  correlates  with  the  5%  tracking  error  at  all  frequencies  that 
was  accepted  in  the  design  setup.  Figures  5.9  and  5.10  show  that  the  time  histories  for  this 
design  do  look  similar  to  the  CLASSIC  (as  intended).  The  stability  margins  for  this  design  are 
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similar  to  the  CLASSIC  design.  This  is  not  surprising  since  the  H2AOA  design  was  supposed  to 
have  the  same  response  (and  therefore  same  handling  qualities)  as  the  CLASSIC. 
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degrees  ,  deg/sec 


Handling  qualities  predictions  for  this  design  were  calculated  using  the  Hoh’s  Bandwidth 
and  Berry  Flight  Path  Overshoot  Criteria.  The  Hoh’s  bandwidth  criteria  predicted  borderline 
level  I/II  with  a  delay  parameter  (tau  p)  of  0.10  and  a  bandwidth  of  2.7  rad/sec.  The  Berry 
Criteria  predicted  level  I  with  a  19%  overshoot.  Consideration  was  given  to  trying  to  tune  this 
design  (by  adjusting  control  use  weights,  ideal  model,  etc.)  to  get  a  better  handling  qualities 

prediction.  However  the  test  team  pilot  comments  from  initial  ground  simulation  indicated  that 

/' 

they  liked  the  response  of  the  design.  Therefore  no  changes  were  made. 

The  dynamics  of  the  H2AOA  control  law  and  verification  of  their  implementation  is 
found  in  Table  D5,  Appendix  D. 

Mixed  H2/H„  Angle  of  Pitch  Rate  Command  Design.  This  design  was  a  tradeoff  of  noise 
rejection/tracking  with  output  stability  margins.  It  was  comprised  of  the  same  H2  problem 
design  formulation  as  for  the  H2AOA  design,  and  the  output  stability  margin  constraint, 

However,  the  ideal  model  for  this  design  was  changed  somewhat  from  the  H2AOA  configuration 
after  some  discussion  with  Calspan  test  pilots.  They  indicated  that  ideal  model  natural 
frequencies  of  as  low  as  2  rad/sec  should  result  in  good  handling  qualities  for  the  landing  task. 
Hence  an  ideal  model  with  a  natural  frequency  of  2  rad/sec  was  chosen  for  this  design.  The 
transfer  function  is  given  by 


Wm  =  —. - 

5  +2s  +  4 


(5.2) 


The  tracking  performance  weight,  Wp,  the  control  use  constraint,  Wcu,  and  the  control 
rate  usage  constraint,  Wcr,  were  reset  to  the  values  in  section  3.3.  This  was  done  because  these 
weights  worked  for  that  pitch  rate  design. 
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The  plant,  actuator,  turbulence,  state  weighting,  and  sensor  noise  models  all  remained 
unchanged  from  previous  designs. 

This  control  design  was  again  given  the  name  MXINI.  H2  and  H„  norms  are  shown  in 
Table  5.8  and  vector  stability  margins  are  shown  in  Table  5.9. 


TABLE  5.8 

NORM  DATA  FOR  MXINI  DESIGN 


Name 

iU 

■an 

mam 

MXINI 

1.63 

261.20 

0.26 

3:15  1 

TABLE  5.9 

MXINI  STABILITY  MARGINS 


Name 

Compliment. 
Sensitivity 
Vector  Gain 
Margins  (dB) 
Input  of  Plant 

Sensitivity 
Vector  Gain 
Margins  (dB) 
Input  of  Plant 

Phase 
Margins 
(deg) 
Input  of 
Plant 

Compliment. 
Sensitivity 
Vector  Gain 
Margins 
(dB) 

Output  of 
Plant 

Sensitivity 
Vector  Gain 
Margins  (dB) 
Output  of 
Plant 

Phase 
Margins 
(deg) 
Output  of 
Plant 

Ii2 

Subproblem 

[-5.86, 28.48] 

[-5.86, 3.47] 

57.53 

[-3.30, 2.39] 

[-2.66, 3.86] 

20.66 

MXINI 

[-12.28,4.96] 

[-5.41, 17.39] 

51.25 

[-8.97, 4.32] 

[-4.87,  12.13] 

44.21 

The  high  norm  for  the  tight  tracking  constraint  ( )  indicates  this  design  does  not 

tightly  track  command  inputs.  However,  recall  that  this  design  has  given  up  some  tracking 
performance  for  better  stability  margins  by  virtue  of  the  mixed  problem  formulation.  A 
comparison  between  the  stability  margins  for  the  H2  portion  of  the  problem  and  the  MXINI  in 
Table  5.9  shows  that  the  output  margins  are  significantly  improved  (as  intended).  Figures  5.9 
and  5.10  show  the  time  histories  for  this  design. 
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ees  ,  deg/sec 


Input  Command  Pulse 


Elevator  Deflection,  5  deg,  4  sec  pulse 


Angle  of  attack,  4  deg,  5  sec  pulse 


Figure  5.9.  Short  Period  Simulations  for  MXINI  -  Noises  off 
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Handling  qualities  predictions  for  this  design  were  calculated  using  the  Hoh’s  Bandwidth  and 
Berry  Flight  Path  Overshoot  Criteria.  The  Hoh’s  bandwidth  criteria  predicted  solid  level  I 
handling  qualities  with  a  delay  parameter  (tau  p)  of  0.025  and  a  bandwidth  of  3 .2  rad/sec. 
However,  the  Berry  Criteria  predicted  level  II  with  a  40%  overshoot. 

This  design  was  produced  in  the  time  period  between  initial  implementation  verification 
in  Buffalo  and  the  final  verification  at  Edwards  AFB  1  day  prior  to  flight  test.  Consequently  the 
implementation  could  not  be  verified  until  the  final  verification  simulation.  It  turned  out  that 
there  was  an  error  in  the  script  file  that  implemented  the  MXINI  control  law.  The  pole-zero 
mismatch  is  shown  in  Table  D6,  Appendix  D.  A  simple  fix  to  the  script  file  within  the 
established  SIMULINK™  control  law  architecture  was  attempted  at  the  end  of  the  ground 
simulation  period  at  Edwards  AFB.  However,  complete  verification  of  the  fix  was  not  made 
prior  to  the  first  HQDT  sortie  due  to  time  constraints. 

5.4  Ground  and  Airborne  Flight  Control  Law  Validation 

Validation  testing  was  accomplished  during  ground  test  by  comparing  the  closed-loop 
time  responses  to  a  PTI  pulse  input  at  the  control  stick.  The  comparison  was  made  between  each 
control  system  coded  on  the  VSS  computer  to  those  generated  in  the  design  process.  Since  the 
design  model  was  a  short  period  approximation  of  the  F-16,  only  the  pitch  rate  and  the  AOA 
responses  were  examined.  The  linear  time  history  response  matches  had  some  discrepancies,  but 
were  considered  an  acceptable  validation  of  the  closed  loop  response.  See  Figures  D1-D5  in 
Appendix  D  for  ground  time  history  validation  plots.  Recall  that  no  time  history  validation  was 
done  for  the  MXINI  configuration  due  to  time  constraints. 
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In  addition  to  comparing  PTI  time  responses  for  validation,  the  test  team  pilots  were  able 
to  “fly”  simulated  approaches  with  each  flight  control  configuration  using  an  integrated  synthetic 
horizon,  Attitude-Direction  Indicator  (ADI),  and  flight  path  marker  display.  This  not  only  gave 
them  a  feel  for  how  each  control  configuration  might  fly,  but  provided  valuable  qualitative 
validation  information  that  was  not  obtained  from  the  time  response  matches.  The  pilots  noticed 
an  apparent  negative  speed  stability  behavior  in  the  aircraft  for  the  H2INI,  KGAIN,  and  H2AOA 
configurations  which  had  not  been  noticed  previously.  Upon  discussion  with  Calspan,  it  was 
explained  that  the  VSS  computer  modified  the  Learjet’s  short  period  dynamics  to  look  like  F-16 
short  period  dynamics,  but  did  not  eliminate  or  modify  the  actual  Learjet  phugoid  that  would  be 
present  inflight.  Consequently,  when  the  Learjet  phugoid  mode  was  included  in  the  simulation, 
which  had  not  been  done  previously,  three  of  the  configurations  appeared  to  have  negative  speed 
stability.  The  Learjet  phugoid  mode  was  not  considered  in  the  design  process  upon  the  advice  of 
the  AFIT  faculty  and  Calspan.  They  felt  this  mode  would  not  have  a  significant  effect  on 
handling  qualities  and  would  add  unnecessary  complexity  and  increased  order  to  the  resulting 
flight  control  designs.  In  actuality,  it  appeared  that  leaving  the  Learjet  phugoid  mode  out  of  the 
design  process  might  have  an  impact  on  the  handling  qualities  of  the  flight  control  designs.  The 
CLASSIC  and  MXAOA  control  configurations  did  not  exhibit  any  negative  speed  stability 
characteristics  during  piloted  ground  simulations.  The  CLASSIC  control  configuration  may 
have  handled  the  additional  dynamics  well  because  it  only  involved  a  small  amount  of  gain  on 
the  AOA  and  pitch  rate  feedback  paths.  Unlike  all  of  the  other  flight  control  configurations, 
there  were  no  controller  dynamics  interacting  with  the  phugoid  mode  in  the  CLASSIC 
configuration.  The  MXAOA  controller  was  dynamic  (i.e.,  a  transfer  function);  however,  it  was 
also  designed  to  handle  model  uncertainty.  Thus,  the  MXAOA  control  configuration  also 
performed  as  expected. 
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Inflight  model  validation  results  closely  matched  the  ground  model  validation  results 
(Figures  D6-D10)  and  were  considered  acceptable.  The  pilots  again  noticed  the  apparent 
negative  speed  stability  characteristic  in  the  H2INI,  H2AIN,  and  H2AOA  configurations.  Some 
of  the  discrepancies  in  the  latter  portions  of  the  time  histories  for  the  H2INI,  H2AIN,  and 
H2AOA  configurations  were  attributed  to  the  apparent  negative  speed  stability.  Again,  the 
MXAOA  and  CLASSIC  configurations  performed  as  expected  with  no  apparent  instabilities. 

The  fix  of  the  MXINI  configuration  was  spot  checked  during  this  flight.  This 
configuration  caused  a  small,  bounded  PIO.  Analysis  of  the  implementation  revealed  the 
intended  fix  had  not  worked.  Furthermore,  a  simple  fix  within  the  existing  VSS  software 
architecture  was  not  possible.  The  decision  was  made  to  drop  MXINI  from  subsequent  flight 
evaluation. 

Post-flight  analysis  using  the  design  controller  and  aircraft  dynamics  revealed  that  one  of 
the  low  frequency  poles  in  the  H2INI,  H2AIN,  and  H2AOA  control  law  configurations  became 
unstable  in  the  closed  loop  dynamics  when  the  Learjet  phugoid  was  included  in  the  aircraft 
dynamics  (See  Tables  D7-D1 1  in  Appendix  D).  This  indicates  that  these  designs  were  not  robust 
to  the  addition  of  unmodeled,  low  frequency  dynamics  such  as  the  phugoid. 

5.5  Summary 

Taking  care  of  the  logistics  and  setting  the  methodology/procedures  to  be  used  in  this 
flight  test  program  was  they  key  to  the  successful  verification  and  implementation  of  the  HAVE 
INFINITY  II  control  laws.  It  was  this  process  that  produced  the  decision  to  make  a  trip  Calspan 
facilities  two  months  prior  to  flight  test  for  the  purpose  of  initial  verification  and  validation  of 
the  control  laws.  The  early  discovery  of  problems  with  the  implementation  of  the  H2INI/H2AIN 
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and  MXINI/MXAIN  control  laws  allowed  time  for  redesign  and  supplemental  implementation  of 
the  new  designs  (CLASSIC,  H2AOA  and  (new)  MXINI). 

Ultimately,  the  verification  and  validation  phase  of  the  test  program  confirmed  the  flight 
control  laws  implemented,  and  flown,  were  the  same  as  those  designed,  except  for  the  MXINI 
configuration.  The  dynamics  of  the  actual  aircraft  were  different  than  the  dynamics  of  the  design 
model  due  to  the  exclusion  of  the  phugoid  mode  from  the  design  model.  As  a  result,  the  H2INI, 
H2AIN,  and  H2AOA  control  configurations  were  negatively  impacted  due  to  a  lack  of 
robustness  to  low  frequency  differences  between  the  design  model  and  the  actual  aircraft.  The 
system  should  have  been  tested  with  the  phugoid  mode  early  in  the  design  process.  Then,  if  a 
stability  problem  was  found,  a  redesign  could  have  been  accomplished  by  adding  additional 
uncertainty  to  the  short  period  model  or  simply  including  the  phugoid  mode  in  the  design.  The 
bottom  line  is  that  early  PCS  ground  tests  should  be  conducted  with  the  highest  fidelity  aircraft 
model  available  to  allow  time  for  required  redesigns  prior  to  flight  test. 

Table  5.10  provides  a  summary  of  the  designs  that  were  successfully  implemented, 
verified,  and  validated  for  flight  test. 


TABLE  5.10 

SUMMARY  OF  VERIFIED  AND  VALIDATED  FLIGHT  CONTROL  DESIGNS  SELECTED 
FOR  HANDLING  QUALITIES  FLIGHT  TEST 


Method  of 
Design 

Type  of 
Design 

Hoh’s  Bandwidth 
HQ  Prediction 

Berry  Flight  Path 
Overshoot  Criteria  HQ 
Prediction 

Feedback 

Gain 

AOA  &  Pitch 
rate  feedback 

Level  1 

Level  1 

H2 

AOA  Cmd 

Level  1 

Level  1 

H2INI 

H2 

Pitch  Rate 
Cmd 

Level  1 

Level  1 

H2AIN 

H2 

AOA  Cmd 

Level  1 

Level  1 

MXAOA 

H2/H„ 

AOA  Cmd 

Level  1 

Level  1 
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VI.  Flight  Test  Results 


This  chapter  reviews  the  handling  qualities  results  of  the  five  configurations  flight  tested 
based  on  the  test  team’s  analysis  of  the  quantitative  and  qualitative  data.  All  of  the  raw 
quantitative  and  qualitative  handling  qualities  data  is  presented  in  Appendix  E.  The  final  section 
is  an  overall  summary  comparing  handling  qualities  flight  test  data  to  their  respective 
predictions.  Table  6.1  summarizes  the  test  team  pilots’  experience. 


TABLE  6.  1 


EVALUATION  PILOT  EXPERIENCE 


Pilot 

Primary  A/C 

Flight  Hours 

Boe  (USAF) 

F-15C 

1,500 

Stevenson  (USN) 

F-14A/B 

1,300 

Cantiello  (ItAF) 

F-3  TORNADO 

1,100 

6. 1  CLASSIC  Flight  Control  Configuration 

The  evaluation  pilots  performed  four  straight-in  and  four  horizontal  offset  landings  with 
the  CLASSIC  flight  control  configuration.  See  Figures  El  through  E4  (Appendix  E)  for  Cooper- 
Harper  (C-H)  and  PIO  histograms.  Level  I C-H  ratings  were  given  for  each  landing  with  one 
exception  that  received  C-H  4.  All  Level  I  ratings  were  associated  with  PIO  1,  and  the  C-H  4 
received  a  PIO  3.  All  the  C-H  and  PIO  ratings  were  unaffected  by  differences  in  wind  conditions 
(Appendix  E). 

All  pilots  agreed  that  the  CLASSIC  flight  control  configuration  provided  a  predictable, 
linear  initial  response,  although  the  response  was  somewhat  slower  than  several  of  the  other 
configurations.  The  majority  of  Level  I  ratings  confirmed  the  satisfactory  handling  qualities  of 
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this  configuration;  however,  while  performing  HQDT,  pilot  two  noted  that  the  aircraft  response 
was  not  linear  with  respect  to  the  size  of  the  input.  If  the  aircraft  response  was  truly  linear,  a 
larger  input  should  result  in  a  larger  output  over  the  same  period  of  time.  In  other  words,  a  faster 
initial  response  should  occur  with  a  larger  input.  Pilot  two  found  that  large  stick  inputs  did  not 
produce  a  faster  pitch  rate  response  than  smaller  inputs.  While  pilot  two  saw  this  problem 
consistently  during  the  HQDT  tasks,  he  did  not  notice  it  during  any  of  the  landing  tasks.  On  one 
landing,  however,  pilot  three  downgraded  the  configuration  to  C-H  4  because  he  had  to  sacrifice 
desired  performance  when  he  encountered  an  unexpected  pitch  response  in  the  flare.  Postflight 
analysis  showed  that  pilot  three  increased  the  size  and  rate  of  his  stick  inputs  at  the  initiation  of 
the  flare,  but  the  aircraft  did  not  respond  faster  to  these  inputs.  This  forced  pilot  three  to  lower 
his  gains  and  accept  a  long  landing.  Both  pilot  two  and  pilot  three’s  comments  about  the  pitch 
response  were  consistent  with  a  rate  limitation  problem  in  the  flight  control  system;  however, 
flight  data  evaluation  did  not  support  a  rate  limitation  problem.  No  other  explanation  for  the 
nonlinear  pitch  response  was  found.  Inputs  large  enough  to  uncover  the  nonlinear  pitch  response 
were  not  required  on  a  majority  of  the  landings  so  the  problem  had  little  impact  on  the  C-H  and 
PIO  ratings.  Overall,  the  CLASSIC  flight  control  configuration  was  rated  Level  1. 

6.2  H2AOA  Flight  Control  Configuration 

The  evaluation  pilots  performed  four  straight- in  and  four  horizontal  offset  landings  with 
the  H2AOA  flight  control  configuration.  See  Figures  E5  through  E8  (Appendix  E)  for  C-H  and 
PIO  histograms.  The  H2AOA  flight  control  configuration  received  the  worst  ratings  of  all  the 
configurations  flown  in  the  landing  tasks.  Level  II  C-H  ratings  were  given  for  all  landings  with 
one  exception  that  received  C-H  3.  Two  landings  were  given  C-H  6  ratings  with  PIO  4.  Another 
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landing  was  given  a  C-H  6  rating  with  a  PIO  3.  In  general,  the  ratings  for  the  horizontal  offset 
landing  task  were  better  than  the  ratings  for  the  straight- in  landing  task,  although  one  of  the  C-H 
6  ratings  was  given  for  a  horizontal  offset  task.  The  ratings  of  this  configuration  also  worsened 
with  higher  winds.  The  only  Level  I  rating  (C-H  3)  was  given  during  very  light  winds 
(Appendix  E). 

All  three  pilots  agreed  the  long-term  pitch  axis  response  was  slowly  divergent.  A 
deviation  in  airspeed  started  a  pitching  moment,  resulting  in  slowly  increasing  stick  forces  that 
continued  to  build  in  the  initial  direction  until  the  pilot  retrimmed  the  aircraft.  The  pitching 
moment  was  opposite  of  what  was  expected.  When  the  airspeed  slowed  from  a  trimmed 
condition,  forward  stick  pressure  was  required  to  keep  the  nose  from  rising.  The  amount  of 
disturbance  required  to  start  a  divergence  was  very  small.  A  1-knot  deviation  in  airspeed  was 
enough  to  initiate  the  slow  buildup  of  stick  force.  The  small  size  of  the  disturbance,  required  to 
initiate  the  divergence,  made  the  direction  unpredictable  to  the  pilot;  such  that  the  pilot  had  no 
indication  whether  the  divergence  would  result  in  a  slow  buildup  of  push  or  pull  force.  One  pilot 
commented  that  flying  the  configuration  “felt  like  trying  to  balance  on  a  bowling  ball.”  This 
problem  increased  pilot  workload  because  the  pilots  had  to  pay  constant  attention  to  airspeed  and 
pitch  attitude,  and  continually  retrim  the  aircraft.  These  characteristics  became  less  apparent 
during  the  offset  landing  tasks  when  the  pilots  were  exercising  tighter  control  of  the  aircraft. 
Under  tight  control,  the  aircraft  was  not  allowed  to  diverge  from  trim,  and  the  C-H  and  PIO 
ratings  improved.  Regardless  of  task,  all  pilots  felt  the  nose  pitch  down  after  main  wheel 
touchdown.  The  pilots  were  unable  to  keep  the  nose  from  dropping  with  a  reasonable  amount  of 
stick  force  and  did  not  attempt  to  use  large  displacements  of  aft  stick  to  arrest  the  nose 
movement.  Pilot  two  felt  the  stick  forces  increase  dramatically  in  the  flare  on  both  a  straight-in 
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task  and  a  horizontal  offset  task,  which  made  the  pitch  attitude  very  hard  to  control.  This 
characteristic  generated  the  two  PIO  4  ratings  given  by  pilot  two.  The  small  amplitude  pitch 
oscillations  caused  pilot  two  to  freeze  the  stick  and  he  achieved  desired  performance  only 
because  the  PIO  occurred  very  late  in  the  flare.  Overall,  the  long-term  divergent  response  of  this 
flight  control  system  had  a  large  impact  on  the  majority  of  the  landings.  The  H2AOA  flight 
control  configuration  was  rated  Level  II. 

The  divergent  nature  of  the  H2AOA  flight  control  configuration  was  not  completely 
unexpected.  As  discussed  in  Chapter  V  (section  5.4),  ground  simulations  completed  Just  prior  to 
flight  test  revealed  that  the  system  had  a  slowly  divergent  first  order  mode.  This  problem  was 
not  predicted  with  the  short  period  approximation  of  the  aircraft  model  used  in  design  and  during 
initial  ground  tests.  The  problem  was  only  found  when  the  phugoid  mode  was  included  during 
the  final  ground  simulation  1  day  prior  to  the  start  of  flight  test. 

6.3  H2INI  Flight  Control  Configuration 

The  evaluation  pilots  performed  eight  straight-in  and  eight  horizontal  offset  landings 
with  the  H2INI  flight  control  configuration.  See  Figures  E9  through  E12  (Appendix  E)  for  C-H 
and  PIO  histograms.  Pilot  one  consistently  rated  this  configuration  Level  II  while  pilot  three 
consistently  rated  this  configuration  Level  1.  Pilot  two  rated  the  configuration  Level  I  for  the 
straight-in  task,  but  gave  both  Level  I  and  Level  II  ratings  for  the  horizontal  offset  task.  The  PIO 
ratings  by  all  the  pilots  fell  between  PIO  1  and  PIO  3  with  the  majority  of  the  ratings  PIO  1  or 
PIO  2.  The  C-H  and  PIO  ratings  were  worse  for  the  horizontal  offset  task  than  they  were  for  the 
straight-in  task.  This  configuration  also  received  poorer  ratings  under  higher  winds  (Appendix 
E). 
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All  the  pilots  liked  the  short-term  pitch  response,  describing  it  as  smooth  and 
predictable.  All  three  pilots  could  change  the  pitch  attitude  rapidly;  however,  pilot  one 
commented  on  one  of  his  four  landings  that  the  response  felt  sluggish.  No  explanation  was 
found  for  this  comment.  All  three  pilots  flew  this  configuration  with  winds  from  8  knots  to  20 
knots.  During  these  flights,  each  pilot  noticed  an  airspeed  sensitivity  and  divergence  similar  to 
the  H2AOA  configuration,  but  not  as  objectionable.  Under  these  windy  conditions,  pilot  three 
also  noticed  “negative  speed  stability”  during  the  go-around  portion  of  a  low  approach.  Pilot 
Two’s  second  look  at  the  configuration  occurred  on  a  day  when  the  winds  varied  from  calm  to  10 
knots.  During  this  sortie,  pilot  two  only  noticed  the  divergent  characteristics  when  the  wind  was 
above  5  knots.  Pilot  two  gave  the  configuration  Level  II  ratings  during  the  windy  conditions  due 
to  the  increased  workload  required  to  maintain  airspeed.  Pilot  one  was  able  to  consistently 
achieve  desired  performance  with  this  configuration  but  he  found  the  divergent  nature  of  the 
system  to  be  objectionable  enough  to  warrant  Level  II  ratings.  All  the  pilots  noticed  a  slight 
tendency  for  the  nose  to  pitch  down  during  the  landing  roll.  This  pitching  moment  was  more 
benign  than  the  H2AOA  configuration.  Like  the  H2AOA  design,  the  slight  divergent  nature  of 
the  aircraft  with  this  flight  control  configuration  was  predicted  during  the  final  ground  test. 
Overall,  the  slight  divergent  nature  of  the  pitch  response  had  an  impact  on  the  landing  tasks. 

The  H2INI  configuration  was  rated  between  Level  I  and  Level  II. 

6.4  H2AIN  Flight  Control  Configuration 

The  evaluation  pilots  performed  six  straight-in  and  six  horizontal  offset  landings  with 
the  H2AIN  flight  control  configuration.  See  Figures  El 3  through  E16  (Appendix  E)  for  C-H  and 
PIO  histograms.  The  C-H  ratings  for  this  configuration  ranged  from  C-H  3  to  C-H  5  with  the 
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majority  of  the  ratings  being  C-H  4.  The  PIO  ratings  for  this  configuration  ranged  from  PIO  1  to 
PIO  3  with  the  majority  of  ratings  being  PIO  2.  Although  each  pilot’s  C-H  and  PIO  ratings 
varied,  the  overall  distribution  of  the  ratings  were  virtually  the  same  for  both  the  straight-in  and 
horizontal  offset  tasks.  None  of  the  ratings  were  affected  by  the  winds  (Appendix  E). 

The  H2AIN  flight  control  configuration  exhibited  many  of  the  same  characteristics  as 
the  H2INI  flight  control  configuration,  although  to  a  greater  degree.  This  similarity  was  logical 
as  the  only  difference  between  H2ArN  and  H2INI  was  a  filter  placed  in  the  command  channel  of 
the  H2ArN  configuration.  This  filter  was  designed  to  make  the  original  pitch  rate  tracking 
system  behave  more  like  an  AOA  tracking  system.  While  this  objective  was  achieved,  the  filter 
made  the  initial  pitch  response  more  sluggish.  The  pilots  did  not  like  this  sluggishness  and 
commented  that  the  pitch  response  “takes  a  while  to  get  going.”  The  sluggishness  was  very 
noticeable  when  the  stick  was  displaced  slightly  and  held  while  the  pitch  response  was  observed. 
The  pitch  rate  accelerated  slowly  until  it  stabilized  at  a  steady  state  value.  The  added  filter  in  the 
H2AIN  design  may  have  caused  the  divergent  nature  of  the  system  to  be  more  noticeable.  All 
the  pilots  felt  that  the  long-term  pitch  axis  response  diverged  in  a  manner  similar  to  the  H2AOA 
configuration,  although  at  a  slower  rate.  This  divergence  made  the  configuration  very  sensitive 
to  airspeed  deviations  from  trim,  therefore  pilot  workload  increased  trying  to  maintain  airspeed. 
Pilot  two  commented  that  he  was  5  KIAS  below  his  trim  airspeed  during  the  approach  and  had  to 
push  on  the  stick;  however,  during  the  flare,  all  the  pilots  commented  that  the  nose  felt  heavy. 
This  made  the  aircraft  hard  to  control  precisely  in  the  flare  and  forced  the  pilots  to  adopt  a  low- 
gain  technique  to  obtain  desired  performance.  The  pilots  were  unable  to  keep  the  nose  from 
dropping  after  landing  due  to  the  heavy  stick  forces.  This  nose  down  motion  was  worse  than 
H2INI,  but  not  as  bad  as  H2AOA.  Like  the  H2AOA  and  H2INI  configurations,  the  H2AIN 
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system  was  found  to  be  divergent  during  the  final  ground  test.  Both  the  sluggishness  of  the  pitch 
response  and  the  divergent  nature  of  the  system  had  an  impact  on  the  landing  tasks.  The  H2AIN 
configuration  was  rated  between  Level  I  and  Level  II. 

6.5  MXAOA  Flight  Control  Configuration 

The  evaluation  pilots  performed  1 1  straight-in  and  10  horizontal  offset  landings  with  the 
MXAOA  flight  control  configuration.  See  Figures  El 7  through  E20  (Appendix  E)  for  C-H  and 
PIO  histograms.  Pilots  one  and  two  consistently  rated  this  configuration  Level  I  while  pilot  three 
consistently  rated  this  configuration  Level  II.  PIO  ratings  ranged  from  PIO  1  to  PIO  3.  The 
winds  did  not  have  an  effect  on  the  C-H  and  PIO  ratings  (Appendix  E). 

The  MXAOA  flight  control  configuration  was  stable  with  no  tendency  for  a  long-term 
pitch  divergence.  Stability  was  predicted  for  this  configuration  during  the  final  ground  test. 
While  the  long-term  pitch  response  was  acceptable,  the  pilots  felt  that  the  short-term  pitch 
response  was  fast  and  too  lightly  damped.  Each  of  the  pilots  noticed  a  small  overshoot  in  the 
aircraft’s  pitch  response  which  they  described  as  a  “pitch  bobble.”  This  bobble  could  be 
eliminated  by  anticipating  the  aircraft’s  response  and  adjusting  the  size  of  the  input  accordingly. 
Pilots  one  and  two  were  able  to  eliminated  this  pitch  overshoot  with  what  they  considered  to  be 
minimal  compensation  and  gave  the  configuration  Level  I  C-H  ratings.  Pilot  three,  however, 
found  the  pitch  bobble  objectionable  during  the  flare  and  had  to  sacrifice  desired  performance  on 
several  landings.  During  these  landings,  the  pilot  cautiously  delayed  the  power  reduction  until 
the  pitch  bobble  damped  out.  This  caused  Pilot  Three  to  carry  excess  power  into  the  flare  and 
float  beyond  the  desired  touchdown  point.  On  the  landings,  when  pilot  three  was  able  to  achieve 
desired  performance;  he  still  felt  the  workload  was  high  enough  to  warrant  Level  II  ratings.  The 
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pitch  sensitivity  was  also  noted  during  HQDT.  Both  pilot  two  and  three  gave  the  configuration 
PIO  4  ratings  during  HQDT  tasks.  The  HQDT  results  indicate  that  this  configuration  was  PIO 
prone  to  large  amplitude,  high  frequency  inputs.  This  tendency  may  have  been  the  cause  of  the 
unwanted  pitch  motions  seen  by  pilot  three.  Overall,  the  pitch  sensitivity  of  this  flight  control 
system  had  an  impact  on  some  of  the  landings  and  the  MXAOA  flight  control  configuration  was 
rated  between  Level  I  and  Level  II. 

6.6  Summary 

Overall  the  handling  qualities  results  were  excellent.  Most  of  the  Level  II  comments 
related  to  the  instability  in  the  H2INI,  H2AIN,  and  H2AOA  designs  were  attributed  to  the 
phugoid  mode.  The  H2INI,  H2AIN,  and  H2AOA  designs  produced  acceptable,  but 
unsatisfactory,  handling  qualities,  and  may  have  been  closer  to  satisfactory  had  the  phugoid 
mode  been  included  in  the  design  aircraft  model.  The  results  of  this  program  simply  reinforce 
that  early  verification  and  validation  of  design  models  through  the  use  of  the  highest  fidelity 
simulation  possible,  must  be  supported  and  required  by  the  Responsible  Test  Organization 
(RTO). 

It  should  also  be  pointed  out  here  that  the  variation  in  the  number  of  landings  for  each 
flight  control  configuration  was  intentional.  The  requirement  for  success  was  that  each  pilot 
evaluate  the  flight  control  configuration  twice  in  the  landing  phase  (1  straight-in,  1  offset),  unless 
the  configuration  received  a  C-H  8,  9,  or  10  (in  which  case  it  would  be  dropped).  Beyond  that, 
the  configurations  evaluated  were  determined  by  the  TC.  The  borderline  Level  I/II 
configurations  got  more  landings  in  an  attempt  to  eliminate  this  ambiguity.  While  not  successful 
in  that  respect,  the  added  data  was  helpful  in  confirming  previous  pilot  comments  and  ratings. 
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Table  6.2  compares  the  flight  test  handling  qualities  results  with  the  designs  predictions. 


TABLE  6.2 

HANDLING  QUALITIES  RESULTS  VS.  PREDICTIONS 


Name 

Type  of 
Design 

Hoh’s  Bandwidth 
HQ  Prediction 

Berry  Flight  Path 
Overshoot  Criteria  HQ 
Prediction 

Flight  Test  Results 

CLASSIC 

AOA  &  Pitch 
rate  feedback 
gain 

Level  1 

Level  1 

Level  1 

H2AOA 

H2  AOA  Cmd 

Level  1 

Level  1 

Level  2 

H2INI 

H2  Pitch  Rate 
Cmd 

Level  1 

Level  1 

Level  1/2 

H2  AOA  Cmd 

Level  1 

Level  1 

Level  1/2 

MXAOA 

H2/H„  AO  A 
Cmd 

Level  1 

Level  1/2 

Since  the  actual  aircraft  included  a  mode  not  present  in  the  design  model,  the 
comparison  in  Table  6.2  is  “apples  to  oranges”.  To  remedy  this,  the  handling  qualities 
predictions  were  recomputed  with  the  phugoid  mode  included  in  the  aircraft  design  model.  The 
plots  of  the  Hoh’s  Bandwidth  Criteria  and  time  history  predictions  are  in  Appendix  F.  A  plot  of 
the  Hoh’s  Bandwidth  predictions,  for  each  configuration  without  the  phugoid  is  included  for  ease 
of  reference.  These  results  are  summarized  in  Table  6.3. 
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TABLE  6.3 


HANDLING  QUALITIES  RESULTS  VS.  PREDICTIONS  WITH  PHUGOID  MODE 


Name 

Type  of 
Design 

Hoh’s  Bandwidth 
HQ  Prediction 
(w/phugoid) 

Berry  Flight  Path 
Overshoot  Criteria  HQ 
Prediction  (w/phugoid) 

Flight  Test  Results 

AOA  &  Pitch 
rate  feedback 
gain 

Level  1 

Level  1 

Level  1 

/f^AOACmd 

Level  2 

Divergent  1 

Level  2 

H2INI 

H2  Pitch  Rate 
Cmd 

Level  2 

Divergent^ 

Level  1/2 

H2AIN 

^2AOACmd 

Level  2 

Divergent  f 

Level  1/2 

MXAOA 

H2/H^  AOA 
Cmd 

Level  1 

Level  1 

Level  1/2 

1.  Handling  qualities  prediction  were  not  possible  from  the  Berry  criteria  for  the  divergent  control  laws. 


The  data  in  Table  6.3  show  good  correlation  between  the  predictions  and  actual  results. 
Moreover,  these  results  would  tend  to  nullify  the  hypothesis  made  earlier  that  the  Hoh’s 
bandwidth  criteria  would  not  do  a  good  job  of  predicting  handling  qualities.  While  it  cannot  be 
concluded  from  these  results  that  the  Hoh’s  criteria  will  work  will  give  good  predictions  in  all 
cases  (history  has  proven  that  [23]),  it  did  do  a  good  job  here.  Of  course  the  borderline  Level 
I/II  ratings  make  it  much  easier  to  correlate  the  flight  test  data  with  predictions! 

The  borderline  ratings  are  important  in  that  they  indicate  there  are  a  lot  of  good 
characteristics  about  the  designs.  They  also  indicate  there  are  some  bad  characteristics,  but  not 
so  bad  such  that  the  design  should  be  scrapped.  The  problems  caused  by  the  negative  speed 
stability  could  be  fixed  by  simply  re-running  the  H2  designs  with  the  phugoid  mode  included  in 
the  aircraft  model.  The  sluggishness  of  the  H2AIN  configuration  was  likely  caused  by  the  filter 
dynamics  (since  it  was  the  only  difference  from  H2INI).  This  problem  could  be  easily  remedied 
by  taking  the  filter  out  which  would  leave  H2INI!  The  split  in  MXAOA  came  down  to  pilot 
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preference.  However,  the  fact  that  a  slight  pitch  bobble  was  noticed  by  all  pilots  increasing  the 
ideal  model  damping  would  probably  eliminate  the  problem. 
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VII.  Conclusions  and  Recommendations 


7.1  Conclusions 

The  overall  objective  of  this  research  was  to  perform  an  inflight  handling  qualities 
evaluation  of  several  multiobjective,  optimal  flight  control  designs.  This  objective  was  met. 

Overall  the  results  were  good  and  indicate  that  the  multiobjective,  optimal  model 
following  design  techniques  used  in  this  research  can  produce  at  least  Level  I/II  (acceptable) 
handling  qualities. 

Specific  objectives  were  categorized  into  two  phases:  development  and  simulation 
phase  (phase  I),  and  the  flight  test  phase  (phase  II). 

Phase  I  specific  objectives  included  the  development  and  simulation  of  H2 ,  H„,  ju, 
mixed  H2lH„ ,  and  mixed  H2/n  problems.  All  of  these  objectives  were  met. 

Selection  of  the  ideal  model  used  for  model  following  in  the  designs  was  based  on 
handling  qualities  predictions  using  Hoh’s  Bandwidth  Criteria  and  academic  knowledge  of  what 
a  good  closed  loop  model  would  be.  Initially,  it  was  thought  that  the  ideal  model  would  be  the 
sole  determinant  of  handling  qualities  for  all  of  the  designs  produced.  However  within  the 
context  of  a  specific  ideal  model,  it  was  found  that  the  selection  of  the  control  rate  usage  and 
tracking  performance  weights  could  significantly  influence  the  Hoh’s  handling  qualities 
predictions  of  a  given  H2  design  by  changing  the  bandwidth  of  the  closed  loop  system. 

Increases  to  control  rate  limiting  weight  degraded  handling  qualities  predictions  by 
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producing  systems  with  smaller  bandwidth  that  rolled  off  rapidly  in  phase  at  frequencies  above 
the  bandwidth  frequency.  This  became  important  because  of  the  selection  of  the  tracking 
weights.  The  higher  the  DC  gain  of  the  tracking  weight,  the  better  the  tracking  performance. 
However,  demanding  tight  tracking  many  times  caused  rate  limiting.  The  solution  to  this  was  to 
increase  the  rate  limit  weight.  The  bottom  line  was  that  there  was  a  tradeoff  between  these  two 
weights  in  the  designs  considered  in  this  study.  This  resulted  in  a  fair  amount  of  iteration  during 
the  design  process  to  get  level  I  handling  qualities  predictions.  While  this  finding  may  be 
specific  to  the  work  done  here,  it  serves  as  a  place  to  start  when  the  designer  is  looking  to 
manipulate  the  handling  qualities  of  a  model  following  design. 

The  final  objective  of  phase  I  was  to  pick  the  best  design  candidates  for  flight  test. 
Designs  were  selected  for  flight  test  based  on  the  following  criteria  and  constraints: 

1 .  The  compensator  was  stable  (All  poles  in  LHP). 

2.  Design  handling  qualities  must  have  been  predicted  Level  I  using  the  Hoh’s 
Bandwidth  Criteria  (Landing  Phase).  The  Berry  Flight  Path  Overshoot 
Criteria  was  also  considered,  but  not  used  as  a  primary  criteria  for  flight  test 
selection. 

3.  Designs  were  desired  to  have  6  dB  of  gain  margin  and  45°  of  phase  margin. 

4.  Designs  should  have  good  command  tracking  characteristics. 

All  of  the  desired  design  types  were  completed  and  evaluated  for  flight  test  in  phase  I. 
The  designs  selected  for  flight  test  had  the  best  overall  combination  of  characteristics  for  the 
criteria  listed.  There  weren’t  any  //  designs  or  single  norm  designs  represented  in  the  flight 
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test  candidates.  The  jii  designs  evaluated  were  disqualified  as  not  implementable  because  of  RHP 
compensator  poles.  All  of  the  formulations  were  used  solely  as  singular  constraints  for  the 
mixed  H2lHa  problems  and  never  intended  to  be  stand  alone  designs.  This  was  done  to  allow  the 
designer  to  better  manage  the  tradeoff  process  thereby  facilitating  the  development  of  designs 
having  the  desired  properties. 

Limitations  associated  with  the  Calspan  Learjet  led  to  the  use  of  the  Euler  integration 
routine  in  MATLAB™  SIMULINK'*’’^  for  inflight  simulation  of  the  control  laws.  This  in  turn 
led  to  changes  in  the  H2INI/H2AIN  configuration  and  disqualified  the  MXTNI/MXAfN 
configuration.  Three  new  designs  were  generated  (CLASSIC,  H2AOA,  and  MXINI)  to  replace 
the  two  that  were  disqualified.  The  optimal  design  ideal  models  were  based  on  dynamics  that 
had  proven  good  handling  qualities  from  previous  Calspan  flight  test. 

Phase  II  objectives  were  to  flight  test  the  control  designs  selected  for  flight  test  from 
phase  I,  collect  Cooper-Harper  handling  qualities  ratings,  compare  those  ratings  with  the  design 
phase  predictions,  and  draw  conclusions  with  regards  to  any  apparent  trends  in  predictions  vs. 
actual  ratings. 

Proper  model  verification,  validation  on  the  ground  and  inflight  was  a  critical  step  that 
allowed  for  a  satisfactory  evaluation  of  the  handling  qualities  of  the  HAVE  INFINITY  II  flight 
control  system  designs.  The  Calspan  Learjet  did  an  excellent  job  of  simulation  on  the  ground 
and  in  the  air.  This  was  confirmed  by  the  correlation  of  the  ground  and  flight  test  time  response 
matches  and  pilot  qualitative  comments.  Handling  qualities  results  indicate  that  the  optimal 
design  methods  used  gave  Level  II  or  borderline  Level  I/II  handling  qualities.  Most  of  the  Level 
II  comments  related  to  the  instability  in  the  H2INI,  H2AIN,  and  H2AOA  designs  caused  by  the 
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phugoid  mode  that  was  unaccounted  for  in  the  design  process.  These  designs  may  have  been 
rated  closer  to  Level  I  if  the  phugoid  had  been  included  in  the  design  aircraft  model. 

Design  handling  qualities  predictions  of  the  Hoh  criteria  and  flight  test  handling  qualities 
ratings  correlated  very  well  after  the  phugoid  mode  was  included  in  the  aircraft  model  (Tables 

6.2  and  6.3). 

The  negative  speed  stability  and  inability  to  properly  implement  MXINI  obscured  any 
possible  comparisons  of  pilot  preferences  between  AOA  and  pitch  rate  command  systems. 

7.2  Recommendations 

•  Further  flight  test  of  multiobjective,  optimal  control  designs  (like  those  produced  in 
this  research)  should  emphasize  testing  some  the  purported  benefits  of  these  design 
techniques  such  as  turbulence  rejection  and  design  robustness  once  Level  I  handling 
qualities  are  established.  This  will  be  the  true  measure  of  whether  or  not  these 
methods  offer  any  advantages  over  more  simple  conventional  methods  like  feedback 
gains  and  root  locus  techniques.  Evaluating  turbulence  rejection  could  be  done  by 
performing  the  same  task  in  calm  and  turbulent  conditions  and  comparing  the 
handling  qualities.  Evaluating  design  robustness  would  be  done  by  systematically 
varying  the  aircraft  dynamics  off  nominal  (for  a  given  task  in  the  same  flight 
conditions)  and  documenting  at  what  level  of  deviation  the  handling  qualities 
degraded. 

•  Involve  the  organization  that  will  be  implementing  and  flying  the  control  designs 
(in  this  case  Calspan)  in  the  concept  phase  of  the  research.  These  organizations 
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usually  have  volumes  of  practical  experience  that  can  save  research  time  and 
positively  influence  the  outcome  of  the  results.  Additionally,  the  contractor  can  get 
a  clear  appreciation  of  what  the  project’s  goals  and  methods  are.  They  can  then 
clarify  implementation  limitations  so  that  research  time  will  not  be  wasted 
producing  something  that  can’t  be  implemented  (as  was  done  with  the 
SIMULINK™  integration  routines).  The  bottom  line  is  that  all  significant  players 
in  the  research  program  should  be  involved  from  the  beginning. 

•  Early  PCS  ground  tests  should  be  conducted  with  the  highest  fidelity  aircraft  model 
available  to  allow  time  for  required  redesigns  prior  to  flight  test.  This  lesson  was 
learned  by  both  the  HAVE  INFINITY  in  the  verification  phase  and  the  HAVE 
INFINITY  II  test  team  in  the  validation  phase. 

•  The  model  following  technique  should  be  used  when  designing  an  aircraft  flight 
control  law  that  will  be  subject  to  handling  qualities  evaluations.  The  experience  in 
this  research  shows  that  the  optimal  multiobjective  methods  do  a  good  job  of  model 
following  when  a  reasonable  tracking  constraint  is  included.  By  taking  advantage  of 
this  fact,  and  following  an  ideal  model  with  known  good  handling  qualities,  the 
handling  qualities  “unknown”  in  the  flight  control  design  can  be  eliminated.  This 
leaves  the  designer  the  flexibility  to  use  weights  and  constraints  to  build  in  other 
capabilities  in  the  control  law  (turbulence  rejection,  robustness  to  uncertainty,  better 
stability  margins,  etc.). 
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Appendix  A 


Cooper-Harper  and  PIO  Rating  Scales 


A.  1  Cooper-Harper  Rating 


A  Cooper-Harper  (C-H)  rating  was  given  for  each  landing  task.  Figure  B4  was  used  by 
the  test  conductor  to  aid  the  pilot  in  determining  the  appropriate  C-H  rating. 


lEtPACY  HORSHKnSDT&k 


AIRCRAFT 


PEWiNDS  ON  THE  PE,OT  IN  SELECTED  PJ 
TaSORREOtmEDOPEBATTON  RA' 


lExcellent 

Highly  desireable 

Klot  compensation  not  a  factor  for 
desired  performance 

Good 

Pilot  compensation  not  a  factor  for 

Neglibable  deficeincies 

desired  performance 

Fair-  Some  mildly 
unpleasant  deficiencies 

Minimal  pilot  compensation  required 
for  desired  performance 

Deficiencies 

k 

Minor  but  anoying 
deficiencies 

Desired  performance  requires 
moderate  pilot  compensation 

warrant 

1  Moderately  objectionable  Adequate  performance  requires  | 

improvement 

deficencies 

considerable  pilot  compensation 

III 

Very  objectionable  but 

Adequate  performance  requires 

tolerable  deficiencies 

extensive  pilot  compensation 

V  ■■  . 

Major 

Adequate  performance  not  attainable  wil 

maximum  tolerable  pilot  compensation. 

1  Deficiencies 

Controllability  not  in  question 

require 

Major 

Considerable  pilot  compensation  is 

^  improvement' 

r 

deficencies 

required  for  control 

Major 

Intense  pilot  compensation  is  required 

deficencies 

to  retain  control 

bb 

' 

Improvement! 

Major 

Control  will  be  lost  during  some 

^  mandatory  1 

h 

deficencies 

portion  of  required  operation 

Figure  Al.  Cooper-Harper  Rating  Scale 


A.2  Pilot-in-the-loop  Oscillation  Rating 

A  pilot-in-the-loop  oscillation  PIO  rating  was  given  for  handling  qualities  during 
tracking  (HQDT)  with  an  C-23target  aircraft,  and  for  each  landing  task.  Figure  A2  was  used  by 
the  test  conductor  to  aid  the  pilot  in  determining  the  appropriate  PIO  rating.  Descriptions  for  the 


A-2 


PIO  ratings  are  shown  in  Figure  A3. 


Figure  A2.  Pilot-in-the-Loop  Oscillation  Rating  Decision  Tree 
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Pilot-in-the-loop  Oscillation  (PIO)  Scale 

Description 

Numerical 

Rating 

No  tendency  for  pilot  to  induce  undesirable  motions. 

1 

Undesirable  motions  tend  to  occur  when  pilot  initiates  abrupt 
maneuvers  or  attempts  tight  control.  These  motions  can  be 
prevented  or  eliminated  by  pilot  techniques. 

2 

Undesirable  motions  easily  induced  when  pilot  initiates  abrupt 
maneuvers  or  attempts  tight  control.  These  motions  can  be 
prevented  or  eliminated,  but  only  at  sacrifice  to  task  performance  or 
through  considerable  pilot  attention  and  effort. 

3 

Oscillations  tend  to  develop  when  pilot  initiates  abrupt  maneuvers 
or  attempts  tight  control.  Pilot  must  reduce  gain  or  abandon  task  to 
recover. 

4 

Divergent  oscillations  tend  to  develop  when  pilot  initiates  abrupt 
maneuvers  or  attempts  tight  control.  Pilot  must  open  loop  by 
releasing  or  freezing  the  stick. 

5 

Disturbance  or  normal  pilot  control  may  cause  divergent 
oscillation.  Pilot  must  open  control  loop  by  releasing  or  freezing 
the  stick. 

6 

Figure  A3.  Pilot-in-the-Loop  Oscillation  Rating  Scale 


A-4 


Appendix  B 


Mathematical  Development 


B.  1  Detailed  H2  Mathematical  Development 

In  mathematical  shorthand,  H2  optimization  can  be  written  as;  Determine  a  compensator  K(s) 
such  that 

^0  ~  inf  NI2  ~ 

K{s)Stabilizmg 

=  inf 

K{s)Stabilizlng 

~  inf  \Pzw  ~  Pyw 

K{s)Stabilizmg 

where  ttp  is  the  minimum  achievable  two-norm. 

A  state  space  realization  of  P  is  represented  by 

x  =  A2X  +  B^w  +  B^u 

z=-C^x  +  D^„w  +  Dji 

y  =  C2X  +  Dy^w  +  Dy^u 

The  “2”  subscript  indicates  that  this  is  H2  optimization  and  will  become  important  later  when 
mixed  H2  /  optimization  is  discussed.  The  following  assumptions  apply  to  H2  optimization: 

(i)  D2yv  =  0 

(ii)  Dyu  =  0 

(iii)  (A2,B2)  is  stabilizable  and  (C2,A2)  is  detectable 

(iv)  =  /  and  D^^Dj^  =  I 


(1) 

(2) 

(3) 


(4) 

(5) 

(6) 
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(V) 


has  full  column  rank  for  all  G) . 


(Vi) 


A^-jo)!  B^' 

^2  ^yw 


has  full  row  rank  for  all  co . 


Condition  (i)  is  imposed  to  ensure  the  two-norm  of  is  finite.  Condition  (ii)  is  included  to 
simplify  the  problem;  however,  it  is  not  formally  required.  Condition  (iii)  is  necessary  for  a 
stabilizing  solution  to  exist.  Condition  (iv)  is  required  for  regularity.  It  ensures  that  there  is  a 
direct  penalty  on  control  use  and  no  perfect  measurements.  This  condition  can  be  relaxed  to  a 
full  rank  condition  through  scaling  [9].  Finally,  conditions  (v)  and  (vi)  are  required  to  ensure 
existence  of  stabilizing  solutions  to  the  two  ARE's  used  in  solving  the  //j  problem. 


The  controller  that  minimizes  (1)  is  unique  and  will  be  denoted  K2opU  with  a  corresponding 
minimum  two-norm  .  It  is  desired  to  parameterize  the  family  of  stabilizing  sub-optimal  i/j 
controllers  for  the  purpose  of  trading  off  performance  for  performance.  This  family  of 

controllers  can  be  parameterized  through  an  LFT  of  a  transfer  function  Jand  a  constrained 
freedom  parameter  Q  BH2,  as  shown  in  Figure  B 1 .  A  particular  form  of  J is  expressed  as 
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Figure  B 1 .  H2  System  with  Parameterized  Controller 


where 


Aj  =  A,  -  -  B,K, 

(8) 

K,=BlX,+DlC, 

(9) 

(10) 

I 

II 

(11) 

Kj,=B, 

(12) 

where  X2  and  Y2  real,  unique,  positive  semidefinite  solution  to  the  following  AREs 
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{A,-B,Dl,CyX,+X,iA,-B,DlC^)-X,B,BlX,  +CjQ  =  0  (13) 


where 

C,=U-D„Dl)C,  (15) 

and 

iA,  -  B,D'^C,fY,  +  r,(A,  -  B,d;.C,)-  +  B^  =  0  (16) 

where 

B„=BM-D^D^)  (17) 

The  family  of  controllers  that  produce  ^  or  can  now  be  parameterized  by 

K(s)  =  F,[J(s),Q(s)]  (18) 

where  Q  can  be  chosen  to  be  any  Q  SH2  such  that 

Ml<a^-al  (19) 

The  optimal  H2  controller  is  obtained  from  the  above  parameterization  when  Q  is  chosen  equal 
to  zero.  The  resulting  optimal  controller  is  given  by 


K. 


lopt 


-K,  0 


(20) 


B-5 


B.2  Detailed  Ha,  Mathematical  Development 


where  Xo  is  the  minimum  achievable  infinity-norm.  Recall  that  the  infinity-norm  of  be 

calculated  by  finding  the  maximum  singular  value,  ct  .  The  mathematical  expression  is 

fcIL  =supa[7;^]  (24) 

(O 

A  state  space  realization  of  (2.13)  is  expressed  as 


JC  =  A^x  +  Bjd  +  B^u 

(25) 

d  =  C^x  +  D,jd  +  D^ji 

(26) 

y  =  C„x  +  D^d  +  D„„u 

(27) 

It  should  be  noted  that  Z)^,here  is  the  same  as  in  the  //j  problem,  if  the  plants  are  the 
same.  This  fact  will  be  very  important  later  when  the  mixed  H2  /  problem  is  discussed. 
The  CO  subscript  indicates  the  problem  is  setup  for  optimization.  As  with  H2  optimization, 
there  are  several  assumptions  that  are  made: 

(i)£>g^=0 


B-6 


(ii)  =  0 


(iii)  (A^,B^)  is  stabilizable  and  ( C„ ,  )  is  detectable 

(iv)  A»A«  =  ^  and  ^ 


(V) 


A  -jo)I 

a 


(vi) 


A 

A' 

A.. 


has  full  column  rank  for  all  (O 


has  full  row  rank  for  all  a 


Conditions  (i)  and  (ii)  are  included  to  simplify  the  problem;  however,  they  are  not  necessary  for 
a  solution  to  exist.  Condition  (iii)  is  necessary  for  a  stabilizing  solution  to  exist.  Condition  (iv) 
is  required  for  regularity.  It  ensures  that  there  is  a  direct  penalty  on  control  use  and  no  perfect 
measurements.  This  condition  can  be  relaxed  to  a  full  rank  condition  through  scaling[9]. 

Finally,  conditions  (v)  and  (vi)  are  required  to  ensure  existence  of  stabilizing  solutions  to  the  two 
ARE's  used  in  solving  the  problem. 


While  the  conditions  above  are  similar  to  H2  optimization,  there  are  major  differences 
associated  with  the  solution  to  the  optimization  problem.  In  general,  the  controller  that 
achieves  in  equation  (21)  is  not  unique.  Furthermore,  y^  is  found  through  an  iterative 
method  based  on  the  solution  to  two  AREs  and  a  coupling  condition.  This  method  is  based  on 
the  parameterization  of  all  suboptimal  controllers  where 


(28) 
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for  some  /  <  Xo  •  Note  that  (28)  excludes  /o  >  but  allows  it  to  be  approached  to  any  desired 
tolerance.  The  family  of  all  admissable  controllers  that  satisfy  (28)  is  illustrated  by  the  LFT 
shown  in  Figure  B2. 


Figure  B2.  H„  System  with  Parameterized  Controller 

Even  though  there  are  major  differences  in  the  solution  to  the  optimization  problem,  the 
optimal  controller,  ,  can  still  be  found  by  applying  equation  (7),  which  is  repeated  below 
for  convenience.  Of  course,  the  elements  in  (7)  are  defined  differently,  so  that 


where 


*^uy  *^ur 

*^vy  *^vr 

4/ 

Jis)  = 

= 

0 

I 

1 

0 
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(29) 


KfC„  -  +r-Xc: (C,  - 


K=<.B.x,+D:.c,xi-r'Xx.r' 

(30) 

k,=y,c1+b,dI, 

(31) 

K^  =  -(r-^D^B;x,*c.xi-r'Xx.r' 

(32) 

K,=r"Y,c:D„+B, 

(33) 

where  and  are  symmetric,  positive  semidefinite  solutions  to  the  following  AREs 


(A  -B  d'^C  /x  +X  (A  -B  d'^C  )+X  (g~^B,B.-B  b'^ )X  +C^C  =0  (34) 

'  m  rr\  o^tt  a '  rr>  rr\  '  rrs  pr\  an  a  ^  rr>  Cl  Cl  00  00  00  €  B 


00  oo  eu  e  00  oo'oo  oo  eu  e 


where 


q=(/-aX)Q 


(35) 


and 


(A^ -B,Dl,C^)Y^+YJA^ -B,Dl,Cy +  YJy-^CjC,-ClC^)Y^+B,Bl  =  0 


(36) 


where 


^di.^  ^yd^yd)  (37) 

Finally,  Q  can  be  chosen  to  be  any  Q  e  such  that 

\Ul<r  (38) 

Notice  that  Q  does  not  have  to  be  chosen  equal  to  zero  to  get  an  optimal  controller  here,  as  was 
the  case  in  H2  optimization. 
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BJ  Detailed  Mathematical  Development  for  the  Structured  Singular  Value. 

The  structured  singular  value  is  a  matrix  function  based  on  the  underlying  structure  of  a 
set  of  block  diagonal  matrices 

A.  =  {rffa45,/ . 5s4,A„...,A,]  |5,  €C,  A  eC''"')  (39) 

where  is  the  rth  scalar  block  of  order  r/  and  Aj  is  the yth  full  block  of  order  mj.  For 

simplicity,  we  will  assume  that  A  is  square,  but  the  theory  applies  for  no-square  perturbations  as 
well.  The  dimension  n  of  A  €  is  given  by 

s  F 

"  =  (40) 

;=1  >1 

We  will  want  to  place  a  limit  on  the  maximum  infinity-norm  that  A  can  have.  To  do  this,  we 
define  a  set  of  block  diagonal  matrices  that  have  an  infinity-norm  bounded  by  as 

BA:={AgA  \a(A)<y~'}  (41) 


The  structured  singular  value  of  a  complex  matrix  M  defined  over  a  set  of  perturbations  A  is 


min|cr(A)  |a  €  A,det(/  -  MA)  =  0} 


(42) 


unless  there  is  no  A  e  A  which  makes  I  —  A/A  singular,  in  which  case  (A/)  =  0 . 


From  the  definition  of  ^^(A/) ,  it  can  be  seen  that  the  maximum  singular  value  of  M  is  always 
an  upper  bound.  However,  this  bound  is  conservative  when  we  have  more  than  one  A  block. 
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One  method  of  reducing  this  conservativeness  of  the  upper  bound  is  to  consider  some 
transformation  that  does  not  affect  the  value  of  {M)  but  does  affect  the  value  of  .  We 
do  this  by  defining  a  set  of  scaling  transfer  functions  D  which  has  the  same  block  structure  as  the 
perturbations  A .  These  transfer  functions  are  given  by 


D:=  {[a . lA  =D;  >Q.dj  >o) 


Now,  an  upper  bound  is  given  by  the  following: 


Theorem  B.l  Assume  Me  is  defined  by  (39),  andH  is  defined  by  (43).  Then 

(M)  <  ^  a(DMD'^ )  (44) 

Proof:  See  [1 1],  Theorem  2.3.3. 

Therefore,  we  can  reduce  the  calculation  of  an  upper  bound  on  //  to  computing  the  maximum 
singular  value  of  a  matrix.  The  next  step  is  to  take  the  system  in  Figure  2.3  and  put  it  into  the  P- 
K  form  that  is  shown  in  Figure  2. 1 .  This  is  shown  in  Figure  B2. 

Next  we  define  the  Mtransfer  function  matrix  that  will  allow  us  to  calculate  // . 

M=FfiP,K)  (45) 
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Figure  B3.  P-K  Diagram  of  System  with  Uncertainty 


Using  equation  (45),  Figure  B2  reduces  to  the  M-  A  system  shown  below. 
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A  eA  = 


(47) 


^mil  ^  ^oul  >^P  ^  } 

The  main  result  of  fi  theory  can  now  be  stated  as  follows. 

Theorem  B.2  (Main  Loop  Theorem)  The  following  are  equivalent: 
i.  jUAQ^)<r 
«•  (^) 

A  (M’^om)]<r 

A„„,6BA„„ 

Proof:  [1 1],  Corollary  4.7. 

The  importance  of  this  theorem  is  that  we  now  have  a  single  test  on  Mthat  will  provide  us  with 

information  about  the  system  response  to  each  perturbation  block.  Let  ( M,, ^  <  y  be  a 

requirement  for  stability  and  ( M)  <y  be  a  desired  measure  of  performance.  Then  from 

the  Main  Loop  Theorem,  <  y  implies  our  performance  condition  is  satisfied  and  the 

system  is  robustly  stable.  The  theorem  also  allows  us  to  exchange  A^,  and  .  Now 

/r^(A/)  <  y  implies  that  our  system  has  robust  performance  for  all  perturbations.  The 

development  to  this  point  has  only  considered  the  case  where  Mis  a  constant  matrix.  Next  we 

expand  these  concepts  to  include  dynamic  matrices. 

B.3. 1  Frequency  Domain  //  -synthesis.  Let  M(s)  be  a  MIMO  transfer  function  with 

inputs  and  outputs.  Assume  A  cr  C"*'’"''  has  the  block  structure  given  in  (39).  Now  define 
the  set  of  all  dynamic  perturbations  which  have  the  desired  diagonal  structure  as 
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M(A):=  {A(^j  g5R//„  I  A(5(,)eA  forall  Sq  eC"^} 


(48) 


where  is  the  extended  closed  right-half  complex  plane. 

Now  we  can  define  the  complex  structured  singular  value  of  the  dynamic  transfer  matrix  M(s) 
over  the  structured  perturbations  A(s)  e  M  {A)  as 

\G(s)\  =supn[G(j(o)]  (49) 

Note  that  the  structured  singular  value  is  not  really  a  norm.  The  notation  above  is  used  only  for 
convenience.  Next  we  define  the  set  D  of  scaling  transfer  functions  which  have  the  same  block 
diagonal  structure  as  A ,  where  each  block  is  a  stable  perturbation.  Then  an  upper  bound  on  the 
structured  singular  value  of  the  transfer  matrix,  G{s),  is  given  by 


\G( <  sup  inf  ar( DMD  ‘ ) 

©691  DeD 

(50) 

=  inf||DMZ)-1 

(51) 

Equation  (51)  converts  the  p  problem  into  an  problem  that  includes  the  selection  of  the 
scaling  matrix  D.  We  can  now  perform  p  -synthesis  by  the  following  process: 

(1)  Determine  the  closed-loop  transfer  function  Af  by  performing  a  standard 
optimization  on  some  open-loop  plant,  P. 

(2)  Select  an  order  of  scaling  matrix  D,  that  best  fits  the  singular  values  of  A 
over  a  large  frequency  range  and  solve  the  optimization  problem  again. 

(3)  Iterate  on  (1)  and  (2)  until  the  value  of  p  stops  decreasing. 
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This  method  is  known  as  D-K Iteration  and  is  mechanized  in  MATLAB™  (see  [12]  for  details). 
This  method  is  not  guaranteed  to  converge  to  the  optimal  scaling  D  or  the  optimal  controller  K, 
but  in  practice  the  method  works  well.  The  D-K  iteration  method  results  in  a  controller  order 
equal  to  the  plant  order  plus  twice  the  order  of  the  D  scales.  This  implies  that  there  is  a  tradeoff 
between  controller  order  and  the  accuracy  of  the  D  scale  fits.  Also  of  note  is  that  since  //  - 
synthesis  is  based  on  designing  an  controller,  it  can  result  in  a  non-strictly  proper  controller 
(ie,  one  with  a  non-zero  D  matrix). 

B.4  Detailed  Mathematical  Development  for  the  Mixed  i/j  /  Problem. 

B.4. 1  State  Space  Formulation.  The  development  that  follows  was  taken  from  [4]. 

The  mixed  problem  is  shown  graphically  in  Figure  B5.  This  system  has  exogenous 

inputs  w  and  d.  The  controlled  outputs  of  the  system  are  z  and  e.  The  measured  output  is  y  and 
the  control  is  u.  These  definitions  are  the  same  as  those  in  Figures  B1  and  B2.  However,  P  is 
now  a  combination  of  the  P  matrices  from  the  and  subproblems. 


Figure  B5.  Block  Diagram  for  the  Mixed  H2  /  Problem 
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The  resulting  state  space  representation  of  P  is  given  by 

x  =  Ax  +  Bjd  +  B^w  +  B^u 
e  =  C^x  +  D^jd  + 
z  =  C^x  +  D^d  +  b^w  +  bji 

y  =  CyX  +  byjd  +  by„W-\-by„U 

The  transfer  function  form  is 

P 


(52) 

(53) 

(54) 


(56) 


We  can  also  represent  (56)  by  its  individual  //j  problems.  These  are  given  by 

equations  (4-6)  and  (25-27)  respectively.  The  assumptions  for  the  mixed  /fj  /  ^co  problem  are: 

(i)  Dyw  —  0 

(ii)  Dy^  =  0 

(iii)  {A2,Bii)  is  stabilizable  and  {Cy,A^  is  detectable 

(iv)  DlD^  =  I  and  =  / 

\A^-jcDl  B^2~\ 

(v)  has  full  column  rank  for  all  co . 


B-16 


has  full  row  rank  for  all  (O . 


(Vi) 


A^-jal 

^yl  ^yw 


Note  that  these  assumptions  are  the  same  ones  that  were  used  for  the  H2  problem.  There  are  no 
assumptions  included  dealing  with  the  part  of  the  mixed  problem  because  it  is  a  constraint. 

No  optimization  will  be  performed  and  as  such,  we  will  not  have  to  solve  any  Riccati 
equations.  We  can  also  allow  perfect  measurements  and  no  control  penalty  in  the  portion  of 
the  problem.  This  will  be  handled  by  the  H2  part  of  the  problem  since  the  controller  energy  will 
be  constrained  by  the  /fj  of  the  problem.  We  are  also  guaranteed  to  have  a  stabilizing 

controller  since  the  plant  is  the  same  for  both  problems.  A  state  space  representation  of  the 


controller  in  Figure  B5  is 

jc,  =  +  B^y  (57) 

u=C,x^  +  D^y  (58) 

Combining  (4-6)  and  (57-58)  produces  the  closed  loop  state  space  equation  for  ^zw 

Xj  =(^A  +  +  ^B^  +  B2D^DyJj  w  (59) 

K  =  (60) 

^  ~  (Az  ^zu^c^l)  ^2  ^zu^c^c  ^zu^c^yw^  (61) 


Notice  that  must  be  equal  to  zero  based  on  assumption  (i),  for  the  two-norm  of 

to  be  finite.  Also  notice  that  this  combined  with  assumption  (iv)  implies  that  Dq  =  0.  Thus,  the 
only  way  this  problem  has  a  solution  is  for  the  controller  K  to  be  strictly  proper.  If  we  close  the 
loop  on  our  mixed  H2  /  system  with  the  controller  represented  in  (59-61),  we  get 
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Xj  =  -^2  ^2  +  ^ 


(62) 


and 


where 


^2 


^2  -^2^ 


K 

BcDy. 


“  [bcD^,_ 

C=[C,  D^C,] 

Q=[Q  A„Q] 


(63) 

(64) 

(65) 

(66) 

(67) 

(68) 

(69) 

(70) 

(71) 

(72) 

(73) 
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(74) 


"^ed  ~  ^ed 

The  following  definitions  are  to  simplify  the  subsequent  development: 

=  the  unique  K(s)  that  makes  =« 

=  a  solution  to  the  /  //«  problem  for  some  y 

^-||TedLwhen^  =  ;i:,^^ 

/  ^\T,j\\^v^^^nK(s)=Kmix 

y  =  the  constraint  which  y*  must  stay  less  than  or  equal  to  (y  *  <  yj 

The  mixed  H2  /  H„  problem  can  now  be  restated  as  follows:  determine  a  K(s)  such  that 

(i)  The  underlying  Hj  and  problems  are  stable,  i.e.,  ^^e  stable 

(ii)  y*  <  y  for  some  given  y  >  y 

(iii)  ||r^„,||2  is  minimized. 

Walker  [4]  introduced  the  following  theorem  to  develop  this  problem: 

Theorem  B.4.1.1  Let  be  given  and  assume  there  exists  a  Q„  =  Ql  >0  satisfying 

e]  +  el  +  ^Pla^'''  +  =  0  (75) 
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where  R  =  {^^I  -  >  0 .  Then  the  following  are  equivalent: 

(i)  is stabilizable 

(ii)  is  stable 
(Hi)  is  stable. 

Moreover,  if  the  above  hold  then  the  following  are  true: 

(v)  the  two-norm  of  the  transfer  function  '^zw  is  given  by 

lTX=lr[c,Q,Cj]  =  lr[Q,CjC,] 

where  Q2  —  Q2  ^  0  w  the  solution  to  the  Lyapunov  equation 

^2  Q2  +  Q2  ^2  +  0 

(vi)  all  real  symmetric  solutions  of  Equation  (75)  are  positive  semidefinite 

(vii)  there  exists  a  unique  minimal  solution  to  Equation  (75)  in  the  class  of  real  symmetric 
solutions 

(viii)  is  the  minimal  solution  of  equation  (75)  if  and  only  if 

«-'<?,+ G.^J  rV,)]sO/or»ff/  (76) 

0^)  W.  <W  V  1#  A'V.)]  (<)  <0  y^here  Q,  is 

the  minimal  solution  to  Equation  (75). 
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The  main  result  from  this  theorem  is  that,  given  a  controller  that  is  closed-loop  stable  for 
the  H2  problem,  we  can  determine  the  minimum  level  of  the  constraint,  y ,  by  determining 

the  minimum  value  of  y  for  which  a  positive  semidefinite  solution  to  (75)  exists. 

Using  Theorem  B.4.1.1  the  mixed  problem  can  now  be  restated  as:  Determine  the  K(s) 
which  minimizes  the  objective  function 

44,B„C,)  =  (r[ac;C,]  (77) 

where  Q2  is  the  real,  symmetric,  positive  semidefinite  solution  to 

=  ^  (78) 

and  satisfying  the  constraint 

<  Q. + Q.  yil + (e.  (?,’■  +  e] + + sX  =  0  (79) 

where  is  the  real,  symmetric,  positive  semidefinite  solution. 

This  is  a  minimization  problem  with  two  equality  constraints  and  will  be  solved 
using  a  Lagrange  multiplier  approach.  An  excellent  discussion  of  this  method  can  be  found  in 
Arora  [13].  The  Lagrangian  for  the  mixed  problem  is 

L = »-[ac;c.] + fr{[4  a + Q^  4  +  bx]x] 

c:+bmr~'(q.  c:+ s.a'/ 

*BX]y\  m 

where  "X  and  ^  are  symmetric  Lagrange  multiplier  matrices. 

The  first  order  necessary  conditions  for  the  minimum  of  this  Lagrangian  are 


B-21 


(81) 


—  = 

dA„ 


=  Xl,Q,,+X,Q,+Y^,Q^,  +  Y,Q,=0 


^  +  Xl,V„  +  X,B^V,  +  Y^Q,Cl  +  Y,Ql,Cl, 

+  +  W,*{y^q.  +  r,e;)CA/+(i;[e^ +r,a)cJi)XAf = 0 


(82) 


ac 


“  ^u2^\Q\1  ^  ^ul^XlQl  ^llQu  ■*■  ^l^cQl  ■*■  ^uoJ^xQab  ^uoa^llQb 

+ + K^Q^ySQ^ + J!ia,i'2a + + ^*aai;je^ 

+ Kc.QLy,A + s.aai'ja + + ^a) + + j'^a) = »  (83) 


— -/^2a +a-^2^+^w^i  "0 


(84) 


.T 


dQ,  ■  ^ 


(85) 


dJi 

d% 


= ^„  e« + a  <+ ^; + + ^’X  =  o 


(86) 


a^ 

aa 


= R-^  ee+  a  el  R-^  elf  % 

+  ^(<+  ^■'  ^■'  ^e)  =  0 


(87) 


where 


M=R-^D,,Dl, 


P,  =  Z)I,M5J 


(88) 

(89) 

(90) 
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These  equations  have  not  been  solved  analytically  but  do  provide  some  insight  into  the  nature  of 
the  solution.  Equation  (86)  implies  either  ^=0  or  +  i?”’ + i?"' is 

neutrally  stable.  The  first  condition  means  that  the  solution  is  off  the  boundary  of  the 
constraint  (where  the  boundary  is  defined  as  the  constraint  being  satisfied  at  equality).  The 
second  condition  implies  that  the  solution  lies  on  the  boundary  of  the  constraint  and  is  a 
neutrally  stabilizing  solution  for  the  Riccati  equation  (78).  Since  Walker  was  not  able  to 
solve  the  problem  analytically,  a  numerical  approach  was  taken. 

Note  that  Equation  (78)  has  the  form 

=  0  (99) 

where  i?''  e,  +  6]  6^  ■ 

There  are  two  theorems  from  Snyders  and  Zakai  [15]  that  are  useful  at  this  point: 

Theorem  B.4.1.2  If  Ay  is  stable,  then  ^  =  0is  the  only  solution  to  Equation  (99). 

Recall  Theorem  B.4.1.1  (ix)  says  that  if  Ay  is  not  neutrally  stable,  is  not  the  minimal 
solution  to  (74).  Since  ^=0,  the  Lagrangian  (80)  reduces  to 

=  tr[Q,ClC,]  +  tr{[^,  Q,  +  Q,  (100) 

This  is  the  Lagrangian  associated  with  the  H2  problem  (note  we  still  have  to  solve  the  Lyapunov 
equation  in  the  constraint  to  evaluate  the  objective  function).  The  other  theorem  applicable  to 
(99)  is: 
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Theorem  B.4.1.3  Let  Ay  be  neutrally  stable.  Then  (79)  has  infinitely  many  ^  >  0  solutions  of 
possibly  varying  ranks. 

However;  theorem  B.4.1.1  (ix)  says  that  is  the  minimal  solution.  In  order  to  relate  this  to  the 
original  y  *  constraint,  Walker  states  the  following  theorem: 

Theorem  B.4.1.4  Assume  /i^is  stable  and  R  =  I  -  ^  If  there  exists  a  Q^>0 

satisfying 

then  the  following  are  equivalent: 

(^)  ll^e</|L=Y 

00  +  ^d^ld  ^.+  2co  ^e]  “  neutrally  stable 

Furthermore,  in  this  case  Q„  is  unique. 

Thus  we  can  discern  two  things  from  Equation  (79).  If  '^  =  0  we  have  an  unconstrained  H2 

problem.  If  Ay  =  +  SJDIj  i?”'  R~^  is  neutrally  stable,  we  are  on  the 

boundary  of  the  original  y  *  ^  y  constraint  and  is  the  neutrally  stabilizing  solution  to  (79). 

At  this  point.  Walker  fixed  the  order  of  the  controller  to  the  order  of  the  H2  subproblem  or 
greater.  This  led  to  the  following  theorem: 

Theorem  B.4.1.5  Assume  >  ^2 .  Then  the  following  hold: 

0)  If  y  <y ,  no  solution  to  the  mixed  problem  exists 
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(ii)  ^  Y  <  Y  <  y  ,  Kfnix  is  such  that  y*  =  Y 

(Hi)  If  y  >y ,  K2opt  is  the  solution  to  the  mixed  H2  /  problem. 

Proof:  See  [4].  ■ 

For  a  controller  with  order  greater  than  equal  to  the  order  of  the  H2  problem,  the  solution  to  the 
mixed  problem  with  y  <  y  <  y  lies  on  the  boundary  of  the  constraint,  y*  —y  . 

Thus,  in  this  region,  a*  is  a  montonically  decreasing  function  of  y  as  shown  in  Figure  B6.  This 
curve  is  called  a  Pareto  optimal  curve.  In  general,  points  on  a  Pareto  curve  represent  an  optimal 
tradeoff  between  an  objective  and  a  constraint.  Unfortunately,  since  the  solution  to  (75)  must  be 
the  neutrally  stabilizing  solution,  (87)  becomes  very  difficult  to  handle  as  we  move  to  the  left  on 
the  curve. 


Figure  B6.  Typical  Mixed  H2  /  y  vs.  a  Curve 
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Flight  Test  Matrix 


TABLE  Cl 


HAVE  INFINITY  II  PROFILE  TEST  MATRIX 


MISSION 

TEST 

CONFIG. 

TASK 

PRIQRITY 

CQMMENTS 

PROFILE 

POINT 

1* 

A 

HQDT-  1000’  Trail 

1 

Each  pilot  flies 

HQDT  -  Formation 

2 

Profile  1  first 

B 

HQDT  -  1000’  Trail 

1 

HQDT  -  Formation 

2 

C 

HQDT  -  1000’  Trail 

1 

*  Airborne  model 

HQDT  “  Formation 

2 

validation 

D 

HQDT  -  1000’  Trail 

1 

will  be  performed 

HQDT  -  Formation 

2 

by  the 

E 

HQDT  -  1000’  Trail 

1 

first  pilot  to  fly 

HQDT  -  Formation 

2 

profile  1 

F 

HQDT  -  1000’  Trail 

1 

HQDT  -  Formation 

2 

2 

1 

A 

Straight-in  Landing 

1 

Profile  Order 

2 

A 

Horizontal  Offset  Landing 

1 

Pilot  1-2,  3, 4 

3 

A 

Vertical  Offset  Landing 

2 

Pilot  2-3, 4, 2 

4 

B 

Straight-in  Landing 

1 

Pilot  3-4,  2,  3 

5 

B 

Horizontal  Offset  Landing 

1 

6 

B 

Vertical  Offset  Landing 

2 

7 

c 

Straight-in  Landing 

2 

8 

c 

Horizontal  Offset  Landing 

2 

9 

c 

Vertical  Offset  Landing 

3 

10 

D 

Straight-in  Landing 

2 

11 

D 

Horizontal  Offset  Landing 

2 

12 

D 

Vertical  Offset  Landing 

3 

E 

Straight-in  Landing 

3 

BH 

E 

Horizontal  Offset  Landing 

3 

■B 

E 

Vertical  Offset  Landing 

4 

F 

Straight-in  Landing 

3 

Bl 

F 

Horizontal  Offset  Landing 

3 

MBm 

F 

4 

Configurations  (In  evaluation  priority  order) 

A  -  Classic  Design  D  -  H2/H-Infinity  Pitch  Rate  Design 

B  -  H2  AOA  Design  E  -  H2  Pitch  Rate  Design 

C  -  H2/H-Infmity  AOA  Design  F  -  H2  AOA  Filter 


Note;  Pilots  will  not  be  informed  of  configurations  so  as  to  minimize  rating  biases. 
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PROFILE 


TEST  CONFIG. 
POINT 


PRIORITY 


COMMENTS 


Straight-in  Landing 
Horizontal  Offset  Landing 
Vertical  Offset  Landing 


Straight-in  Landing 
Horizontal  Offset  Landing 

Vertical  Offset  Landing _ 

Straight-in  Landing 
Horizontal  Offset  Landing 
Vertical  Offset  Landing 


Straight-in  Landing 
Horizontal  Offset  Landing 
Vertical  Offset  Landing 


Straight-in  Landing 
Horizontal  Offset  Landing 
Vertical  Offset  Landing 


Straight-in  Landing 
Horizontal  Offset  Landing 
Vertical  Offset  Landing 


Straight-in  Landing 
Horizontal  Offset  Landing 
Vertical  Offset  Landing 


Straight-in  Landing 
Horizontal  Offset  Landing 
Vertical  Offset  Landing 


Straight-in  Landing 
Horizontal  Offset  Landing 
Vertical  Offset  Landing 


Straight-in  Landing 
Horizontal  Offset  Landing 
Vertical  Offset  Landing 


Straight-in  Landing 
Horizontal  Offset  Landing 

Vertical  Offset  Landing _ 

Straight-in  Landing 
Horizontal  Offset  Landing 
Vertical  Offset  Landing 


Configurations  (In  evaluation  priority  order) 

A  -  Classic  Design  D  -  H2/H-Infinity  Pitch  Rate  Design 

B  -  H2  AOA  Design  E  -  H2  Pitch  Rate  Design 

C  -  H2/H-Infmity  AOA  Design  F  -  H2  AOA  Filter 


Note:  Pilots  will  not  be  informed  of  configurations  so  as  to  minimize  rating  biases. 
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Ground  and  Flight  Test  Verification  and  Validation  Data 


D.  1  Pole-Zero  Format 


The  pole-zero  notation  used  in  Chapter  FV  and  Tables  C1-C7  is  standard  Test  Pilot 
School  (TPS)  shorthand  for  expressing  the  poles  and  zeros  of  a  transfer  function.  This  notation 
lists  first  order  poles  and  zeros  as  positive  if  they  are  stable,  and  negative  if  they  are  unstable. 

For  example;  (2.0)  expresses  the  s-domain  monomial  s+2.  Setting  this  expression  equal  to  zero 
results  in  a  pole  of  s  =  -2.  Second  order  pole  and  zero  pairs  are  expressed  as  an  ideal  damping 
ratio  (Q  and  the  pole  or  zero  natural  frequency  (©n)  in  brackets.  For  example;  [0.5  4]  would 
express  the  s-domain  polynomial  s2+4s+16,  which  is  of  the  form  s2+2(^C0ii+«)ii^.  The  gain  is  the 
result  setting  s=0  in  the  numerator  and  denominator  of  the  transfer  function.  The  optimal  control 
transfer  functions  are  divided  into  the  three  input  channels  of  the  control  law;  command,  pitch 
rate  (q),  and  angle  of  attack  (a). 


TABLE  D1 

POLE-ZERO  VERIFICATION-AIRCRAFT 


Design  Model 

Implementation  Model 

Poles  (rad/sec) 

(0.3029) 

(0.3028) 

(-1.4971) 

(-1.4970) 

Zeros  (rad/sec) 

(-0.8329) 

(-0.8329) 

(-70.2309) 

(-70.2310) 

D-2 


TABLE  D2 


POLE-ZERO  VERIFICATION  -  H2INI/H2AIN 


Design  Model 

Implementation  Model 

Poles  (rad/sec) 

(6.0) 

(6.0) 

Command  Channel 

(0.8405) 

(0.8405) 

[0.9999,31.6083] 

[0.9999,31.6083] 

Zeros  (rad/sec) 

(21.4843) 

(21.4843) 

Command  Channel 

Gain 

Command  Channel 

-107.42 

-107.42 

Zeros  (rad/sec) 

(6.0) 

(6.0) 

Pitch  Rate  Channel 

(4.0809) 

(1.4402) 

(1.4402) 

Gain  -  q  Channel 

58.71 

58.71 

Zeros  (rad/sec) 

(6.0) 

-0.8329 

AOA  Channel 

(3.9061) 

-70.2310 

(1.5840) 

Gain  -  a  Channel 

0.50 

0.4999 

D-3 


TABLE  D3 


POLE-ZERO  VERIFICATION  -  MXAOA 


Design  Model 

Implementation  Model 

Poles  (rad/sec) 

Command  Channel 

(85.4428) 

[0.6456, 25.3192] 

[0.9375,  5.1344] 

[0.9454, 0.3811] 

(85.4428) 

[0.6456,25.3192] 

[0.9375,  5.1344] 

[0.9454, 0.3811] 

Zeros  (rad/sec) 

Command  Channel 

(53.0559) 

(1.3192) 

(0.8195) 

(0.4175) 

[0.7064, 24.8286] 

(53.0559) 

(1.3192) 

(0.8195) 

(0.4175) 

[0.7064, 24.8286] 

Gain ' 

Command  Channel 

-36.8823 

-36.8823 

Zeros  (rad/sec) 

Pitch  Rate  Channel 

(0.8936) 

(0.4391) 

[0.7898, 34.0126] 

[0.8267,  5.3258] 

(0.8936) 

(0.4391) 

[0.7898,34.0126] 

[0.8267, 5.3258] 

Gain  -  q  Channel 

43.4441 

43.4441 

Zeros  (rad/sec) 

AOA  Channel 

(60.5737) 

(4.9107) 

[0.6775, 24.5839] 

[0.9705, 0.4375] 

(60.5737) 

(4.9107) 

[0.6775, 24.5839] 

[0.9705, 0.4375] 

DC  Gain  -  a  Channel 

3.5357 

3.5357 

TABLE  D4 

POLE-ZERO  VERIFICATION-CLASSIC 


Design  Model 

Implementation  Model 

Gain  -  Kq 

Pitch  Rate  Feedback  Loop 

0,441 

0.441 

Gain  -  Ka 

AoA  Feedback  Loop 

1.60 

1.60 

NOTES  -  There  was  no  dynamic  compensation  used  in  this  design.  Only  the  feedback  gains  shown  above 
were  used  to  produce  the  Classic  design. 
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TABLE  D5 


POLE-ZERO  VERIFICATION  -  H2AOA 


Design  Model 

Implementation  Model 

Poles  (rad/sec) 

All  Channels 

Zeros  (rad/sec) 

(6.6167) 

Command  Channel 

(1.5171) 

(0.2990) 

(0.2990) 

Gain 

Command  Channel 

-14.5252 

-14.5252 

Zeros  (rad/sec) 

(1.5461) 

(1.5461) 

Pitch  Rate  Channel 

[0.53, 2.30] 

[0.53,2.30] 

Gain  -  q  Channel 

11.0656 

11.0656 

Zeros  (rad/sec) 

(1.4962) 

(1.4962) 

AOA  Channel 

[0.53, 2.30] 

[0.53, 2.30] 

Gain  -  a  Channel 

0.1054 

0.1054 
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TABLE  D6 


POLE-ZERO  VERIFICATION  -  MXINI 


Design  Model 

Implementation  Model 

Poles  (rad/sec) 

Command  Channel 

(1.9051) 

(84.9036) 

[0.6906, 29.4591] 

[0.1563,  1.3921] 

(1.9051) 

(84.9036) 

[0.6906, 29.4591] 

[0.1563,  1.3921] 

Zeros  (rad/sec) 

Command  Channel 

(1.8081) 

(0.7184) 

(216.1494) 

[0.3888, 28.1514] 

(1.8081) 

(0.7184) 

(216.1494) 

[0.3888,28.1514] 

Gain  . 

Command  Channel 

-0.5776 

-0.5776 

Zeros  (rad/sec) 

Pitch  Rate  Channel 

(2.4557) 

[0.8242,  55.7059] 

[0.6217, 1.1901] 

(2.4557) 

[0.8242,  55.7059] 

[0.6217, 1.1901] 

Gain  -  q  Channel 

15.6905 

15.6905 

Zeros  (rad/sec) 

AoA  Channel 

(6.3463) 

[0.2511,  1.6603] 

[0.6896, 44.3650] 

[0.2511, 1.6603] 

[0.6896, 44.3650] 

Poles  (rad/sec) 
a  poles 

(1.9051) 

(84.9036) 

[0.6906, 29.4591] 

[0.1563, 1.3921] 

(0.8015) 

(2.4179) 

(1.0) 

[0.6905, 29.4951] 

Gain  -  a  Channel 

1.5090 

1.5090 
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Test  Aircraft:  Lear24,  N101VS 
Configuration:  Ground  Simulation 
Flight  Control  Config:  CLASSIC 


Task:  PTI  Pulse 
(1/2  inch  defi,  3  sec  dur.) 
Date:  01  October  97 
Source:  Ground  Test 


Figure  Dl.  Ground  Test  Closed-Loop  q  and  a  Time  Responses  -  CLASSIC 


Test  Aircraft:  Lear  24,  N1 01  VS  Task:  PTI  Pulse 

Configuration:  Ground  Simulation  (1/2  inch  defI,  3  sec  dur.) 

Flight  Control  Config:  H2AOA  Date:  01  October  97 

Source:  Ground  Test 


Figure  D2.  Ground  Test  Closed-Loop  q  and  a  Time  Responses  -  H2AOA 
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Amplitude  q-deg/s,  alpha-deg,  pti-inches  Amplitude  q-deg/s,  alpha-deg,  pti-inches 


Test  Aircraft:  Lear  24,  N1 01  VS 

Task:  PTI  Pulse 

Configuration:  Ground  Simulation 

(1/2  inch  defI,  3  sec  dur.) 

Flight  Control  Config:  H2INI 

Date:  01  October  97 

Source:  Ground  Test 

Time,  seconds 

Figure  D3.  Ground  Test  Closed-Loop  q  and  a  Time  Responses  -  H2INI 


Test  Aircraft:  Lear24,  N101VS  Task:  PTI  Pulse 

Configuration:  Ground  Simulation  (1/2  inch  defi,  3  sec  dur.) 

Flight  Control  Config:  H2AIN  Date:  01  October  97 

Source:  Ground  Test 


Figure  D4.  Ground  Test  Closed-Loop  q  and  a  Time  Responses  -  H2AIN 
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Amplitude  q-deg/s,  alpha-deg,  pti-inches 


4  5  6  7  8  9  1G  11 

Time,  seconds 

Figure  D6.  Flight  Test  Closed-Loop  q  and  a  Time  Responses  -  CLASSIC 
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Amplitude  q-deg/s,  alpha-deg,  pti-inches  Amplitude  q-deg/s,  alpha-deg,  pti-inches 


Test  Aircraft:  Lear24,  N101VS 

Task:  PTI  Pulse 

Configuration:  140  KIAS,  8K  MSL 

(1  inch,  3  secdur.) 

Gear  Down,  Flaps  20® 

Date:  02  October  97 

Flight  Control  Config:  H2AOA 

Source:  Flight  Test 

’igure  D7.  Flight  Test  Closed-Loop  q  and  a  Time  Responses  -  H2AOA 


Test  Aircraft:  Lear  24,  N101VS 

Task:  PTI  Pulse 

Configuration:  140  KIAS,  8K  MSL 

(1  inch,  3  sec  dur.) 

Gear  Down,  Flaps  20® 

Date:  02  October  97 

Flight  Control  Config:  H2INI 

Source:  Flight  Test 
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Test  Aircraft:  Lear  24,  N1 01  VS 

Task:  PTI  Pulse 

Configuration:  140  KIAS,  8K  MSL 

(1  inch,  3  sec  dur.) 

Gear  Down,  Flaps  20° 

Date:  02  October  97 

Flight  Control  Config:  H2AIN 

Source:  Flight  Test 

Time,  seconds 

Figure  D9.  Flight  Test  Closed-Loop  q  and  a  Time  Responses  -  H2AIN 
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TABLED? 


CLOSED  LOOP  POLE  ANALYSIS  -  CLASSIC 


Closed  Loop  System  Poles 
without  Phugoid  Mode 

Closed  Loop  System  Poles 
with  Phugoid  Mode 

[0.7001,59.1105] 

[0.5322, 2.2772] 

[0.7002,  59.1573] 

[0.5192, 2.2702] 

[0.1424, 0.1698] 

TABLE  D8 

CLOSED  LOOP  POLE  ANALYSIS  -  H2AOA 


Closed  Loop  System  Poles 
without  Phugoid  Mode 

Closed  Loop  System  Poles 
with  Phugoid  Mode 

[0.7038, 60.2400] 
[0.6921,7.4210] 

(1.5176) 

[0.5300, 2.3000] 

(0.3013) 

[0.7038,  60.2399] 

[0.6921,  7.4256] 

(1.5127) 

[0.5300, 2.3000] 

[0.7308, 0.3161] 

(-0.1123) 

TABLE  D9 

CLOSEDLOOP  POLE  ANALYSIS  -  H2AIN 


Closed  Loop  System  Poles 
without  Phugoid  Mode 

Closed  Loop  System  Poles 
with  Phugoid  Mode 

[0.72, 277.41] 

[0.72, 277.41] 

[0.70, 60.00] 

[0.70,  60.00] 

[0.92, 33.43] 

[0.92,  33.43] 

(6.00) 

(6.00) 

(1.52) 

(1.51) 

(0.84) 

(0.88) 

(0.30) 

[0.66,  0.35] 

(-0.14) 
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TABLE  DIO 


CLOSED  LOOP  POLE  ANALYSIS  -  H2INI 


Closed  Loop  System  Poles 
without  Phugoid  Mode 

Closed  Loop  System  Poles 
with  Phugoid  Mode 

[0.72, 277.41] 

[0.72, 277.41] 

[0.70, 60.00] 

[0.70, 60.00] 

[0.92,  33.43] 

[0.92,  33.43] 

(6.00) 

(6.00) 

(1.52) 

(1.51) 

(0.84) 

(0.88) 

(0.30) 

[0.66, 0.35] 

(-0.14) 

TABLE  Dll 

CLOSED  LOOP  POLE  ANALYSIS  -  MXAOA 


Closed  Loop  System  Poles 
without  Phugoid  Mode 

Closed  Loop  System  Poles 
with  Phugoid  Mode 

(84.3732) 

[0.7115,59.5679] 

[0.6186,24.3357] 

(6.8060) 

[0.8230, 2.6992] 

(2.0367) 

(0.7097) 

(0.4358) 

(84.3735) 

[0.7115, 59.5672] 
[0.6186,24.3363] 

(6.8240) 

[0.8247, 2.6746] 

(2.0477) 

(0.6733) 

(0.4345) 

[0.6605,  0.0697] 
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Appendix  E 


Handling  Qualities  Rating  Data 


Number  of  Occurrences  Number  of  Occurrences 


Test  Aircraft:  Lear  24,  N1 01  VS  Task:  Straight  In  Landing 

Configuration:  Gear  Down,  Flaps  20%,  FCS  CLASSIC  Dates:  6  tolO  October  1997 

Altitude/Airspeed:  4,300  ft  MSL/1 30  KIAS  Data  Source:  Flight  Test 


Cooper-Harper  Rating 

Figure  El.  Straight-In  Landing  Cooper-Harper  Ratings  —  CLASSIC 


Test  Aircraft:  Lear  24,  N1 01  vs  Task:  Straight  In  Landing 

Configuration:  Gear  Down,  Flaps  20%,  FCS  CLASSIC  Dates:  6  tolO  October  1997 

Altitude/Airspeed:  4,300  ft  MSL/1 30  KIAS _  Data  Source:  Flight  Test _ _ 


_ 

■  Pilot  1 
□  Pilot  2 

■  Pilots 

■ 

■ 

1  2  3  4  5  6 

PIO  Rating 

Figure  E2.  Straight-In  Landing  PIO  Ratings  --  CLASSIC 
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Number  of  Occurrences  Number  of  Occurrences 


Figure  E3.  Horizontal  Offset  Landing  Cooper-Harper  Ratings  ~  CLASSIC 


Test  Aircraft:  Lear  24,  N101  vs  Task:  Horizontal  Offset  Landing 

Configuration:  Gear  Down,  Flaps  20%,  FCS  CLASSIC  Dates:  6  to10  October  1997 

Altitude/Airspeed:  4,300ft  MSL/1 30  KIAS  Data  Source:  Flightiest 


I  ■  Pilot  1 

j  □  Pilot  2 

I  H  Pilot  3 

1 

— 

m 

.1 _ 1 _ ^ _ 1 _ 1 _ 

1  2  3  4  5  6 


PIO  Rating 


Figure  E4.  Horizontal  Offset  Landing  PIO  Ratings  —  CLASSIC 
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Number  of  Occurrences  Number  of  Occurrences 


Figure  E7.  Horizontal  Offset  Landing  Cooper-Harper  Ratings  —  H2AOA 


7 

6 

5 

4 

3 

2 

1 

0 

1  2  3  4  5  6 

PIO  Rating 

Figure  E8.  Horizontal  Offset  Landing  PIO  Ratings  —  H2AOA 


Test  Aircraft:  Lear  24,  N101VS  Task:  Horizontal  Offset  Landing 

Configuration:  Gear  Down,  Flaps  20%,  FCS  H2AOA  Dates:  6  tolO  October  1997 

Altitude/Airspeed:  4,300  ft  MSL/130  KIAS  Data  Source:  Flight  Test 


- 

■  Pilot  1 

□  Pilot  2 
■  Pilots  i 
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Number  of  Occurrences  Number  of  Occurrences 


Test  Aircraft:  Lear  24,  Nl’OIVS  Task:  Straight  In  Landing 

Configuration:  Gear  Down,  Flaps  20%,  FCS  H2INI  Dates:  6  tolO  October  1997 


1  2  3  4  5  6  7  8  9  10 

Cooper-Harper  Rating 

Figure  E9.  Straight-In  Landing  Cooper-Harper  Ratings  —  H2INI 


PIO  Rating 

Figure  ElO:  Straight-In  Landing  PIO  Ratings  —  H2INI 
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Number  of  Occurrences  Number  of  Occurrences 


1  2  3  4  5  6  7  8  9  10 

Cooper-Harper  Rating 


Figure  El  1 .  Horizontal  Offset  Landing  Cooper-Harper  Ratings  —  H2INI 


P!0  Rating 

Figure  E12.  Horizontal  Offset  Landing  PIO  Ratings  -  H2INI 
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Number  of  Occurrences  Number  of  Occurrences 


Cooper-Harper  Rating 

Figure  El 3.  Straight-In  Landing  Cooper- Harper  Ratings  —  H2AIN 


Test  Aircraft:  Lear  24,  N101VS  Task:  Straight  In  Landing 

Configuration:  Gear  Down,  Flaps  20%,  FCS  H2AIN  Dates:  6  to10  October  1997 

Altitude/Airspeed:  8,000  ft  MSL/140  KIAS  Data  Source:  Flight  Test 


■  Pilot  1  I 

□  Pilot  2 

S  Pilot  3 

■ 

1 

.  M  . _ . 

1  2  3  4  5  6 


PIO  Rating 

Figure  El 4.  Straight-In  Landing  PIO  Ratings  -  H2AIN 
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Number  of  Occurrences  Number  of  Occurrences 


1  2  3  4  5  6  7  8  9  10 


Cooper-Harper  Rating 

Figure  E  15.  Horizontal  Offset  Landing  Cooper-Harper  Ratings  —  H2AIN 


Test  Aircraft:  Lear  24,  N101VS 

Configuration:  Gear  Down.  Flaps  20%,  FCS  H2A1N 

lAltitude/Airspeed;  8,000  ft  MSL/140  KIAS _ 


Task:  Horizontal  Offset  Landing 
Dates:  6  to1 0  October  1 997 
Data  Source:  Flight  Test  _ _ 


■  Pilot  1  i 

■ 

□  Pilot  2  i 
■  Pilots  ; 

1 

■ 

.--.---^1111111-, .  1 

H 

3  4 

PIO  Rating 


Figure  El 6.  Horizontal  Offset  Landing  PIO  Ratings  —  H2AIN 
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Number  of  Occurrences  Number  of  Occurrences 


Cooper-Harper  Rating 

Figure  El 7.  Straight-In  Landing  Cooper- Harper  Ratings  -  MXAOA 


Task:  Straight  In  Landing 
Dates:  6  to10  October  1997 
Data  Source:  Flight  Test 


■  Pilot  1 
□  Pilot  2 

■  Pilot  3 


1  2  3  4  5  6 

PIO  Rating 

Figure  El  8.  Straight-In  Landing  PIO  Ratings  —  MXAOA 
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Number  of  Occurrences  Number  of  Occurrences 


Test  Aircraft:  Lear  24,  N1 01  VS  Task:  Horizontal  Offset 

Configuration:  Gear  Down,  Flaps  20%,  FCS  MXAOA  Dates:  6  tolO  October  1997 

Altitude/Airspeed:  4,300  ft  MSL/130  KIAS  Data  Source:  Flight  Test 


Cooper-Harper  Rating 

Figure  El 9.  Horizontal  Offset  Landing  Cooper-Harper  Ratings  --  MXAOA 


PIO  Rating 

Figure  E20.  Horizontal  Offset  Landing  PIO  Ratings  -  MXAOA 


Flight#!  Mission  date:  6  Oct  97  Eval  Pilot:  Pilot  #l(Boe) 

Winds:  240/10 


Appr 

Task 

Cotin^u  ration 

Laiulin^ 

Zone 

C-M  ,P1()  (.."ooj  men  IS 

1 

St-In 

H2INI 

Desired 

4 

1  minor  pitch 

deficiencies/heavy  nose  after 
touchdown 

2 

St-In 

CLASSIC 

Desired 

3 

1  could  put  the  aircraft  where 

the  pilot  wanted 

3 

St-In 

MXAOA/Hoo  AOA  Cmd 

Desired 

2 

1  easy  to  set  pitch  picture/very 

good  flying  qualities 

4 

Horiz 

H2INI 

Desired 

5 

2  stick  force  change  for  a  given 

displacement 

5 

Horiz 

CLASSIC 

.  Desired 

no  grade  -  pilot  unsure  of 
workload 

6 

Horiz 

CLASSIC 

Desired 

3 

1  good  flying  qualities 

7 

Horiz 

MXAOA/Hoo  AOA  Cmd 

Desired 

2 

1  best  of  all  3  flight  control 

systems  tested  today/ 
negligible  deficiencies 

Flight  #2 

Mission  date:  6  Oct  97 

Winds:  240/1 5G25 

Eval  Pilot:  Pilot  #3  (Cantiello) 

Appr 

lask 

<  'onli^nration 

Landing 

Zone 

OH  Pk 

()  ('omnuMiTs 

1 

St-In 

H2AIN 

Desired 

4  2 

trimming  continuously  /  pitch 
sensitivity  in  the  flare/ 
nose  heavy  after  touchdown 

2 

St-In 

H2AOA 

Desired 

5  3 

hard  to  trim/very  sensitive  in 
pitch/heavy  nose  after 
touchdown 

3 

St-In 

MXAOA 

Desired 

4  3 

pitch  “bobble”  in  the  flare 

4 

St-In 

H2INI 

Desired 

2  1 

not  affected  by  gusts/ 
very  good  flying  qualities 

5 

Horiz 

H2AIN 

Desired 

4  2 

higher  workload  due  to  gusty 
winds 

6 

Horiz 

H2AOA 

Desired 

4  3 

trimming  continuously/ 
pitch  sensitivity  in  the  flare/ 
nose  heavy  after  touchdown 

7 

Horiz 

MXAOA 

Adequate 

5  3 

pitch  sensitive  in  the  flare/ 
light  turbulence 

8 

Horiz 

H2INI 

Desired 

4  2 

mild  undesirable  motion  and 
pitch  sensitivity  in  the  flare/ 
nose  heavy  after  touchdown 
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Flight  #3  Mission  date:  7  Oct  97  Eval  Pilot:  Pilot  #2  (Stevenson) 

Winds:  280/18 


Appr 

Tiisk 

Laruiiriij; 

Zone 

C-H 

PH) 

ConinR'ois 

1 

St-In 

MXAOA 

Desired 

4 

2 

tendency  to  overshoot  desired 
pitch  attitude 

2 

St-In 

H2AOA 

Desired 

6 

4 

small  amplitude  pilot-in-the- 
loop  oscillation  (PIO)  in  flare/ 
heavy  nose  on  landing 

3 

St-In 

H2AIN 

Adequate 

5 

2 

sluggish  initial  response 

4 

Horiz 

MXAOA 

Desired 

3 

2 

predictable/no  speed  stability 
feedback/easy  to  fly 

5 

Horiz 

H2AOA 

Desired 

4 

1 

divergent  when  off  trim 
airspeed 

6 

Horiz 

H2AIN 

Desired 

4 

2 

nose  pitch  up  below  trim 
airspeed/higher  workload  due 
to  constant  trim 

7 

Horiz 

H2AOA 

Desired 

6 

4  . 

heavy  stick  in  flare  led  to  PIO/ 
sensitive  to  pilot  bandwidth 

Flight  #4 

Mission  date:  7  Oct  97 

Winds:  270/20 

Eval  Pilot:  Pilot  #2  (Stevenson) 

Appr 

Task 

Conllgu  nUion 

Liiniiing 

Zone 

OH 

P!<) 

( T)in  ments 

1 

St-In 

CLASSIC 

Desired 

3 

1 

predictable  and  responsive 

2 

St-In 

H2INI 

Desired 

3 

1 

predictable/rapid  initial  pitch 
response 

3 

St-In 

MXAOA 

Desired 

3 

1 

linear  response/slight  pitch 
overshoot  in  the  flare 

4 

Horiz 

CLASSIC 

Desired 

3 

1 

rapid  initial  response 

5 

Horiz 

H2INI 

Desired 

4 

2 

sensitive  to  airspeed/high 
workload  to  maintain 
airspeed 

6 

Horiz 

MXAOA 

Desired 

3 

2 

small  pitch  overshoot  in  flare 
(turbulence) 

7 

Horiz 

H2INI 

Desired 

4 

3 

airspeed  pitch  sensitivity 
required  higher  workload 

8 

Horiz 

MXAOA 

Desired 

3 

1 

precise  pitch  attitude  changes 
required  lower  gain  pilot 
technique 
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Flight  #5  Mission  date:  8  Oct  97  Eval  Pilot:  Pilot  #3  (Cantiello) 

Winds:  Calm 


Appr 

Task 

Configuration 

laiiHlind 

Zone 

C-H 

P!() 

Com  mciits 

1 

St-In 

CLASSIC 

Desired 

2 

1 

.  beautiful,  very  good  flying 
qualities 

2 

St-In 

H2INI 

Desired 

2 

1 

•  very  good  flying  qualities/ 
easy  to  trim 

3 

St-In 

MXAOA 

Desired 

4 

3 

more  sensitive  in  the 
flair/tendency  to  float 

4 

St-In 

H2AIN 

Desired 

4 

3 

pitch  sensitive  in  the  flair/had 
to  stay  low  gain 

5 

St-In 

H2AOA 

Desired 

3 

2 

heavy  nose  in  the  flare 

6 

Horiz 

H2INI 

Desired 

3 

1 

'  very  nice  to  fly/n ice  control 
harmony  and  response 

7 

Horiz 

CLASSIC 

Desired 

4 

3 

moderate  \vorkload  due  to 
sensitivity  in  pitch  and  mild 
undesirable  motion 

8 

Horiz 

MXAOA 

Desired 

4 

3 

mild  undesirable  motion/ 
tendency  to  float  in  flare 

9 

Horiz 

H2AIN 

Desired 

3 

2 

slightly  sensitive  in  the  flare 

Flight  #6 

Mission  date: 
Winds:  Calm 

8  Oct  97 

Eval  Pilot:  Pilot  #1  (Boe) 

Appr 

Task 

Con  Induration 

1  .anding 
Zone 

C-H 

PIG 

Com  ments 

1 

St-In 

H2AIN 

Desired 

3 

1 

no  trim  problems 

2 

St-In 

H2AOA 

Desired 

6 

3 

divergent/worst  flown 
configuration/  heavy  stick 
forces  for  small  speed  changes 

3 

St-In 

CLASSIC 

Desired 

2 

1  . 

very  good  flying  qualities 

4 

St-In 

MXAOA 

Desired 

3 

2 

pitch  “bobble”  in  flare 

5 

Horiz 

H2AOA 

Desired 

4 

2 

hard  to  trim/heavy  nose  after 
touchdown 

6 

Horiz 

H2AIN 

Desired 

4 

2 

continuously  trimming/  heavy 
nose  in  flare 

7 

Horiz 

MXAOA 

Desired 

3 

1 

no  trim  problems 

8 

Horiz 

CLASSIC 

Desired 

** 

no  grade  -  pilot  unsure  of 
workload 

9 

Horiz 

CLASSIC 

Desired 

3 

1 

good  flying  qualities/light 
turbulence 
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Flight  #9  Mission  date:  9  Oct  97  Eval  Pilot:  Pilot  #3  (Cantiello) 

Winds:  230/10 


Appr 

Tusk 

Con  n<^u  ratio 

l.aii(iin<4 

Zone 

C-H 

PK) 

ComincHis 

1 

St-In 

H2INI 

Desired 

2 

2 

very  good  fly  qualities 

2 

St-In 

MXAOA 

Desired 

4 

1 

pitch  sensitivity  in  the  flare/ 
light  stick  forces/tendency  to 
float 

3 

St-In 

MXAOA 

Desired 

4 

3 

undesirable  motions  which 
compromised  the  task 

4 

St-In 

H2AIN 

Desired 

3, 

2 

hard  to  trim  precisely/small 
pitch  “bobble”  in  flare 

5 

Horiz 

H2INI 

Desired 

3 

1 

good  control  harmony 

6 

Horiz 

MXAOA 

Desired 

2 

2 

no  trim  problems/easy  to  fly 

7 

Horiz 

H2AIN 

Desired 

3 

2,  ' 

slight  trim  compensation 
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Appendix  F 


Handling  Qualities  Prediction  Data 


Figure  FI .  CLASSIC  Configuration  Time  Histories  with  Phugoid  Mode-Noise  Off 


F-2 


Input  Command  Pulse 


Pitch  Rate,  5  deg,  4  sec  pulse 


time-seconds 

Elevator  Deflection,  5  deg,  4  sec  pulse 


time-seconds 

Angle  of  attack,  4  deg,  5  sec  pulse 


time-seconds 

Pitch  Angle,  5  deg,  4  sec  pulse 


/ 

'  T - 

.../j . 

i 

time-seconds 

Flight  Path  Angle,  5  deg,  4  sec  pulse 


1 

1 

1 

1 

1 

1  1 

r 

-J- _ i _ , 

S  10 

time-seconds  time-seconds 

Figure  F2.  H2A0A  Configuration  Time  Histories  with  Phugoid  Mode-Noise  off 


deg/sec 


Input  Command  Pulse 


Pitch  Rate,  5  deg,  4  sec  pulse 


rau  p-sec 


Figure  F7.  Hoh’s  Bandwidth  Handling  Qualities  Predictions  with  Phugoid 
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